WL-TR-94-2036 


AD-A283  952 


ELECTRON  INTERACTIONS  WITH 
NON-LINEAR  POLYATOMIC  MOLECULES 
AND  THEIR  RADICALS  AND  IONS 


Ashok  Kumar  Jain 
Charles  A.  Weatherford 


Florida  A  &:  M  University 
Physics  Department 
Tallahassee,  Florida  32307 


DECEMBER  1993 

FINAL  REPORT 


FOR  06/07/90  —  06/06/93 


DTIC 

ELECTE 

?SEP.0H994' 

'8 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


0^*94-28330 

IIII11BIHI 

AERO  PROPULSION  AND  POWER  DIRECTORATE 
WRIGHT  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT  PATTERSON  AFB  OH  45433-7251  1 


94  8  31 


149 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  Government-related  procurement,  the  United 
States  Government  incurs  no  responsibility  or  any  obligation  whatsoever.  The  fact  that 
the  government  may  have  formulated  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication,  or  otherwise  in  any 
manner  construed,  as  licensing  the  holder,  or  any  other  person  or  corporation;  or  as 
conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 

This  report  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


Research  Physicist 

Advanced  Plasma  Research  Section 

Power  Components  Branch 

Aerospace  Power  Division 

Aero  Propulsion  and  Power  Directorate 


FERRELL  M. TURNER 

Chief,  Advanced  Plasma  Research  Section 

Power  Components  Branch 

Aerospace  Power  Division 

Aero  Propulsion  and  Power  Directorate 


MICHAEL  D.  BRAYD1QH,  Lt  Col,  USAF 


Deputy  Chief 

Aerospace  Power  Division 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  notify  WL/POOC-3, 
WPAFB,  OH  45433-7919  to  help  us  maintain  a  current  mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OM8  No.  0704-0188 


Public  reporting  burden  to r  this  collection  of  information  i$  estimated  to  average  ?  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services.  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway.  Suite  1204.  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

December  1993 


3.  REPORT  TYPE  AND  DATES  COVERED 

FINAL  06/07/90—06/06/93 


4.  TITLE  AND  SUBTITLE 

ELECTRON  INTERACTIONS  WITH  NON-LINEAR 
POLYATOMIC  MOLECULES  AND  THEIR  RADICALS 
AND  IONS 


6.  AUTHOR(S) 

I  ASHOK  KUMAR  JAIN 

CHARLES  A.  WEATHERFORD 


5.  FUNDING  NUMBERS 

C  F336 1 5-90-C-2032 
PE  62203 
PR  2003 
TA  12 
WU  03 


i  7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

i 

3  PHYSICS  DEPARTMENT 
!  FLORIDA  A&M  UNIVERSITY 
TALLAHASSEE,  FL  32307 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

AERO  PROPULSION  AND  POWER  DIRECTORATE 
WRIGHT  LABORATORY 
AIR  FORCE  MATERIAL  COMMAND 
WRIGHT  PATTERSON  AFB,  OH  45433-7251 


|  12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

!  APPROVED  FOR  PUBLIC  RELEASE; 
3  DISTRIBUTION  IS  UNLIMITED 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 

WL-TR-94-2036 


12  b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

rhis  is  the  final  report  for  work  accomplished  during  the  period  of  June  1990 
ugh  1993  under  Air  Force  Contract  F33615-90-C-2032. 

rhe  work  involved  the  computation  of  collision  cross  sections  for  electron 
:ollisions  with  non-linear  polyatomic  molecules  and  their  radicals  and  ions.  A 
reliable  set  of  cross  section  data  for  elastic  and  inelastic  processes  in  these  systems 
tias  been  produced  without  resorting  to  a  fitting  procedure.  The  computational 
model  employed  the  adiabatic-nuclei-approximation  for  electron  scattering  from 
molecules  in  the  gas  phase. 

[n  addition,  a  method  for  calculating  the  required  cross  sections  by  solving  the 
chrodinger  equation  as  a  partial  differential  equation  using  a  finite  difference 
iscretization  has  been  developed. 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 


Standard  form  298  (Rev  2-89) 

Prescribed  by  ANSI  Std  239-18 
298-102 


TABLE  OF  CONTENTS 


Section  Page  No. 

LIST  OF  FIGURES  . . v 

PREFACE  AND  ACKNOWLEDGEME.  LS  . vii 

1.  INTRODUCTION . 1 

2.  SUMMARY  OF  WORK  PERFORMED  (June  1990  -  1993) . 4 

3.  LIST  OF  PUBLICATIONS  (June  1990  -  1993)  7 

4.  REVIEW  OF  THEORETICAL  APPROACH . 9 

4.1  Single-Center-Expansion  (SCE)  Iterative  Method . 9 

4.2  Static  and  Exchange  Interactions . 11 

4.3  Polarization  Interaction . 13 

4.4  Cross-Section  Formulae . 15 

5.  REVIEW  OF  NUMERICAL  DETAILS . 19 

5.1  Target  Wavefunctions . 

5.2  SCE  Parameters  from  Multicenter  Target  Wavefunctions  .  . 

5.3  An  Optimized  Iterative  Scheme . 

5.4  Computer  Codes . 


□  □ 


6.  RESULTS  AND  DISCUSSION . 29 

6.1  Electron-CH4  Cross  Sections . 29 

6.2  Electron-SiH4  Cross  Sections . 30 

6.3  Electron-GeH4  Cross  Sections . 35 

6.4  Electron-NH3  Cross  Sections . 37 

6.5  Electron-H20  Cross  Sections . 39 

6.6  EIectron-H2S  Cross  Sections . 40 

7.  ELECTRON  SCATTERING  FROM  C2H2  MOLECULES  . 42 

7.1  On  The  2II3  Shape  Resonance . 42 

7.2  Low  Energy  (0.01-20  eV)  e-C2H2  Cross  Sections . 46 


8.  EFFECT  OF  GAS  TEMPERATURE  ON  THE  CROSS  SECTIONS  .  69 

9.  ELECTRON-MOLECULE  SCATTERING  AT  HIGH  ENERGIES  .  .  74 


10.  CONCLUSIONS . 81 

11.  REFERENCES . 83 


APPENDIX:  PARTIAL  DIFFERENTIAL  EQUATION  APPROACH  TO 
ELECTRON,  POSITRON-MOLECULE  SCATTERING  BY  CHARLES  A. 
WEATHERFORD  . 89 


IV 


List  of  Figures 


Figure  Page  No. 

7.1.1  Eigenphase  sums  for  the  II?  e-C2H2  scattering  in  SCE  (curve  A),  AAFEGE 

(curve  B)  and  FEGE  (curve  C)  models . 52 

7.1.2  Eigenphase  sums  for  the  11^  e-C2H2  scattering  in  MSCEP  (curve  A),  MSCE 

(curve  B)  and  TFEGE  (curve  C)  models . 53 

7.1.3  Eigenphase  sums  for  the  IIfl  e-C2H2  scattering  in  MSCEP  (curve  A),  MSCE 

(curve  B)  and  TFEGE  (curve  C)  models.  Results  are  at  the  static-exchange 
(without  polarization)  level . 54 

7.2.1  DCS  for  the  e~— 2h2  scattering  at  1  eV.  Solid  line,  present  calculations; 

Crosses  and  open  boxes  (multiplied  by  a  factor  of  two),  experimental  data 
of  Kochem  et  al  (1985) .  55 

7.2.2  DCS  for  the  e~-2h2  scattering  at  2  eV.  Solid  line,  present  calculations; 

Open  boxes  (multiplied  by  a  factor  of  two),  experimental  data  of  Kochem 
et  al  (1985) .  56 

7.2.3  DCS  for  the  e~— 2h2  scattering  at  2.5  eV.  Solid  line,  present  calculations; 

Open  boxes  (multiplied  by  a  factor  of  two),  experimental  data  of  Kochem 
et  al  (1985) .  57 

7.2.4  DCS  for  the  e~-2h2  scattering  at  3  eV.  Solid  line,  present  calculations; 

Open  boxes  (multiplied  by  a  factor  of  two),  experimental  data  of  Kochem 
et  al  (1985) .  58 

7.2.5  Present  DCS  (solid  line)  as  a  function  of  energy  at  40°  angle  are  compared 

with  experimental  data  (multiplied  by  a  factor  of  two)  of  Kochem  et  al 
(1985) .  59 


7.2.6  Present  DCS  (solid  line)  as  a  function  of  energy  at  60°  angle  are  compared 

with  experimental  data  (multiplied  by  a  factor  of  two)  of  Kochem  et  al 
(1985) .  60 

7.2.7  Present  DCS  (solid  line)  as  a  function  of  energy  at  90°  angle  are  compared 

with  experimental  data  (multiplied  by  a  factor  of  two)  of  Kochem  et  al 
(1985) . 61 

7.2.8  Present  DCS  as  a  function  of  energy  at  angles  of  120°  and  150°.  ...  62 

7.2.9  DCS  for  the  e— 2h2  scattering  at  5  eV.  Solid  line,  present  calculations; 

Crosses,  experimental  data  of  Khakoo  et  al  (1993) .  63 

7.2.10  DCS  for  the  e-2h2  scattering  at  10  eV.  Solid  line,  present  calculations; 
Crosses,  experimental  data  of  Khakoo  et  al  (1993).  The  dash  curve  is 
the  calculations  of  Thirumalai  et  al  (1981),  while  open  circles  are  relative 


measurements  of  Hughes  and  McMillen  (1933) .  64 

7.2.11  DCS  for  the  e-2h2  scattering  at  15  eV.  Solid  line,  present  calculations; 
Crosses,  experimental  data  of  Khakoo  et  al  (1993) .  64 


7.2.12  DCS  for  the  e-2h2  scattering  at  20  eV.  Solid  line,  present  calculations; 

Crosses,  experimental  data  of  Khakoo  et  al  (1993) .  65 

7.2.13  Integral  cross  sections  for  the  e-2h2  collisions  in  the  0.01-20  eV  range. 

Present  calculations,  solid  line;  x,  experimental  data  of  Sueoka  and  Mori 
(1989);  o,  experimental  data  of  Briiche  (1929) .  66 

7.2.14  Present  results  on  the  momentum  transfer  cross  sections  for  the  e~-2h2 

scattering  in  the  0.01-20  eV  energy  range . 67 


VI 


PREFACE  AND  ACKNOWLEDGEMENTS 


This  work  was  accomplished  during  the  period  June  1990  through  1993  un¬ 
der  the  Air  Force  Contract  F33615-90-C-2032.  The  Air  Force  contract  man¬ 
ager  was  Lt  Patrick  D.  Kee  and  the  contract  officer  was  Rebecca  B  O’Kelley. 
This  technical  report  entitled  “ Electron  Interaction  with  Non-Linear  Polyatomic 
Molecules  and  their  Radicals  and  Ions"1'1  was  prepared  at  the  Physics  Department, 
Florida  A&  M  University,  Tallahassee,  Florida. 

We  would  like  to  thank  the  Florida  State  University  Supercomputer  Re¬ 
search  Institute  (SCRI)  for  providing  supercomputer  time  on  the  CRAY-YMP 
machine.  In  addition  supercomputer  time  on  the  NMFECC  CRAY-2  machine 
at  Livermore  was  also  made  available  to  us  through  SCRI,  which  is  partially 
funded  by  the  US  Department  of  Energy  through  contract  no.  DE-FC05- 
85ER250000.  We  extend  our  thanks  to  SCRI  staff  for  their  help  throughout 
the  execution  of  this  work.  We  would  like  to  take  this  opportunity  to  thank 
especially  Joe  Lannutti  (Director,  SCRI),  Chris  Lacher  (Chairman,  Resource 
Allocation  Committee),  Robert  Holden  (Manager,  Supercomputing  group), 
and  many  others  at  the  SCRI  and  innovation  park  facilities. 

A  part  of  the  present  research  work  was  performed  in  collaboration  with 
Professor  F.  A.  Gianturco  (Department  of  Chemistry,  University  of  Rome, 
Italy)  and  Dr.  D.G.  Thompson  (Department  of  Applied  Mathematics  and 
Theoretical  Physics,  The  Queen’s  University  of  Belfast,  N  Ireland).  A  travel 
grant  from  NATO  (CRG  890470)  is  thankfully  acknowledged.  A  research  col¬ 
laboration  with  Dr.  V.  Di  Martino  (Italy)  and  Nico  Sanna  Italy)  is  thank¬ 
fully  acknowledged  here.  It  is  also  a  pleasure  to  thank  our  colleagues  at  the 
Physics  Department,  FAMU.  We  are  grateful  to  Drs.  L.  Boesten,  T.  W.  Shyn, 
O.  Sueoka,  H.  Tanaka,  M.  Hayashi,  Steve  Buckman,  M  A  Khakoo  and  C. 
Szmytkowski  for  sending  their  experimental  data  prior  to  publication.  I  am 
also  thankful  to  Avesh  Jain  for  helping  in  preparing  some  tables  and  figures  of 
the  present  research  data  . 


vii 


1.  INTRODUCTION 


This  final  report  describes  the  work  done  on  the  project  entitled  “  Electron  Inter¬ 
actions  with  Non-Linear  Polyatomic  Molecules  and  their  Radicals  and  Ions  ” 
during  the  entire  period  of  June  1990  to  June  1993.  The  present  theoretical  research  was 
carried  out  under  the  auspices  of  Air  Force  contract  No.  F33615-90-C-2032.  The  work  was 
performed  at  the  Physics  Department,  Florida  A&M  University,  Tallahassee,  Florida  with 
the  direct  supervision  and  participation  of  Dr.  Ashok  K.  Jain,  Principal  Investigator  (PI), 
and  Dr.  Charles  A.  Weatherford  (Co-PI).  A  part  of  this  project  was  performed  in  collab¬ 
oration  with  Professor  F.  A.  Gianturco  (Department  of  Chemistry,  University  of  Rome, 
Italy)  and  Dr.  D.  G.  Thompson  (Department  of  Applied  Mathematics,  The  Queen’s  Uni¬ 
versity  of  Belfast,  N.  Ireland),  and  Dr.  Aaron  Temkin  (NASA-Goddard,  Greenbelt,  MD). 
The  collaboration  with  the  University  of  Rome  was  partially  funded  by  the  NATO  (under 
contract  No  CRG  890470)  in  terms  of  a  renewed  2-year  travel  grant.  The  work  with  Dr. 
Aaron  Temkin  was  performed  with  partial  support  of  NASA  grants  NAG-5307  and  NAGW- 
2930.  The  computational  work  was  fully  supported  by  the  FSU  Supercomputer  Research 
Institute  (SCRI)  in  terms  of  supercomputer  time  on  the  CRAY-YMP  (at  Tallahassee)  and 
on  the  NMFECC  CRAY-2  machine  (at  Livermore). 

The  overall  goal  of  this  project  was  to  develop  state-of-the-art  computer  programs 
and  to  perform  calculations  on  various  elastic  and  inelastic  cross  sections  for  low  energy 
electron  collisions  with  non-linear  polyatomic  molecules  and  their  ions  and  radicals.  It 
is  a  well  known  fact  that  such  a  project  involves  a  large  amount  of  computing  because 
of  complexities  involved  in  the  dynamics  of  the  collision  system.  Our  main  purpose  has 
been  to  produce  a  reliable  set  of  useful  cross-section  data  for  elastic  as  well  as  inelastic 
processes  without  involving  any  fitting  procedure.  When  electrons  interact  with  neutral 
molecular  targets,  time  considerations  of  the  electronic  and  nuclear  motions  allow  us  to 
treat  rotational,  vibrational,  and  electronic  processes  rather  independently.  This  is  the  so- 
called  adiabatic-nuclei-approximation  (ANA)  which  has  made  ab  initio  electron-molecule 
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calculations  a  dream  come  true.  Even  under  the  great  simplifications  introduced  by  the 
ANA  formulation,  the  treatment  of  electron-polyatomics  is  still  a  very  difficult  problem 
owing  to,  for  example,  the  multicenter  nature  of  the  interaction,  the  presence  of  non-local 
electron  exchange,  charge  correlation  and  polarization  effects  and  the  opening  of  several 
other  rearrangement  channels  . 

For  the  simple  case  of  elastic  (including  rotational  excitation  processes)  scattering  at 
the  static-exchange  level  (neglecting  short  range  correlation  and  long  range  polarization 
effects),  the  solution  of  inhomogenous  integro- differential  coupled  equations  iteratively  is 
an  arduous  task  for  the  case  of  non-linear  targets.  Consequently,  it  is  not  surprising 
that  just  a  few  years  ago,  most  of  the  electron-polyatomic  calculations  were  carried  out  by 
employing  model  potentials  for  exchange  and  polarization  effects  [1].  For  the  first  time,  our 
group  developed  scattering  codes  [2-7]  which  solve  electron-polyatomic  scattering  equation 
iteratively  to  treat  electron  exchange  effects  exactly.  To  implement  short  range  correlation 
and  long  range  polarization  effects  correctly  is  still  an  open  problem.  We,  however,  include 
these  charge  distortion  effects  non-empirically  using  the  perturbative  techniques. 

In  brief,  the  electron-molecule  problem  is  set  up  in  the  single-center-expansion  (SCE) 
scheme  under  the  body-fixed  (BF)  coordinate  system  in  the  ANA  close-coupling  formula¬ 
tion.  The  continuum  electron  function  is  obtained  numerically  for  each  irreducible  repre¬ 
sentation  of  the  molecule  point  group.  The  convergence  of  bound  and  continuum  functions 
with  respect  to  SCE  size  is  tested  properly  to  make  this  approach  meaningful.  We  have 
recently  developed  SCE  orbitals  from  multicenter  molecular  wavefunctions  [8].  Thus,  given 
a  molecule  of  any  symmetry,  we  should  be  able  to  investigate  low  energy  electron  scattering 
at  the  ab  initio  level. 

In  Section  2,  we  provide  a  summary  of  work  completed  so  far,  and  Section  3  gives  a  list 
of  published  (or  accepted)  and  submitted  papers  during  the  entire  period  of  the  contract. 
The  theory  of  the  present  calculations  is  presented  in  Section  4  and  Section  5  provides  the 
corresponding  numerical  details  and  a  list  of  computer  programs  employed  for  this  work. 
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In  Section  4,  we  also  discuss  our  recent  method  [6]  on  a  most  efficient  use  of  an  iterative 
scheme  in  terms  of  computer  time  without  sacrificing  any  numerical  accuracy. 

In  Section  6,  v*  summarise  our  main  calculations  on  the  rotationally  elastic,  inelastic 
and  summed  cross  sections  for  electron  collisions  from  several  molecules  such  as  the  CH4 , 
SiH4,  GeH4,  NH3,  H2O,  and  H2S  molecules  in  the  energy  range  of  0-20  eV.  A  detailed  set 
cross  sections  values  (Figures  and  Tables)  has  already  been  provided  in  previous  progress 
reports  [9-10].  For  all  these  species,  experimental  data  on  the  differential  (DCS),  integral 
and  momentum  transfer  cross  sections  are  available  for  comparison.  For  the  e-E^O  case, 
there  are  experimental  data  on  the  rotational  excitation  process  as  well.  In  Section  7, 
we  present  our  low  energy  electron  scattering  study  with  C2H2  molecules  by  employing 
a  similar  (as  mentioned  above)  technique  but  using  a  different  set  of  computer  codes  for 
linear  molecules.  The  present  e-C2H2  cross  sections  are  parameter-free  and  compare  very 
well  with  recent  measurements. 

In  Section  8,  we  discuss  our  main  findings  on  the  effect  of  gas  temperature  on  the  cross 
sections  in  connection  with  several  theorems  by  Shimamura  [11-17].  Section  9  summarises 
our  present  results  on  the  total  cross  sections  for  a  large  variety  of  diatomic  and  polyatomic 
molecules  from  about  20  eV  to  up  to  several  thousand  eV. 

Finally,  the  concluding  remarks  axe  made  in  Section  10.  The  last  section  is  presented 
as  an  Appendix.  In  this  Appendix,  Dr.  Charles  A  Weatherford  reports  his  work  on  the 
PDE  technique  in  electron-molecule  scattering. 
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2.  SUMMARY  OF  WORK  PERFORMED  (June  1990  -  June  1993) 

1.  The  electron  -  SiH4  scattering  at  0.01-20  eV  was  investigated  [5].  The  exchange 
(exactly)  and  polarization  (approximately  but  parameter-free)  effects  were  included 
to  make  a  direct  comparison  with  the  experimental  data.  This  is  a  detailed  study 
on  the  low-energy  behaviour  of  electrons  in  silane  gas.  Our  ab  initio,  parameter-free 
calculations  clearly  exihibit  the  cross  section  minimum  below  1  eV  and  a  shape- 
resonance  feature  around  2-4  eV;  this  is  in  agreement  with  measured  data.  No  any 
other  ab  initio  calculation  (with  polarization  effects)  on  the  e-SiH4  low-energy  cross 
sections  with  such  detr;ls  and  amount  of  cross  section  data  is  available  in  the  literature. 

2.  We  have  suggested  and  tested  a  new  optimized  iterative  scheme  to  solve  the  integro- 
differential  equations  more  economically  In  order  to  treat  exchange  effects  exactly  in 
an  efficient  way.  The  new  scheme  is  more  suitable  for  the  case  of  polar  molecules 
where  convergence  problems  are  severe.  In  the  optimized  scheme,  the  continuum 
wavefunction  of  the  ele  rion  with  full  iterative  procedure  is  calculated  only  once  at 
a  selected  energy,  say  E ,  i.e.,  fw(r\  E).  The  fu'(r\  E)  is  employed  as  a  starting 
point  for  iterative  scheme  at  other  energies.  We  found  for  the  H2O  case  (and  very 
rcently  for  the  NH3  case  also)  that  the  number  of  iterations  required  in  the  optimised 
scheme  is  smaller,  by  a  factor  of  two  to  three,  than  the  usual  one  [6].  Further,  we 
noticed  that  the  exact-exchange  treatment  is  necessary  for  low  partial  waves  only.  The 
non-penetrating  higher  partial  waves  can  be  described  by  a  simple  model  exchange 
potential  under  the  density  functional  theory.  Thus,  the  new  scheme,  which  is  very 
economical  and  essential  to  treat  polar  molecules,  can  be  described  as  follows.  For  low 
partial  waves  (say  t  —  £%),  we  consider  full  exact-exchange  calculation;  for  t>eE, 
say  up  to  £  =  £m,  we  employ  a  model  exchange  (for  example  the  free-elect ron-gas 
exchange  )  approximation;  and  finally  for  very  high  partial  waves  £  >  £m,  say  l  —  1\jb, 
we  use  unitarised  Born  approximation.  A  paper  on  the  optimized  iterative  scheme 
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with  test  calculations  on  the  e-H20  system  has  been  published  recently  [6]. 

3.  The  effects  of  gas  temperature  were  investigated  on  the  rotational  excitations  in  a 
molecule  (CH4  and  SiH4)  by  slow  electron  impact.  We  have  found  that  the  rotational 
excitation  process  is  sensitive  to  gas  temperature  in  an  experimental  situation  and 
one  should  be  careful  in  comparing  theoretical  rotational  excitation  cross  sections 
with  measurements.  In  particular,  for  example,  the  theory  predicts  zero  cross  sections 
at  zero  angle  for  the  0  — ►  3  rotational  transition  in  a  spherical  top  molecule,  while 
measurements  (on  the  CH4  molecule)  give  non-zero  DCS  in  the  forward  direction  [18]. 

4.  We  studied  [7]  various  approximations  to  include  polarization  effects  in  the  e-SiH4  col¬ 
lisions  at  low  energies.  In  this  calculation,  we  have  employed  two  types  of  parameter- 
free  model  polarization  potentials  alongwith  exact-exchange  effects.  Our  conclusion 
is  that  a  polarization  potential  determined  from  polarised-orbital  type  approach  is  a 
better  model  than  the  one  from  the  target  density  functional  theory. 

5.  The  low-energy  electron  scattering  with  H2O,  H2S  and  NH3  gases  was  investigated 
under  the  same  theory  as  mentioned  above.  We  tested  convergence  of  the  differential 
and  total  cross  sections  for  these  polar  gases,  for  which  the  ANA  theory  fails  in  the 
forward  direction.  The  final  detailed  results  on  the  rotationally  elastic,  inelastic  and 
summed  cross  sections  (differential,  integral,  momentum  transfer  and  energy  loss)  for 
both  the  polar  gases  (H2O  and  H2S)  were  completed.  This  work  is  being  written  up 
in  separate  papers. 

6.  We  also  developed  a  single-center  computer  program  [8]  to  expand  molecular  quan¬ 
tities  (orbitals,  density,  static  potential  etc.)  around  the  center-of-mass  (COM)  of 
the  target  from  the  calculations  of  Quantum  Chemistry  codes  such  as  the  HONDO  or 
GAUSSIAN.  This  work  was  done  in  collaboration  with  the  group  at  Chemistry  De¬ 
partment,  University  of  Rome,  Italy.  We  made  attempts  to  study  some  big  molecules 
such  as  the  GeH4,  CF4,  SF6,  etc..  The  first  calculations  on  the  e-GeH4  collisions  were 
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carried  out  and  results  were  published  [19].  We  have  also  finished  some  calculations 
on  the  electron-SF6  scattering  at  low  energies. 

7.  In  order  to  see  the  behavior  of  electron-molecule  cross  sections  at  higher  energies 
(from  about  20  eV  to  up  to  several  thousand  eV),  we  have  adopted  a  simple  but 
accurate  scheme  to  determine  total  (elastic  plus  all  inelastic  channels)  cross  sections 
for  a  large  class  of  molecules.  A  comparison  is  made  with  experimental  data  for  all 
the  gases  [20]. 

8.  We  carried  out  extensive  calculations  on  the  low  energy  behaviour  of  electrons  in  C2H2 
molecule.  The  cross  sections  (differential  as  well  as  integral)  compare  very  well  with 
experimental  data,  particularly,  the  shape  resonance  phenomenon.  Two  papers  were 
published  on  this  work  [21-22]. 

9.  Finally,  our  state-of-the-art  computer  codes,  some  modified  during  this  period  of 
contract,  are  in  a  situation  to  carry  out  low  energy  electron  scattering  calculations  for 
any  molecule  of  reasonable  size  and  symmetry.  There  are  a  large  number  of  molecules 
for  which  no  theoretical  study  exists,  while  the  cross  sections  are  needed  in  many 
applied  sciences. 

The  computing  facilities,  provided  by  the  FSU  Supercomputer  Research  Institute 
(SCRI),  have  been  excellent  in  order  to  carry  out  the  above  programs  effectively.  We  have 
made  use  of  supercomputers  (CRAY-2  at  the  NMFECC,  Livermore  and  CRAY-YMP  at 
the  SCRI,  Tallahassee)  through  the  FSUCC  allocation  committee.  Most  of  the  graphical 
work  and  small  calculations  were  carried  out  at  the  FSU  and  SCRI  VAX  machines  which 
serve  as  front-end  system  for  the  supercomputers. 
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4.  REVIEW  OF  THEORETICAL  APPROACH 


When  low-energy  electrons  interact  with  molecular  targets,  severed  processes,  elastic 
and  inelastic,  can  take  place  in  the  present  energy  region  (  E  <  20  eV  ): 

e  AjiSjji  *  e  +  AjiBfjx  ( clastic ) 

— ►  e~  -f  ( AnBmy  ( inelastic ;  vibrational  and/or  rotational ) 

(There  are  also  other  inelastic  processes  such  as  the  electronic  excitation,  dissociative, 
associative,  or  ionization  channels.)  The  basic  Schrodinger  equation  involves  an  interaction 
Hamiltonian  and  the  corresponding  wave  function  of  the  total  system: 

=  (4.1) 

where  the  electron-molecule  Hamiltonian  is  given  by: 

Hint  —  Htitcij^i't  R)  d-  VJnt(r,  r,-,  R)  +  Hnuc( R).  (4-2) 

Here  r  represents  the  projectile  electron  coordinate  and  and  R  collectively  represent  the 
target  electronic  and  nuclear  coordinates,  respectively.  It  is  not  possible  to  solve  equation 
(1)  exactly  for  any  of  the  channels  or  collision  systems.  We  first  describe  the  iterative 
techniques  to  solve  integro-differential  (see  below)  equation  under  the  SCE  formalism. 

4.1.  Single-Center-Expansion  (SCE)  Iterative  Method 

In  order  to  obtain  various  elastic  and  inelastic  cross  sections  (differential  as  well  as 
integral),  we  make  use  of  the  adiabatic-nuclei-approximation  (ANA)  [23]  and  set  up  our 
scattering  equations  in  the  SCE  scheme  under  the  close-coupling  formalism.  The  FNA  is 
generally  a  very  good  approximation  except  in  certain  cases:  for  example,  near  threshold 
and  sharp  shape-resonance  energy  regions,  for  polar  molecules,  etc. 
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In  the  body-fixed  (BF)  frame  of  reference,  the  time-independent  Schrodinger  equa¬ 
tion  of  the  electron-molecule  system  reads  as  (omitting  the  parametric  dependence  on 


intemucleax  coordinates  R): 

[-^V2  +  tfelecCr)  +  Mnt  -  £]*t( r,n)  =  0,  (4.1.1) 

where  the  operators  in  the  bracket  axe  respectively  the  kinetic  energy  operator  of  the 
scattered  electron,  the  unperturbed  molecular  target  Hamiltonian,  the  interaction  potential 
of  the  electron-molecule  complex  and  the  total  energy  of  the  system.  The  total  wave 
function  \&t  can  be  expanded  such  that: 


=  A'22$i(ri,r2,...rz)Fi(r),  (4.1.2) 

«=i 

where  are  the  target  electronic  states  and  any  other  suitable  pseudostates  that  can 
represent  the  target  response  function  to  the  polarization  of  the  molecule  because  of  the 
incoming  electron.  The  continuum  function  l'i(r)  describes  the  projectile  motion  in  molec¬ 
ular  state  i  and  A  is  the  usual  antisymmetrization  operator  to  ensure  Pauli  principle.  In 
the  expansion  (4.1.2),  we  have  neglected  the  correlation  functions  of  the  ( N  +  l)-electron- 
molecule  complex. 

We  can  assume  that  the  center  of  the  electron-molecule  system  is  at  the  center-of- 
mass  (COM)  of  the  target.  The  spin  of  the  scattered  electron  is  coupled  with  the  total  spin 
of  the  molecule  to  form  the  eigenstate  of  S2  and  its  projection  along  the  symmetry  axis  ( 
say  Sz  )  corresponding  to  the  total  spin  quantum  number  5  and  its  substate  Ms  which  are 
constants  of  motion  during  the  collision.  Each  of  these  eigenstates  can  be  expanded  over 
the  set  of  angular  functions  according  to  the  prescriptions  of  a  BF  frame  of  reference  which 
provides  the  basis  of  an  irreducible  representation  (IR)  pp  of  the  molecular  symmetry  point 
group.  Thus  we  can  write  the  total  wave  function  as: 
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(4.1.3) 


K“SUs  =  E  E  r2  . . .  rz,  f,  (r), 

1  =  1  th 

where  a  is  the  spin  variable  of  the  scattered  electron.  The  functions  $  in  equation  (4.1.3) 
are  defined  as: 


pifiiSiMsj 


(r1,r1„..rzW(f'),(i .m.jcj  4 


2 

S;  m, 


S 


(4.1.4) 


where  r/’s  are  electron  spin  functions  for  the  scattered  electron  and  C’s  are  the  Clebsch- 
Gordon  coefficients.  The  bound  orbitals  <f>i  describe  Z-electron  target  wave  function  used 
in  expansion  (4.1.4),  each  transforming  like  a  particular  IR  of  the  molecular  point  group 
and  with  a  specific  total  spin  eigenstate.  The  main  molecular  symmetry  axis  defines  the 
Ms  direction. 

As  a  standard  technique,  we  now  project  the  scattering  equation  onto  the  channel 
functions  and  obtain  the  well  known  coupled  integro-differential  equations  for  each  scat¬ 
tered  electron  channel  function  ( pfiS ): 


U{li  +  1) 

9*2 


(4.1.5) 


4.2  Static  and  Exchange  Interactions 

The  direct  potential  matrix  Vp,tS  defines  the  coupling  between  two  channel  functions 
through  the  Coulomb  operator  (for  details  see  Ref.  1).  The  exchange  interaction  is  repre¬ 
sented  by  the  non-local  term,  WpmS,  which  gives  rise  to  the  integro-differential  nature  of 
the  coupled  scattering  equations.  We  do  not  give  all  the  details  here  for  obvious  reasons 
(see  Ref.  1  for  full  details).  For  a  closed-shell  molecule,  belonging  to  the  totally  sym¬ 
metric  1  A\  IR  of  the  molecular  point  group,  equation  (4.1.5)  does  not  depend  on  the  spin 
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quantum  numbers.  The  equation  (4.1.5)  is  known  as  the  static-exchange  interaction  of  the 
electron-molecule  system  if  we  retain  only  the  ground  state  (i  =  1)  of  the  target.  From 
a  numerical  point  of  view,  it  is  convenient  to  include  the  target  distortion  (polarization) 
effects  via  a  model  potential  added  to  the  direct  potential  matrix  in  equation  (4.1.5).  At 
short  distances,  the  correlation  effects  cam  be  treated  in  a  similar  fashion.  The  exchange 
term  in  equation  (4.1.5)  is  basically  of  the  following  type: 

W"”‘=E/ «<*)!«•  (4.2.i) 

Recently  a  program  has  been  developed  to  treat  WPft  exactly  under  the  iterative  pro¬ 
cedure  [2].  Under  the  static-exchange  plus  (model)  polarization  (SEP)  model,  the  integro- 
differential  coupled  equation  to  be  solved  for  a  closed-shell  system  for  each  symmetry  pp 
is  rewritten  from  Eq.  (4.1.5)  as  (suppressing  the  pfi  symbol  for  simplicity): 

[|1  -  £!£+i!  +  fc2]F/A(r)  =  2£[Vg>„,(r)  +  Vt\%,  +  wn,„,(r)1,  (4.2.2) 

ar  T  t'h 1 

The  potential  matrix  t,h,  in  Eq.  (4.2.2)  is  determined  from  the  following  expression 
for  the  static  potential: 

/Z  M 

|*o|2{£  |r  -  rj  r<fri<ir2  •  •  •  drz}  -  £  Z> |r  -  Rif1 .  (4.2.3) 

j=i  1=1 

where  $o  Is  the  target  ground  state  wavefimction  given  as  a  single  Slater  determinant  of 
one-electron  Z  spin  orbitals  <£,-( r)  and  M  the  number  of  nuclei  in  the  molecule. 

In  a  model  exchange  problem,  the  non-local  exchange  term  in  the  right-hand  side  of 
Eq.  (4.2.2)  can  be  replaced  by  a  local  potential,  say  V(£tlh,(r).  A  very  popular  form  of 
Vex(r)  has  been  the  free-elect ron-gas-exchange(FEGE)  approximation  which  was  mod¬ 
ified  by  Hara  [24]  (to  be  denoted  by  HFEGE)  for  electron-molecule  collisions.  In  the 
HFEGE,  the  Vex  reads  as: 


12 


V“(r)  =  |Mr){i  +  iXrl“lT±!l}.  (4-2.4) 

7T  2  4  7?  1  -  V 

where,  kp( r)  =  (37r2p0(r))»  and  T)  =  (k2  +21  +  k2F)%  jkp.  Here  I  is  the  ionization  potential 
of  the  target  and  po(T)  is  the  unperturbed  charge  density  of  the  molecule.  Several  versions 
of  HFEGE  are  available  by  varying  the  value  of  I  [25].  Thus,  writing  the  integro-differential 
equation  (4.2.2)  in  a  convenient  matrix  form,  LF  =  WF  where  WF  is  the  exchange  term, 
the  iterative  scheme  is  LF'  =  WF*-1,  where  i  =  0,1, . . .  .  In  order  to  start  the  solution, 
we  chose  F°  to  be  the  solution  obtained  from  the  HFEGE  potential  (4.2.4).  Our  earlier 
scattering  code  (POLY,  see  Section  5.3)  has  been  modified  for  the  solution  of  Eq.  (4.2.2). 
The  new  code  (POLYEX)  is  now  working  successfully  on  the  CRAY-2  and  FSU’s  VAX 
and  CRAY-YMP  systems. 


4.3  Polarization  Interaction 

It  is  very  hard  to  include  polarization  effects  accurately;  however,  the  asymptotic  form 
is  known  exactly,  i.e.: 


VpdM,  i)  =  +  <«(!*)*  SS1  +  <4(^)4  Sf1],  (4.3.1) 

where  the  S™q  is  a  real  spherical  harmonic  (see  Ref.  1  for  its  definition  and  various 
properties  ),  (r,0,<j>)  are  the  coordinates  of  the  projectile  referring  to  the  center  of  the 
target  and  the  spherical  (ao)  and  nonspherical  (q2  and  a2  )  polarizabilities  are  expressed 
in  terms  of  the  polarizability  tensor  of  the  target,  namely: 


1/  ,  2.  1  1  .  , 

Q° =  3^ai1  +  0122  +  a 33':  02  =  jr  33  —  2au  ~  ~2°L22  ’  0:2  =  ai1  —  a22’ 

In  order  to  represent  the  Vpo1  at  short  distances,  one  has  to  introduce  some  kind  of 
adjustable  parameter(s)  in  the  cut-off  function  which  is  multiplied  to  the  right-hand  side 
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of  equation  (4.3.1).  The  value  of  the  unknown  parameter  is  determined  by  the  fitting  pro¬ 
cedure.  An  important  development  was  made  by  our  group  in  evaluating  an  approximate 
parameter-free  polarization  potential  [26]  based  on  the  second-order  perturbation  theory 
[27]  and  the  criterion  of  Temkin  [28].  In  this  method,  the  molecule  is  supposed  to  be  in  an 
electric  field  E  =  Er  (  produced  by  the  incoming  electron)  at  r(r0<£),  which  results  in  an 
extra  potential  energy  V,  the  leading  term  of  which  is  given  by: 

z 

V(ri,r2,...rz;E)  =  E 

1=1 

z 

^E^riTTi  (4.3.2) 

i=i 

The  energy  of  the  molecule  in  the  electric  field  can  be  expanded  as  a  Taylor  series  in 
Ei(i  =  1,2,3),  the  components  of  E  along  the  axes  of  our  coordinate  system.  The  second 
order  energy  reads  as: 


W2  =  l'£Ei<-iiEi,  (4.3.3) 

•J 

In  the  method  of  Pople  and  Schofield  [27],  the  first  order  wave  function  is  expanded 
in  terms  of  ground-state  $o  functions  as: 

Z 

*l -*«£/(*)•»,  (4.3.4) 

«=1 

where  unknown  coefficients  /’s  are  determined  variationally  from  the  following  expression 
of  W2,  i.e.: 


W2  =  (*,| Hmol  -  W0|$i)  +  2<$oM$i>-  (4-3-5) 

The  success  of  this  parameter-free  model  for  the  polarization  potential  (to  be  denoted 
as  the  JT  polarization  approximation)  along  with  an  exact  treatment  of  exchange  has  been 
exceptional  in  our  recent  investigation  of  the  low-energy  electron-CH.*  interactions  [3]. 
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4.4  Cross-Section  Formulae 


We  now  summarize  our  cross  section  formulae.  As  a  special  case,  we  will  provide 
explicit  expressions  in  the  simple  case  of  spherical  molecules  such  as  the  CH4,  SiH4  etc. 
For  full  details  we  refer  to  our  review  article  [1].  In  this  project,  we  have  investigated 
mainly  the  rotationally  elastic,  inelastic  and  summed  channels  for  differential,  integral, 
momentum  transfer  and  energy-loss  cross  sections.  As  mentioned  above,  we  solve  our 
scattering  problem  in  the  fixed-nuclei  approximation  under  the  BF  frame  of  reference.  In 
order  to  obtain  physical  parameters,  we  transform  the  BF  scattering  amplitude,  /( k  •  r), 
(  k  and  r  are  respectively  the  initial  and  final  directions  of  the  projectile  )  into  the  space- 
fixed  amplitude  /( k  •  r';  ctfii)  (  where  r'  now  refers  the  direction  of  scattered  electron  with 
respect  to  SF  coordinate  system  and  (a/fy)  are  the  three  Euler  angles  )  by  making  use  of 
rotation  matrices.  The  SF  scattering  amplitude  is  given  by: 

/( k.r':o/37)=  £  i*(2<  +  l 

Ihmt'h'm'  Xpti 

x  dL(«WDL(«/57)(S?;,„'  -  ««.&».),  (4.4.1) 

where  b's  are  defined  earlier1,  D’s  are  the  rotation  D-matrices  and  S  is  the  scattering 
S-matrix.  In  the  impulse  approximation  [29],  the  scattering  amplitude  for  a  particular 
rotational  transition  JKM  — ►  J'K'M1  is  written  as: 

f(JKM  -  J’K’M')  =  tyjKM  !/(&•?';  a0j)  O’ J'K'M'),  (4.4.2) 

where  \ip jkm)  Is  the  rotational  eigenfunctions  of  a  molecule  (asymmetric,  symmetric  or 
spherical  tops)  which  are  given  in  standard  texts.  Expressions  for  a  general  asymmetric  top 
are  discussed  in  Ref.  1.  Here  we  discuss  formulae  for  spherical  tops  only.  The  differential 
cross  section  for  the  J  —*  J'  transition  is  obtained  by  summing  over  all  final  magnetic 
substates  K'M'  and  averaging  over  all  initial  substates  KM,  i.e.: 
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(4.4.3) 


±<j_  J')  =  - _ ! _ 

d(i{  >  *(2j  +  i)2KMi?:M, 


\f(JKM  J'Jf'M')!2, 


where  is  the  wavenumber  of  the  scattered  electron  given  by  the  relation: 


2k'2  =  2 k2  +  Ej  -  Ej.  .  (4.4.4) 

For  a  spherical  top  molecule,  the  energy  of  the  Jth  level  is  given  by  Ej  —  BJ(J  +  1), 
where  B  is  the  rotational  constant  of  the  molecule  in  question.  It  is  convenient  to  express 
the  DCS  in  terms  of  Legendre  polynomials: 

|j(J -■/•)  =  £  53+  Al(J  J')PUcost).  (4.4.5) 

L 

Here  6  is  the  scattering  angle  between  the  vectors  k  and  r'.  The  expansion  Ai 
coefficients  are  found  to  be  [30]: 


Ai{j  -  J')  =  (2J'+4^((22j^.1i))(  £K2i  +  W21'  +  W2'  +  1X2p  +  l)]*i'_,'(-0I_r 


ll'll 

j+r 


(o  0  o)(o  0  o)  E  (-iyw( irrt-,iL)MpMp'  (4.4.6) 


where  the  M-matrix  is  defined  as: 


M’P=  53  (-1  rvdl  l  TS>V  ■  (4.4.7) 

mm'  hh'pfi  '  ' 

Here,  as  visual,  the  T-matrix  is  defined  in  terms  of  the  S-matrix  as  T  =  (S  —  1)  where 
S  =  (1  +  iK)(l  —  iK)-1.  The  scattering  K  matrix  for  each  symmetry  (pp).  In  equation 
is  a  3  —j  symbol,  W{abcd\  ef )  is  the  well  known  Racah  coefficient  and 
the  tfhm  coefficients  are  the  expansion  terms  in  the  definition  of  symmetry-adapted  basis 
functions  in  terms  of  real  spherical  harmonics  [1].  From  equation  (4.4.5),  it  is  easy  to  find 
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simple  forms  of  the  total  (a()  and  momentum  transfer  (<rm)  cross  sections  in  terms  of  Al 
coefficients,  namely: 


aJtJ'  =  4ttA0(7  -4  J')  ;  oJmJ' =An{AQ(J ->  J')).  (4.4.8) 

It  is  easy  to  determine  A0  and  Ax  coefiicients  for  at  and  om  rather  than  the  full 
expansion  (4.4.6)  for  the  DCS.  For  example,  the  Ao  and  A\  are  given  by  the  following 
simpler  expressions: 

Ao  =  Y1  (-1)'+,'| M’p  \\  (4.4.9) 

tt'jmj 

and: 

-4.  =  £  (-1)'+''+)l«<  +  +  l)}*WVrt  +  W  +  1; 

"'jmj 

+  {£{£'  +  1  )}±W{££'£  -  1£'  + 

+  {(£  +  1  )£'}X>W(££'£  +  l£'  - 

+  (££')±W(££'£  -  W  -  MJt%_x]  (4.4.10) 

where  g(JJ')  represents  the  statistical  (2 7'  +  l)/(2 7  +  1)  factor.  The  vibrationally  elastic 
(rotationally  summed)  cross  sections  me  obtained  by  summing  over  7'  for  7  =  0.  The 
selection  rule  for  the  transition  7  — *  7'  can  easily  be  worked  out  from  the  asymptotic  form 
of  the  static  potential  (4.2.3)  which  transforms  as  1  Ax  symmetry  of  the  molecular  point 
group.  For  any  specified  values  of  7  and  7',  the  angular  momentum  transfer  j  (Eq.  4.4.6) 
takes  the  following  values: 

|7-  j'\<j<  7  +  7'  ;  <l  +  l',  (4.4.11) 

In  addition,  calculations  on  the  polar  polyatomic  molecules  (  H2O  and  NH3  )  need 
further  attention  in  order  to  remove  the  deficiencies  of  the  ANA  formalism  towards  conver¬ 
gence  of  the  differential  and  total  integral  cross  sections  [see  Refs.  31-32].  We  employ  the 
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multipole-extracted- adiabatic- nuclei  (MEAN)  approximation  [32].  In  the  MEAN  method¬ 
ology,  a  maximum  use  has  been  made  of  the  first-Bom- approximation  (FBA)  for  the  higher 
symmetries  and  partial  waves.  So  far,  such  a  program  exists  only  for  the  linear  molecules 
and  we  have  modified  our  previous  code  [see  Section  5.3].  Finally,  in  the  MEAN  approxi¬ 
mation,  the  differential  cross  sections  (DCS)  for  a  rotational  transition  J  — ►  J'  are  written 
as  [32]: 


dvJJ'  _  A, ri»A 
dQ.  dQ 


i  -^m»x 

+  DC(J,'7';00>f'l5  £(*V.  -  BlT)Px(c°s«), 

J  i  \ _ ft 

(4.4.12) 


A=0 


where  the  first  term  is  the  usual  closed  form  for  the  rotational  excitation  ( J  — ►  J')  DCS 
in  the  space-fixed  FBA;  kj  and  kj<  are  respectively  the  wave  vectors  in  the  initial  and 

final  channels;  C(  .  )  is  a  Clebsch-Gordan  coefficient;  lt  is  the  angular  momentum 

transfer  during  the  collision;  B\ti,  are  the  DCS  expansion  coefficients  ,  expressed  in  terms 
of  scattering  matrix,  and  are  the  corresponding  quantities  in  the  FBA  evaluated  in 

the  BF  coordinate  system 

The  expression  for  the  space-fixed  FBA  quantities  are  given  in  Ref.  1.  The  important 
quantity  is  the  FBA  differential  and  integral  cross  sections  to  be  employed  in  equation 
(4.4.12).  Here  we  give  these  expressions  explicitly: 

k'  (2 J'  +  1) 

TTk  (2 J  +  1)  ^ 

and  the  integrated  quantity  is  expressed  as: 


—(j 

dir 


rrx> 

J'  =  -t  ~r ~  ^2(2i+  l)"1  J  drr2vtm(r)jt(qr ) 


(4.4.13) 


2(2  J'  + 1) 


rk+k'  ft  oo 

V'  (2^  f  l)-1  /  dqq  /  drr‘2vem(r)jt(qr ) 


(4.4.14) 


where  dipole  and  quadrupole  terms  correspond  respectively  from  J  =  0  to  J'  =  1  and 
2.  It  is  to  be  noted  here  that  except  at  very  low  energy  (below  0.1  eV)  the  short-range 
form  of  V(m(r)  coefficients  is  very  important.  We  have  included  these  coefficients  exactly 
at  all  radial  distances,  rather  than  employing  the  asymptotic  form  in  terms  of  dipole, 
quadrupole,  etc.,  moments. 
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5.  REVIEW  OF  NUMERICAL  DETAILS 


Here,  we  provide  a  brief  review  of  numerical  and  computational  aspects  of  the  present 
approach  employed  for  electron  scattering  from  molecules.  The  basic  ingredients  are  the 
molecular  wavefunction  and  the  Coulomb  interaction  to  describe  bound  and  bound-free 
quantities.  A  near-Hartree-Fock  description  of  the  target  is  sufficient  for  the  present 
purpose.  For  full  details  we  refer  our  earlier  work  [1,  9-10].  For  completeness  purpose, 
here  we  summarize  the  target  bound  functions,  iterative  procedure,  and  finally  a  fist  of 
computer  programs  employed  for  the  present  calculations. 

5.1  Molecular  Wavefunctions 

For  any  polyatomic  molecule,  the  multi-center  target  wavefunctions  are  easy  to  deter¬ 
mine  using  standard  Quantum  Chemistry  codes.  In  a  single-center-expansion  approach, 
however,  determining  actual  electronic  densities  around  the  bound  nuclei  still  rquires  a 
large  computational  effort  in  order  to  represent  correctly  or  with  high  precision  the  total 
charge  distribution.  Relativistic  effects  are  neglected  in  our  present  work  and  the  sepa¬ 
ration  of  electronic  and  nuclear  motions  is  assumed  under  the  usual  Bom-Oppenheimer 
approximation  (BOA).  To  understand  the  nature  of  polyatomic  electrons,  it  is  necessary  to 
consider  the  symmetry  properties  in  terms  of  symmetry-adapted  wave  functions  (SAWF). 

For  our  purpose  to  study  molecular  species  like  CH4,  SiH4,  GeH4,  SF6,  CF4,  H2O, 
H2S,  NH3,  etc.,  the  SCE  approach  is  very  successful.  We  employed  a  computer  code 
MOLMON  (see  Section  5.4)  to  generate  bound  orbitals  at  the  fixed  equilibrium  geometry 
of  the  target.  In  general,  the  molecular  orbitals  axe  written  in  terms  of  LCAO  (linear 
combination  of  atomic  orbitals)  centered  at  the  COM  of  the  molecule.  In  order  to  represent 
molecular  states,  we  follow  the  nomenclature  of  Herzberg  [33].  In  brief,  the  capital  letter 
A  is  used  to  represent  one-dimensional  IR  with  character  +1,  while  for  character  —  1, 
the  letter  B  is  used.  For  a  two-dimensional  IR,  the  letter  E  is  used,  while  a  three- 
dimensional  IR  is  denoted  by  letter  T.  Suffixes  1,2,3..  are  used  to  distinguish  different 
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IR’s  with  the  same  dimensions.  The  subscripts  g  or  u  differentiate  between  even  (gerade) 
or  odd  (ungerade)  representations  of  the  same  class  if  inversion  symmetry  also  exists  (for 
example  in  the  case  of  SFg  molecule).  Lower  case  letters  (a,  b ,  e,  t  etc.)  axe  used  to  denote 
individual  orbitals.  For  example  the  ground  state  of  the  H2O  molecule  is  written  as  follows: 
la\  2a\  Za\  lb\  lb* 

For  the  present  set  of  molecules,  we  have  included  up  to  £max  =  7  in  the  SCE  expansion 
of  bound  and  continuum  orbitals.  The  molecular  quantities  (such  as  the  total  energy, 
various  multipole  moments,  etc.)  obtained  with  the  present  SCE  basis  compare  reasonably 
well  with  the  experimental  and  other  sophisticated  calculations  (see  Ref.  9). 

5.2  SCE  Parameters  from  Multicenter  Target  Wavefunctions 

Recently  our  group  [8]  has  developed  an  interface  computer  code  to  generate  SCE 
orbitals  from  multicenter  target  wavefunctions.  First,  the  bound  and  continuum  orbitals 
around  the  COM  of  the  target  axe  written  as  follows: 

Mn)  =  £  r->uS,(r)X'r-(t),  (5.2.1) 

hi 

and: 


**'(*)  =  £  r-'K^Xi (5.2.2) 

hi 

Here  a  labels  a  specific  multicenter  MO.  In  order  to  perform  expansion  (5.2. 1-2),  we 
need  to  start  from  the  multicenter  wavefunction  which  describes  the  molecule  and  then 
generate  the  Uki(r)  coefficients  by  quadrature,  from  these  coefficients  we  can  determine 
the  corresponding  total  density  in  its  expansion  form: 

p(r)  =  r~lPhi(r)Xhi(r),  (5-2.2) 

hla 


20 


In  the  present  work  we  started  with  four  molecular  systems,  belonging  to  Tj  symmetry. 
These  are  the  SiH4,  GeH4,  SF6  and  CF4  molecules.  The  SFg  molecule  has  extra  inversion 
symmetry  also.  The  presence  of  the  large  electronic  density  away  from  the  COM  in  CF4, 
makes  the  study  of  numerical  convergence  of  the  expansion  particularly  interesting.  The 
initial  HF-SCF  calculations  were  carried  out  using  the  HONDO  package  [34]  and  the 
general  features  of  the  relevant  optimized  geometries  and  basis  set  expansions  were  listed 
in  Ref.  8.  We  see  that  the  triple-zeta  plus  polarization  (TZ+P)  functions  were  employed 
for  SiH4  and  CF4,  thus  yielding  a  fairly  accurate  description  of  the  ground  state  electronic 
charge  density,  while  the  GeH4  target  was  treated  at  the  split-valence  level  of  accuracy 
and  involved  66  contracted  functions.  Each  total  charge  density  was  then  employed  in  a 
quadrature  scheme  that  produced  the  symmetry  adapted  coefficients  of  Eq.  (5.2.2)  and 
each  of  the  MO  contributions  of  Eq.  (5.2.1).  The  numerical  procedure  starts  from  a  set 
of  points  in  the  space  r  =  (1,  y,  z)  which  is  the  one  generated  by  the  HONDO  algorithm, 
and  uses  them  to  perform  a  transformation  to  polar  coordinates  (r,  0,  <f>).  It  is  the  number 
of  latter  points  which  is  employed  to  construct  the  integration  grids  for  Eqs.  (5.2.1)  and 
(5.2.2)  that  yields  the  required  coefficients.  The  Gauss-Laguerre  integration,  for  each 
coefficient  in  both  equations,  was  carried  out  over  a  range  of  500  r  values,  from  0.001  a<>, 
up  to  about  15  a0,  with  a  variable  grid  that  became  denser  around  the  origin  and  around 
the  outer  nuclear  positions.  In  the  radial  region  with  r  <  1.0a0,  and  for  i  <  12,  the 
^-integration  and  ^-integration  were  carried  out  on  a  grid  of  20  x  20  =  400  pivots.  As  the 
r-value  increased,  the  grid  size  went  up  to  its  maximum  value  of  9216  pivots  (96  x  96). 
The  test  on  the  original  values  of  either  the  orbitals  or  of  the  densities  showed  agreement 
of  ~  0.1%  for  the  expansion  with  l  <  20  in  the  multipolar  coefficients,  at  typical  distances 
either  around  the  mid-bond  or  at  the  outer  molecular  nuclei. 

A  comparison  of  the  behavior  of  the  radial  density  coefficients  for  the  three  systems 
examined  were  shown  in  Fig.  5.2.1  of  Ref.  10,  where  the  radial  region  around  the  outer- 
atom  position  was  shown  in  each  case.  Our  conclusions  are  the  following  : 
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1.  The  lightest  molecule,  silane,  exhibits  quite  fast  convergence,  even  over  the  rather 
difficult  region  of  the  hydrogen  cusps  around  the  outer  nucleus  position  at  ~  2.787 a0. 
This  is  due  to  the  fairly  compact,  nearly  spherical  structure  of  a  system  in  which  a 
rather  low  electron  density  is  involved  in  the  bonding  and  around  the  H  atoms  away 
from  central  Silicon  atom. 

2.  The  feature  of  rapid  convergence  is  also  exhibited  by  the  germane  molecule  (bottom 
part  of  the  figure  5.2.1)  where,  being  at  even  larger  distance  ( R  =  2.905ao)  than  silane, 
the  location  of  the  hydrogen  atoms  affects  little  the  essentially  spherical  structure  of 
the  central  atom,  where  more  than  90%  of  the  charge  density  is  located. 

3.  The  situation,  on  the  other  hand,  is  quite  different  in  the  case  of  CF4  molecule: 
large  electron  density  portions  exist  around  the  outer  Florine  atoms  and,  therefore, 
the  convergence  in  t  is  much  slower.  We  see  that  at  distances  corresponding  to  the  F 
atom  positions  (2.45ao),  that  all  coefficients  are  much  stronger  than  their  counterparts 
for  silane  (note  the  much  larger  scale  for  the  CF4  case)  and  values  up  to  i  =  20  are 
still  contributing  in  a  sizeable  way  around  the  F  atoms. 

The  specific  behavior  of  the  electronic  densities  is  obviously  to  influence  the  features 
of  the  different  contributions  to  the  electron  molecule  interaction  (4.2.3),  since  they  are 
going  to  depend  rather  closely  on  the  electronic  density  of  the  target.  In  particular,  the 
use  of  expansion  (5.2.1-2)  for  the  bound  and  continuum  electrons  allows  us  to  write  the 
CC  equations  as  coupled,  radial  integro-differential  equations  for  each  FN  configuration 
(see  Eqs.  4.1.5  and  4.2.2).  The  coupling  matrix  of  Eq.  (4.1.5)  is  given  by: 


VZ&'I-W  =  WWir.fHrJTOtr)).  (5.2.3) 

If  we  write  the  SCE  expansion  of  interaction  potential  in  terms  of  SAWF: 


=  (5-2.4) 

hi 
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where  we  have  assumed  a  closed-shell  molecule  with  A\  IR  for  its  total  density  and  inter¬ 
action  potential,  we  have  to  go  through  a  numerical  quadrature  in  order  to  evaluate  the 
multipolar  expansion  of  the  potential.  The  results  for  SiH4,  GeH4  and  CF4  targets  were 
shown  in  Figs.  5.2.2,  5.2.3  and  5.2.4  of  Ref.  10,  respectively. 

From  Fig.  5.2.2  of  Ref.  10,  we  see  for  the  silane  molecule  that  the  radial  region  of 
the  coefficients  around  the  second  nuclear  cusp  on  the  H  atom.  The  Si  cusp  converges 
very  rapidly,  as  expected,  and  the  static  potential  at  r  ~  0  has  its  largest  contribution 
given  by  the  £  —  0  coefficient.  It  is  clear  from  seeing  the  increasingly  narrower  cusps  as  £ 
increases,  that  several  terms  are,  however,  still  needed  to  correctly  describe  quantitatively 
the  localized  anisotropy  of  interaction  due  to  the  four  H  atoms,  placed  at  the  four  comers 
of  the  usual  cube,  describing  tetrahedral  molecules.  On  the  other  hand,  the  values  of  the 
various  coefficients  at  R  ~  2.787ao  (  the  position  of  the  H  atom)  decreases  rather  rapidly 
with  increasing  £  and,  therefore,  one  sees  clearly  in  the  bottom  part  of  the  Fig.  5.2.2  of 
Ref.  10  that  by  the  time  £  =  20  the  contribution  to  the  static  potential  is  less  than  0.2 
atomic  units  and  it  is  so  over  rather  limited  radial  region  (AR  ~  0.3ao).  As  this  value 
should  be  compared  with  the  values  of  lower  coefficients  over  the  whole  range  of  action  of 
the  static  potential,  it  is  reasonable  to  expect  that  such  high  multipoles,  in  the  present 
case,  would  contribute  rather  small  to  the  scattering  process. 

The  similar  coefficients  for  GeH4  molecule  were  depicted  in  Fig.  5.2.3  of  Ref.  10. 
They  show  even  more  dramatically  the  effect  of  the  compact  charge  distribution  around 
the  central  Ge  atom  and  the  more  removed  position  of  the  H  atoms  (R  =  2.905a0);  the 
atomic  charge  at  the  COM  of  the  system  has  increased  here  from  14  to  32  and,  therefore, 
the  corresponding  static  potential  with  an  incoming  electron  is  given  almcst  totally  by 
the  £  =  0  term.  The  coefficients  for  higher  £  values  are,  therefore,  nearly  negligible  both 
around  R  ~  0  and  all  along  the  bond  region  until  one  reaches  the  H  atom  environment, 
where  they  contribute  more  or  less  the  same  amounts  as  in  the  silane  case  of  Fig.  5.2.2 
of  Ref.  10.  The  higher  order  contributions  with  £  >  12  as  shown  in  Fig.  5.2.3  (see  Ref. 
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10),  indicate  once  more  the  fairly  rapid  convergence  of  the  static  potential  expansion  for 
the  GeH4  target;  nearly  all  coefficients  vanish  within  the  examined  radial  region,  except 
in  the  narrow  area  around  ~  2.9a0,  where  they  never  contribute  more  than  ~  1  au.  All 
the  higher  contributions  with  £  >  20  and  up  to  £  ~  30  are  markedly  smaller  and  were  not 
shown  in  the  figure. 

Thus  we  preliminarily  concluded  that  both  silane  and  germane  are  nearly  spherical 
molecules  with  compact  electronic  charge  distribution  which  are  given  quite  realistically 
by  the  first  few  coefficients  in  the  multipolar  expansion  of  Eq.  5.2.4. 

Finally,  for  the  CF4  case  (Fig.  5.2.4  of  Ref.  10),  The  £  —  0  term  is  still  the  largest 
around  the  origin  and  higher  order  terms  are  markedly  smaller  near  the  origin.  On  the 
other  hand,  the  radial  region  where  the  F  atoms  are  located  (i?  =  2.45a0)  show  now  much 
larger  contributions  from  the  higher  multipolar  terms:  the  scale  used  in  Fig.  5.2.4  of  Ref. 
10  is  one  order  of  magnitude  bigger  than  that  of  the  previous  cases  and  yet  we  see  that 
£  .=  10  is  still  contributing  rather  substantially.  The  same  is  true  for,  say,  i  =  19  that  is  now 
about  10  times  larger  than  its  value  for  H  atoms  of  the  previous  cases  (silane  and  germane 
molecules)  and  still  20  %  of  the  cusp  value  for  the  £  =  10  contribution.  One,  therefore, 
sees  here  a  nearly  linear  scaling  of  the  strength  of  the  static  potential  with  Z  value  of  the 
off-center  atom  that  is  being  considered.  Such  a  result  could  be  used  to  estimate  the  rate 
of  convergence  of  the  coefficients  and  the  £ma,x  value  which  should  be  included  in  (5.2.4)  to 
reach  any  predetermined  level  of  accuracy  in  the  static  potential.  In  the  present  example 
of  CF4  one  clearly  sees  that  higher  order  terms  need  to  be  included  and  £  values  up  to 
£  =  30  had  to  be  calculated  in  order  to  reach  the  same  negligible  cusp  behavior  around  F 
atoms  as  those  shown  before  by  the  other  two  molecules  discussed  above. 

Thus  with  the  above  SCE  codes,  we  have  made  preliminary  attempts  to  develop  SCE 
quantities  for  much  heavier  targets  (SF6,  CC14  etc.  )  including  the  CF4  one.  We  will  show 
some  of  our  cross  section  calculations  for  the  electron-SFg  system  at  low  to  high  impact 
energies. 
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5.3  An  Optimized  Iterative  Scheme 


The  full  account  of  this  section  is  already  given  in  the  previos  progress  reports  [9- 
10].  At  the  static-exchange  level  [1],  the  Schrodinger  equation  for  the  continuum  electron 
function  Fpli( r)  for  any  irreducible  representation  (pp)  in  the  BF  coordinate  system  can 
be  written  as  (see  also  Section  4.1): 

+  Kt(r)  -  ifc2]F<">(r)  =  W  «r')|r  -  r'fV <<->(r')*V.(r),  (5.3.1) 

a  ^ 

where  k  is  the  incident  electron  wavevector  and  the  static  potential  14t  is  given  by  Eq. 
(5-2.4).  Here  4>0  is  the  target  ground  state  wavef unction  given  as  a  single  Slater  determinant 
of  one-electron  N  spin  orbitals  <f>i(r)  and  M  is  the  number  of  nuclei  in  the  molecule.  The 
right  hand  side  of  Eq.  (5.3.1)  is  the  non-local  exchange-correlation  term  arising  from 
antisymmetrization  of  the  e-molecule  wave  function.  In  the  SME  approach,  this  term 
is  replaced  by  a  local  exchange  potential  and  the  corresponding  scattering  calculation 
becomes  quite  easy  computationally  (see  Ref.  1).  The  function  Mr)  in  the  following  will 
represent  the  radial  part  of  the  continuum  function  with  channel  indices  i  and  j. 

Finally,  after  projecting  the  integro-differential  equation  (5.3.1)  onto  the  symmetry- 
adapted  angular  basis  functions  of  F^pp\t)  (5.2.2),  we  obtain  a  set  of  coupled  integro- 
differential  equations  (4.2.2)  which  cam  be  written  in  a  convenient  matrix  form,  = 

WF^,  where  WF^^  is  the  exchange  term.  The  iterative  scheme  is  LF^P#1^  =  WF^*, 
with  i  =  1,2,...  .  The  optimized  scheme  is  summarized  as  follows. 

1.  We  need  to  solve  the  full  integro-differentinal  Eq.  (4.2.2)  only  for  a  few  low  partial 
waves;  We  have  earlier  shown  [6]  this  scheme  of  low-,  intermediate-  and  high-£  values 
which  implies  that  for  low  partied  waves  (say  £ma,x  =  Le ),  we  need  the  full  SEE 
scheme,  for  £  >  Le  (say  fmax  =  Lm),  we  can  employ  the  SME  model  and  finally  for 
£  >  Lm  (say  £max  =  Lg),  the  Unitarised-Bom-Approximation  (UBA)  [35]  will  be  a 
good  approximation.  This  partition  is  similar  to  the  angular-frame-transformation 
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(AFT)  scheme  of  Norcross  and  Collins  [31],  but  differs  in  implementation  in  an  actual 
calculation. 

2.  Suppose  we  want  to  determine  the  SEE  K-matrix  at  energy  E\.  When  FqME(£j)  is 
used  to  start  the  iterative  scheme,  it  leads  to  the  final  F|ee(£i)  converged  continuum 
function  after  n  iterations.  However,  if  we  want  to  apply  the  same  iterative  scheme 
at  another  energy,  say  £2,  we  cam  use  either  FqME(£2)  (requiring  approximately  the 
same  number  of  n  iterations)  or  F^ee(£j  )  (requiring  only  n  iterations)  to  obtain  the 
final  F|ee(£2)  at  energy  £2-  We  have  shown  [9]  that  n  is  much  less  than  n  and 
the  final  results  converge  to  the  same  values,  i.e.,  F^ME(£2)  =  F|ee(£2).  If  F(£) 
is  a  slowly  varying  function  with  respect  to  the  impact  energy  £,  we  can  use  a  well 
converged  scattering  function  at  an  energy  E\  to  start  the  iteration  at  the  next  energy 
£2 .  Thus,  by  using  the  criterion: 

PpE(£„)  =  BtEE(£,-1).  (5.3.2) 

to  start  iteration  at  a  given  energy  Ep,  we  can  save  a  significant  amount  of  computer 
time  without  losing  any  accuracy. 

3.  Also,  the  number  of  required  iterations  depends  almost  linearly  on  the  size  of  the 
K-matrix.  For  example,  for  the  electron-H20  scattering  in  the  Aj  symmetry,  a 
16-channel  problem  requires  34  iterations  as  compared  to  only  24  in  a  12-channel 
calculation  [9].  With  the  new  optimized  scheme  (Eq.  5.2),  as  described  above,  only  a 
much  smaller  number  of  iterations  are  required  to  obtain  desired  accuracy. 
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5.4  Computer  Codes 


In  order  to  carry  out  the  present  calculations,  we  used  the  following  computer  codes 

developed  by  our  group  for  more  than  the  last  15  years.  All  these  programs  are  now  set  up 

on  the  VAX,  CRAY-YMP  and  CRAY-2  machines  at  the  FSU  supercomputer  facilities. 

1.  MOLMON  generates  SCE  target  wavefunctions  in  terms  of  STO’s.  At  present,  this 
code  can  include  £miX  =  7  in  the  SCE  of  each  orbital.  This  code  is  made  available  to 
us  by  Professor  F.  A.  Gianturco. 

2.  MOPE  calculates  expansion  coefficient  of  the  single  center  expansion  of  static  po¬ 
tential  and  electronic  density  from  the  orbitals  generated  by  the  MOLMON  program. 
Dr.  D.  G.  Thompson  wrote  this  code. 

3.  FEGE  determines  expansion  coefficients  of  the  SCE  of  FEGE  or  correlation  polar¬ 
ization  potential.  This  code  was  written  by  S.  S.  Salvini. 

4.  ASYMP  provides  rotational  eigenenergies  and  eigenfunctions  of  a  symmetric  top 
molecule.  The  code  was  written  by  Ashok  Jain,  the  PI. 

5.  POPLE  calculates  the  polarization  potential  (see  Sec.  4.3)  by  employing  the  Pople- 
Schoffield  method.  The  program  was  written  by  Ashok  Jain,  the  PI. 

6.  POLY  is  the  main  scattering  code  written  by  N.  Chandra,  modified  for  polyatomic 
case  by  D  G  Thompson  and  later  modified  by  Ashok  Jain.  This  code  calculates  scat¬ 
tering  K-matrix,  eigenphase  sums  and  partial  cross  sections  for  a  particular  scattering 
state. 

7.  POLYEX  is  a  modified  version  of  POLY  to  implement  the  iterative  scheme.  This 
code  was  developed  by  Paul  McNaughten.  Later,  we  modified  this  code  to  implement 
optimised  iterative  scheme  and  also  interfacing  it  with  the  HONDO  code. 

8.  EROTVIB  was  written  by  Ashok  Jain.  It  evaluates  various  cross  sections  for  asym¬ 
metric,  symmetric  and  spherical  top  molecules  from  the  K-matrix  input  of  POLYEX 
code. 
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9.  RESON  is  a  small  code  which  determines  resonance  parameters  from  the  given  eigen- 
phase  sums. 

10.  INTERFACE  is  an  interfacing  program  between  our  POLYEX  code  and  the  Quan¬ 
tum  Chemistry  HONDO  program.  INTERFACE  calculates  SCE  quantities  from  mul¬ 
ticenter  electronic  density,  orbited  and  static  potential  quantities.  The  program  was 
written  by  Vi  D.  Martino  and  Ashok  Jain. 

11.  BORN  is  a  general  program  which  determines  FBA  quantities  in  the  body-fixed 
frame  of  references  (Bom  K-matrices,  Unitarised  Bom  T-matrix  etc.)  and  cross 
sections  (differential,  integral  and  momentum  transfer)  in  the  space-fixed  coordinate 
system.  This  program  has  the  provision  to  use  numerical  potential  (rather  than  the 
asymptotic  form  only).  This  code  is  written  by  Ashok  Jain. 

12.  POLDEN  determines  the  polarized  charge  density  of  the  target  by  employing  the 
Pople  and  Shofeld  method.  This  code  was  developed  by  Ashok  Jain. 

13.  VIBAVEG  calculates  vibrational  matrix  element  from  the  given  R-dependent  scat¬ 
tering  parameter  (amplitude,  potential  matrix,  multipole  moments  etc.  )  and  vibra¬ 
tional  eigenfunctions. 

14.  DERTV  is  a  small  program  which  determines  first  or  second  derivative  of  a  given 
R-dependent  quantity  (energy,  multipole  moments,  cross  sections  etc.  ). 
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6.  RESULTS  AND  DISCUSSION 

6.1  Electron-CH4  Cross  Sections 

The  total  elastic  cross  sections  (DCS,  <rt,  ffm)  on  the  e-CH4  were  published  [3-4]  before 
this  contract.  However,  in  this  project,  we  presented  results  on  the  rotationally  elastic, 
inelastic  and  summed  DCS  integral,  momentum  transfer,  and  energy-loss  cross  sections 
[see  Tables  6. 1-6.2  of  Ref.  9].  These  results  were  compared  with  recent  measurements  on 
the  DCS,  integral,  and  momentum  transfer  parameters.  We  also  carried  out  calculations 
in  the  Ramsauer-Townsend  (RT)  minimum  region,  in  the  7-8  shape-resonance  regime 
and  also  at  very  low  energies  (below  0.1  eV).  No  any  other  ab  initio  calculation  has  been 
extended  to  such  a  wide  range  of  energy  and  several  types  of  cross  sections. 

A  large  set  of  differential,  integral  and  momentum  transfer  cross  sections  for  any 
J  — *  J'  transition  (with  any  value  of  initial  J  value  to  a  final  allowed  J'  value)  are 
available  in  tabular  form  from  the  PI.  The  energy  range  is  0.001-20  eV.  In  addition,  we 
have  generated  single  center  target  basis  functions  for  a  large  value  of  £max  =  20.  In  the 
reported  results  (Tables  6.1  and  6.2  of  Ref.  9),  we  employed  our  earlier  SCE  wave  function 
with  only  ^max  =  7.  The  corresponding  surfaces  for  electronic  density,  static  potential  and 
individual  bound  orbitals  are  available  on  a  radial  mesh  from  the  author  of  this  report. 
The  CH4  molecule  has  always  served  a  prototype  molecule  for  testing  a  theoretical  model 
for  polyatomic  targets.  Our  ab  initio  calculations  with  exact  exchange  and  parameter- 
free  correlation-polarization  potential  have  shown  that  observed  behavior  of  DCS  can  be 
described  very  well  by  theory  even  at  very  low  energies. 

The  cross  sections  reported  in  Tables  6.1  and  6.2  of  Ref.  9  were  determined  from  the 
theory  described  in  the  previous  sections.  In  brief,  exchange  correlation  is  implemented 
exactly  and  polarization  effects  are  included  approximately  via  a  local  real  potential  based 
on  the  method  of  Pople  and  Shofield  [27]  and  the  non-penetrating  criterion  of  Temkin  [28], 
in  which  the  incoming  electron  is  not  allowed  to  penetrate  the  target  charge  cloud.  In  this 
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non-empirical  calculation,  our  results  at  very  low  energies  (RT  minimum  region)  compare 
very  well  with  experiment.  In  addition,  our  value  of  —2.9  au  for  the  scattering  length  for 
the  e-CH4  system  compares  well  with  the  experimental  value  of  —2.48  [36]. 


6.2  Electron-SiH4  Cross  Sections 

We  determined  the  low  energy  electron  collisional  cross  sections  from  silane  molecules 
in  the  exact-static-exchange  plus  polarization  (ESEP)  model  [9].  For  comparison  purpose, 
we  also  included  static-model-exchange  plus  polarization  (SMEP)  results;  this  comparison 
gave  us  an  idea  of  the  importance  of  exchange  effects  in  this  low  energy  region.  The  results 
without  polarization  (ESE)  were  also  included  in  the  discussion  [9], 

The  electron-silane  cross  sections  are  very  important  from  the  point  of  view  of  plasma 
physics  [37].  There  are  several  experimental  [37-46]  and  theoretical  studies  [47-52]  on  the 
low  energy  electron  interactions  with  SiH4.  Similar  to  CH4  case,  the  e-SiH4  low  energy 
cross  sections  are  characterised  by  RT  minimum  below  1  eV  and  a  shape  resonance  feature 
around  2-4  eV.  Very  recently,  low  energy  DCS  have  been  measured  by  Tanaka  et  al.  [46]. 

We  discussed  [9]  the  qualitative  features  of  low  energy  e-SiH4  cross  sections  such  as  the 
RT  minimum  below  1  eV  and  a  3-4  eV  shape  resonance  phenomenon.  Our  ESEP  model 
was  found  successful  for  (nearly)  the  correct  position  of  the  minimum  (around  0.25  eV  in 
a  recent  swarm  study  of  Kurachi  and  Nakamura,  Ref.  44  )  because  of  the  Aj  scattering 
symmetry.  Note  that  our  earlier  SMEP  (static,  FEGE  and  the  JT  polarisation  potentials) 
results  (Jain  and  Thompson,  Ref.  47)  or  the  present  ESE  (without  polarization)  curves 
are  incorrect  in  this  low  energy  region.  The  value  of  the  position  of  this  minimum  as 
obtained  by  Yuan  [51]  in  the  spherical  approximation  is  quite  low  (around  0.11  eV).  Other 
calculations,  where  exchange  is  treated  via  a  model  local  potential  [48-49],  do  not  give 
reliable  results  in  this  rather  low  energy  region.  The  ESEP  model  clearly  exihibited  a 
shape  resonance  phenomenon  around  3.75  eV;  this  structure  is  due  to  the  d-wave  of  the 
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incoming  electron’s  partial  waves.  Therefore,  the  E  symmetry  (which  also  have  a  d-wave 
component)  also  exihibits  a  weak  shape  resonance  feature  at  about  4  eV. 

We  calculated  [9]  rotationally  summed  (vibrationally  and  electronically  elastic)  differ¬ 
ential  cross  sections  in  the  energy  range  from  1.8  to  20  eV,  where  recent  absolute  measured 
values  were  available  [46].  This  set  of  data  and  their  discussion  gave  us  the  idea  of  the 
importance  of  polarization  effects  in  this  energy  range  for  various  angular  regions.  We 
found  that  the  polarization  interaction  changes  the  shape  (and  magnitude  too)  of  the  DCS 
at  all  energies  considered  here.  It  was  seen  [9]  that  the  ESEP  model  agrees  much  better 
with  the  experimental  data  as  compared  to  the  ESE  model.  We  also  noticed  that  the 
polarization  effects  are  more  important  in  the  resonance  region,  where  the  structure  in 
the  DCS  reflects  the  d-wave  dominance.  This  is  expected  since  polarization  interaction 
influences  the  shape-resonance  region  considerably.  In  general,  the  agreement  with  exper¬ 
imental  data  was  found  reasonable.  The  dip  at  about  55°  observed  at  1.8  and  2.15  eV  by 
Tanaka  ei  al.  is  shifted  in  our  results  towards  lower  angle  by  about  10°  difference.  The 
second  dip  occuring  in  our  DCS  around  130°  at  1.8  and  2.15  eV,  seems  to  be  shifted  in 
the  same  fashion  as  the  lower  angle  dip  with  respect  to  experimental  data  (however,  the 
measured  values  are  not  available  above  130°  angle). 

At  somewhat  higher  energies  (E  >  3  eV),  the  two  dips  in  the  DCS  are  still  visible,  while 
such  a  structure  in  the  measured  angular  functions  is  quite  weak.  It  may  be  interesting 
to  compare  our  DCS  with  another  set  of  absolute  experimental  data  (not  available  right 
now).  The  present  predicted  DCS  values  may  be  quite  useful  for  normalizing  purpose  in 
future  measurements.  No  other  theoretical  calculations  were  available  on  the  DCS  for  the 
e-SiH4  collisions  in  this  energy  range.  We  presented  our  ESEP  DCS  values  in  tabular  form 
(see  Tables  6.2.1  and  6.2.2  in  ref.  9). 

The  integral  cross  sections  (total,  at\  momentum  transfer,  erm;  viscosity,  av  and 
energy-loss,  07  )  were  calculated  in  the  RT  and  resonance  region.  Our  value  of  the  cross 
section  near  the  RT  minimum  is  about  2.8  x  10-16  cm2  for  the  crt  and  0.4  x  10-16  cm2  for 
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the  am.  The  swarm  data  of  Kurachi  and  Nakamura  [44]  gives  a  minimum  value  of  about 
1.1  x  10~16  cm2  in  their  am  curve,  while  Hayashi  [45]  estimates  this  minimum  value  to  be 
0.5  x  10"1<5  cm2  .  Finally,  we  concluded  [9]  that  the  agreement  between  present  theory 
(ESEP  model)  and  experimental  results  is  rat  \er  encouraging  keeping  in  mind  that  this 
low  energy  region  is  quite  difficult  to  study  both  in  theory  and  experiment. 

We  compared  our  at  and  am  values  in  the  1-20  eV  region  with  the  experimental 
data  [38,  42].  There  is  significant  discrepancy  between  the  two  sets  of  measured  values. 
Our  present  ESEP  curve  seems  to  be  in  fair  agreement  with  recent  measurements  of  Wan 
et  al  [38].  The  position  of  the  shape  resonance  in  ESEP  theory  (around  3.75  eV)  and 
experiment  (around  3.1  eV)  differs  about  20%,  while  the  magnitude  of  the  cross  sections 
in  this  energy  region  also  has  significant  discrepancy;  for  example,  the  theoretical  (ESEP) 
peak  "'alue  of  the  crt  is  about  20%  higher  than  the  measured  value  of  Wan  et  al..  Note 
that  the  inclusion  of  polarization  effects  shifts  the  shape  resonance  position  from  about  6.5 
eV  to  about  3.75  eV  and  increases  its  magnitude  by  about  45%.  Thus,  the  importance  of 
polarization  interaction  in  this  energy  region  can  not  be  overlooked. 

For  the  am  values,  we  found  significant  discrepancy  bewteen  theory  and  experiment; 
however,  there  seems  to  be  agreement  in  the  general  trend  of  the  cross  sections  as  a  function 
of  impact  energy.  In  the  resonance  energy  region,  our  calculations  are  much  higher  than 
the  swarm  data.  The  accuracy  of  the  swarm  analysis  above  1  eV  may  be  questionable. 
We  also  realized  that  a  better  polarization  potential  may  be  required  in  the  present  energy 
region.  The  momentum  transfer  cross  sections  are  dominated  by  large  angle  scattering, 
therefore,  there  may  be  other  reasons  too  for  the  observed  discrepancy  between  present 
theory  and  swarm  results. 

At  very  low  energies,  the  energy  loss  of  the  incoming  projectile  because  of  rotational 
channel,  is  a  very  important  mechanism  in  e-molecule  collisions.  It  is  now  possible  to 
measure  indirectly  the  rotational  excitation  cross  sections  for  the  electron-molecule  system 
(Muller  et  al.  [53]).  With  this  possibility  in  mind,  we  discussed  [9]  rotationally  elastic  and 
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inelastic  cross  sections.  The  elastic  0  — *  0  channel  dominates  all  other  inelastic  processes 
at  all  angles  except  at  those  angles  where  a  deep  minimum  occurs  in  the  0  — >  0  DCS.  From 
qualitative  point  of  view,  the  general  shape  of  rotationally  inelastic  DCS  for  the  e-SiH4 
system  is  quite  similar  to  the  corresponding  CH4  curves  [3].  The  0  — >  4  transition  depends 
weakly  on  the  scattering  angle,  while  the  0  — *  3  one  has  considerable  angular  dependence. 
The  0  — *  6  cross  sections  are  generally  smaller;  however,  as  the  energy  is  increased,  they 
become  larger  than  the  other  transitions  at  certain  scattering  angles.  The  total  energy-loss 
DCS  are  almost  flat  with  respect  to  angular  variation.  It  is  also  clear  from  this  work  that  a 
spherical  approximation  may  not  be  an  appropriate  model  for  the  total  DCS  in  this  energy 
range. 

We  found  an  effect  of  polarization  potential  on  the  rotational  excitation  process.  Sym¬ 
metry  of  the  polarization  term  (there  exists  only  the  spherical,  £  =  0,  term)  for  the  T 
point  group,  the  first  Bom  approximation  (FBA)  will  give  zero  cross  section  from  the  po¬ 
larization  potential  term  only.  In  the  present  calculation,  we,  however,  saw  a  substantial 
difference  in  the  DCS  with  and  without  polarisation  effects.  We  demonstrated  [9]  that  a 
polarised  target  is  more  efficient  to  excite  its  rotational  modes  through  the  angular  mo¬ 
mentum  transfer  from  the  projectile.  We  will  see  later  the  same  comparison  for  integral 
cross  sections  at  all  energies  considered  here. 

We  also  discussed  the  integral  values  of  various  rotational  excitation  cross  sections.  As 
expected,  the  pure  elastic  process  has  all  the  main  features  (RT  minimum  and  the  3-4  eV 
shape-resonance)  of  the  total  scattering  process.  The  rotational  channels  are  small,  but 
exihibit  interesting  features  in  the  energy  distribution  of  their  cross  sections.  For  example, 
the  0  — ►  3  curves  have  enhancement  in  the  cross  sections  around  8  eV,  while  the  3-4  eV 
region  shows  a  weak  structure  (remember  there  is  a  3-4  eV  shape  resonance  here  in  the 
elastic  and  total  cross  sections).  The  0  — ►  4  cross  sections  are,  however,  enhanced  in  the 
3-4  eV  region,  it  is  interesting  to  notice  that  this  structure  is  produced  mainly  because  of 
polarisation  effect.  The  0  — »  6  cross  sections  were  also  characterised  by  an  enhancement 
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around  10  eV.  Here  also,  we  see  a  significant  change  in  the  inelastic  cross  sections  when 
polarisation  effects  are  included. 

Finally,  we  presented  integrated  energy-loss  cross  sections  in  all  the  three  (A J  = 
3, 4, 6)  transitions  with  and  without  polarisation  effects.  We  saw  clearly  that  the  inclusion 
of  polarisation  effect  changes  the  behavioir  of  the  energy-loss  cross  sections  with  respect 
to  the  projectile  energy.  At  higher  energies  (E  >  15  eV),  the  0  — *  6  transition  dominates 
over  other  lower  order  rotational  excitations.  The  0  — ►  4  cross  sections  are  always  higher 
than  the  0  — ♦  3  ones  except  at  lower  energies  (below  2  eV),  where  the  0  — »  3  transition 
dominates  over  all  higher  excitations.  The  total  energy-loss  or  the  stopping  cross  sections 
were  found  independent  of  initial  J  value  of  rotational  state  of  the  target.  This  agrees 
with  a  general  theorem  (per  Shimamura  [11-17])  that  the  stopping  cross  sections,  when 
summed  over  all  final  rotational  states,  are  independent  of  the  initial  rotational  state  of 
the  molecule  in  a  situation  where  the  adiabatic-nuclei-rotation  approximation  is  valid. 

We  also  extended  our  calculations  at  further  lower  energies  towards  the  zero-energy 
limit.  For  example,  at  0.001,  0.0025,  0.005,  0.01,  0.025,  and  0.05  eV  energies  the  values 
of  the  crt(am)  are  (in  units  of  1(T16  cm2  )  47.82  (46.45),  42.53  (40.4),  37.42  (34.64),  30.88 
(27.63),  20.68  (16.63),  and  12.83  (8.58)  respectively.  At  0.01  eV,  the  am  values  given  by 
variou.  swarm  studies  are  56.0  (Hayashi  [45]),  66.0  (Ohmori  et  al.  [43]),  and  32.0  (Kurachi 
and  Nakamura  [44]).  It  is  rather  frustating  to  see  such  a  large  discrepancy  between  these 
swarm  values  for  the  a,n  value  at  this  low  energy.  The  value  of  scattering  length  was 
estimated  to  be  —4.2  au  in  our  present  ESEP  calculation. 

It  is  now  a  well  established  fart  that  in  low  energy  regime,  the  use  of  model  potentials 
for  both  the  exchange  and  polarization  forces  is  rather  misleading  because  to  make  theory 
and  experiment  closer,  these  interactions  do  compensate  for  each  other.  In  the  situation 
where  exchange  is  included  exactly,  In  this  sub  section,  we  made  a  relative  comparison  with 
our  SEP(JT)  results  (previous  6.2  Sections)  with  another  model  SEP(CP)  calculation. 

Finally,  our  SEP(JT)  model  appears  to  be  better  than  the  SEP(CP)  one,  in  particular 
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in  the  RT  minimum  region.  This  conclusion  is  consistent  with  our  similar  investigation 
on  the  e-CH4  system  [3].  However,  this  conclusion  will  further  be  supported  if  differential 
measurements  are  available  for  the  e-SiH4  case.  There  is  plenty  of  room  to  improve  upon 
the  JT  potential  by  actually  carrying  out  a  full  polarised-orbital  calculation  for  the  e-SiH4 
system. 


6.3  Electron— GeH4  Cross  Sections 

This  is  a  rather  too  heavy  spherical  top  molecule  for  which  very  recently  DCS  have 
been  measured  in  the  range  of  2-100  eV  [54].  No  previous  calculation  exists  for  this  collision 
system.  Molecular  wave  functions  and  consequently  the  density  and  various  interaction 
potentials  were  determined  using  the  HONDO  and  INTERFACE  programs  (see  Section 
4.3).  This  part  of  the  calculation  was  done  at  a  simpler  level  than  the  one  described  above 
for  other  gases.  A  spherical  approximation  [19]  was  employed  under  the  model  exchange 
and  polarization  terms.  The  full  results  were  reported  [9]. 

We  presented  [9]  low  energy  (at  2,  3,  5,  7.5,  and  10  eV)  DCS  together  with  the 
experimental  data  [54].  The  validity  of  the  spherical  model  depends  on  the  fact  that 
the  rotational  excitation  channel  is  very  weak  compared  to  the  elastic  process.  This  is 
not  true  below  10  eV,  where  rotationally  inelastic  DCS  are  significant.  The  importance 
of  rotationally  inelastic  channel  is  increased  when  strong  interference  effects  make  the 
elastic  DCS  become  either  smaller  or  bigger  at  a  specific  angle  for  each  energy.  Hence  the 
ensuing  DCS  are  characterized  by  a  dip  in  the  angular  distribution.  We,  therefore,  expect 
some  significant  discrepancy  to  appear  between  our  theoretical  DCS  and  the  measured  ones 
below  10  eV  whenever  the  above  feature  may  be  present.  Another  source  of  approximations 
below  10  eV  DCS  is  the  treatment  of  exchange  and  polarization  effects  within  the  present 
model.  Nevertheless,  we  found  a  satisfactory  agreement  between  theory  and  observation. 
However,  the  spherical  model  provides  at  times  only  a  qualitative  picture  of  the  observed 
data  since  we  see,  for  example,  a  very  pronounced  dip  around  120°  at  2  eV  which  is  instead 
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very  shallow  in  the  experimental  data. 

Further,  we  found  that  at  lower  angles,  the  experimental  DCS  are  very  small  as 
compared  to  computed  values.  As  discussed  by  Boesten  and  Tanaka  [55],  the  lowest  angle 
DCS  are  rather  too  small  in  their  measurements.  It  is  also  interesting  to  note  that  the 
background  scattering  is  quite  strong  at  2  eV. 

We  also  displayed  results  [9]  at  intermediate  and  high  energies  where  our  present 
model  is  likely  to  work  markedly  better  [50].  We  found  that  there  is  indeed  very  good 
agreement  between  theory  and  experiment  and  that,  in  particular,  the  structure  observed 
in  the  experimental  DCS  is  clearly  reproduced  in  our  calculations  at  all  energies  above  10 
eV. 

Finally,  we  displayed  [9]  the  computed  integral  and  momentum  transfer  cross  sections. 
We  saw  a  shape-resonance  feature  around  5  eV  which  is  caused  by  the  d-wave  component 
of  our  potential  scattering  mechanism.  Note  that  in  CH4  and  SiH4  molecules  this  enhance¬ 
ment  in  the  total  cross  section  occurs  around  7-8  eV  and  2-3  eV  respectively.  Thus  the  5  eV 
structure  seen  here  in  the  e-GeH4  case  becomes  quite  interesting  in  comparison  with  other 
T^-symmetry  targets.  Below  1  eV,  we  have  also  evidence  of  a  minimum  around  0.2  eV, 
but  such  a  structure  is  very  sensitive  to  the  approximations  involved  in  including  exchange 
and  polarization  effects  and  therefore  we  expect  that  it  will  need  possible  confirmation 
from  more  extensive  calculations  that  are  presently  planned  in  our  laboratory. 

In  order  to  further  confirm  the  5  eV  feature  in  the  e-GeH4  elastic  scattering,  we 
plotted  DCS  as  a  function  of  energy  [9].  The  experimental  DCS  of  Boesten  and  Tanaka 
[54]  clearly  exihibit  this  enhancement  around  5  eV;  this  shape  resonance  phenomena  was 
also  present  in  our  calculations.  It  is  still  to  be  seen  if  measured  total  cross  sections  and 
more  accurate  calculations  reveal  this  structure. 

It  is  still  interesting  to  note,  however,  that  the  present  calculations,  which  effectively 
decouple  each  contributing  partial  wave  in  the  spatial  region  where  the  potential  exists,  are 
indeed  presenting  a  fairly  good  description  of  the  scattering  process.  Since  the  coherent 
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sum  of  angular  momenta  carried  out  to  yield  the  DCS  is  a  result  of  dynamical  interference 
between  trajectories,  we  are  able  to  say  that  the  present  spherical  potential  is  realistically 
treating  such  interference  effects.  The  lack  of  angular  anisotropy  in  the  chosen  interac¬ 
tion  proven  f  ”  us  from  correctly  including,  however,  the  potential  coupling  which  causes 
addition^  erference  in  the  multichannel  S-matrix  of  the  full  problem  [1].  However,  it 
appears  from  the  present  results  that,  at  least  above  10  eV,  such  further  effects  are  rather 
negligible  and  that  the  heavy  GeH4  target  can  be  considered  as  being  essentially  still  dur¬ 
ing  the  interaction  (rcou  <<  rrot)  and  is  approximately  spherical  when  probed  by  the 
impinging  electron.  Thus,  the  possible  adiabatic  convolution  of  the  scattering  amplitude 
over  rotational  eigenstates  can  be  approximately  treated  via  its  spherical  component  only: 
i.e.,  we  can  write, 

«,/ 

*,/ 

-  D*  E  )\^f) 
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*>/ 

~]T/oo(ft)%  (6.3.1) 
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where  Q  is  the  space-fixed  orientation,  r^+i  is  the  body-fixed  electron-molecule  orienta¬ 
tion  vector  and  the  | tpi)  is  rotational  eigenfunction  of  the  molecule.  Because  of  time  scale 
considerations,  vibrational  excitations  are  even  less  likely  to  occur. 

6.4  Electron-NH3  Cross  Sections 

There  are  recent  measurements  on  the  e-NH3  total  [55-56]  and  differential  [57]  cross 
sections  in  the  present  energy  region;  however,  no  ab  initio  calculation  (with  exact  exchange 
plus  polarization  effects)  exists  in  order  to  make  a  proper  comparison  with  the  experiment. 
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The  Schwinger  multichannel  variational  (SMV)  calculations  of  Pritchard  et  al.  [58]  are 
carried  out  only  at  the  fixed-nuclei  static-exchange  level;  consequently,  in  the  absence  of 
polarization  effects,  it  may  not  be  proper  to  compare  these  SMV  results  with  experimental 
cross  sections.  More  recently,  Parker  et  al.  [59]  have  employed  complex  Kohn  variational 
principle  to  study  e-NH3  problem  at  the  ab  initio  level.  The  only  other  fixed-nuclei  close¬ 
coupling  calculations  employing  model  exchange  and  polarizations  potentials  are  due  to 
Jain  and  Thompson  [60]  and  Giant urco  [61]. 

Encouraged  by  the  success  of  our  present  calculations  on  the  CH4,  SiH4,  and  GeH4 
molecules,  we  extended  our  ESEP  theory  (Section  4)  for  the  e-NH3  case  where  recent 
measured  data  are  available.  As  is  well  known,  for  a  polar  molecule,  the  ANA  theory 
predicts  infinite  cross  section  in  the  forward  direction.  Here  we  restricted  ourselves  only 
to  the  eigenphase  sum,  partial  cross  sections,  and  momentum  transfer  cross  sections.  We 
also  present  DCS  and  integral  parameters,  but  these  numbers  can  be  considered  only 
prelimenary  at  this  time. 

As  mentioned  above,  the  effects  of  polarization  force  is  significant  at  all  energies.  We 
presented  our  crt and  erm cross  sections  along  with  experimental  data.  We  saw  a  qualitative 
agreement  between  our  theory  and  measurements.  The  enhancement  in  both  the  inte¬ 
gral  and  momentum  transfer  cross  section  around  10  eV  is  faithfully  reproduced  by  our 
calculations  [9]. 

Finally,  we  reported  our  preliminary  angular  functions  at  8.5  and  15  eV  along  with 
measured  points  and  the  static-exchange  (without  polarization)  SMV  calculations  of 
McKoy  and  coworkers  [58].  Unfortunately,  the  experimental  DCS  on  the  low  energy  e- 
NH3  scattering  were  not  available  properly.  The  measured  values  of  Shyn  [57]  are  only 
preliminary  and  definitive  conclusion  about  the  agreement  between  theory  and  experiment 
is  very  hard  to  draw  at  this  moment. 
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6.5  Electron-H2 O  Cross  Sections 


We  calculated  converged  cross  sections  (DCS,  <7t,  <rm)  for  elastic  and  rotationally 
inelastic  processes  [10].  The  H2O  molecule  is  polar  molecule  with  above  critical  value  of 
the  dipole  moment.  Thus  the  dipole  moment  of  water  molecule  is  capable  of  supporting 
infinite  number  of  bound  states  of  the  incoming  electron.  We  have  modified  our  existing 
codes  in  order  to  include  full  contribution  of  the  higher  partial  waves  which  is  important  in 
the  forward  angle  scattering.  There  are  several  recent  measurements  on  the  total  [61-64] 
and  the  DCS  [65-69]  parameters  for  the  e-H20  system.  We  made  a  comparison  between 
our  theory  and  these  measured  data.  In  addition,  other  calculations  [70-74]  were  also 
included  in  the  present  comparison. 

First,  we  reported  our  values  of  eigenphase  sums  (Table  6.1.1  of  Ref.  10)  and  partial 
cross  sections  (Table  6.1.2  of  Ref.  10)  for  all  the  four  symmetries  Ai,  A2,  Bj,  and  B2.  In 
these  tables,  we  included  our  calculations  with  and  without  polarization  effects.  Again  we 
found  considerable  effect  of  polarization  force  in  the  e-H20  scattering.  In  this  case,  there 
is  no  resonance  or  RT  effects.  This  is  mainly  due  to  the  dominance  of  dipole  scattering. 
The  inclusion  of  polarization  force  in  the  ESE  results  is  to  decrease  the  cross  sections. 
We  saw  a  very  weak  hump  structure  in  the  ESEP  Ai  symmetry  results  around  6  eV. 
The  total  DCS  at  various  energies  (1,  2,  3,  5,  6,  7.5,  10  and  15  eV)  were  presented  along 
with  several  recent  experimental  data.  We  found  a  reasonable  agreement  between  present 
theory  and  experiment  at  all  energies  [10].  There  were  no  pronounced  structures  in  all  the 
DCS  presented  here  for  the  electron-H20  system. 

Further  we  presented  our  integrated  cross  sections  for  the  total  (<rt)  and  momentum 
transfer  (<7m)  quantities.  The  experimental  points  were  also  included.  Again  the  agreement 
between  theory  and  experiment  is  fair  considering  the  fact  that  the  measured  values  [65-69] 
have  been  obtained  by  extrapolating  their  DCS  at  lower  and  higher  angles;  this  procedure 
can  introduce  serious  error  in  the  integral  cross  sections.  All  these  cross  sections  were 
tabulated  [10]. 
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6.6  Electron— H2S  Cross— Sections 


The  total  cross  sections  for  the  electron-H2S  system  have  been  measured  by  several 
groups  [63,  75].  The  low  energy  differential  cross  sections  have  also  been  investigated 
in  the  laboratory  [76-77].  The  theoretical  calculations  on  the  DCS  and  integral  cross 
sections  have  been  carried  out  by  our  group  [78-79]  using  model  potentials  for  exchange 
and  polarization  effects.  The  Khon  complex  method  has  also  been  applied  to  study  low 
energy  electron  scattering  with  H2S  molecules  [80].  The  electron-H2S  collision  system  is 
known  to  have  a  shape  resonance  in  the  total  cross  section  behavior  [75-76].  It  was  shown 
by  us  [78]  that  this  shape  resonance  is  due  to  the  B2  symmetry  of  the  collision  system.  The 
occurrence  of  a  resonance  in  any  electron  molecule  event  is  very  sensitive  to  the  exchange 
and  polarization  effects.  In  other  words,  the  resonance  parameters  (width  and  position) 
have  a  strong  dependence  on  the  charge  correlation  effects.  We  found  here  that  the  B2 
symmetry  scattering  is  quite  sensitive  with  respect  to  the  treatment  of  exchange  effect. 

The  effect  of  polarization  is  dramatic  in  shifting  the  resonance  position  from  about 
5  eV  to  the  correct  position  around  2  eV.  The  width  is  changed  greatly  too.  The  model 
(exchange)  potential  calculation  [78]  differes  significantly  with  the  present  ESEP  model. 
The  exact  values  (in  units  of  eV)  of  resonance  parameters  (Er,  T)  in  these  model  are  as 
follows:  ESEP  (  1.82,  0.838),  ESE  (4.2,  2.5)  and  MSEP  (2.21,  1.3).  The  experimental 
values  of  the  position  of  this  shape  resonance  vary  between  1.8  to  2.5  eV. 

The  total  cross  section  for  the  electron-H2S  system  exhibit  two  peaking  regions  around 
2-3  eV  and  6-10  eV  [75].  The  structure  around  6-10  eV  is  definetely  due  to  Ai  symmetry 
as  is  evident  from  the  partied  cross  section  of  this  symmetry.  Our  theoretical  values  for 
differential  as  well  as  integral  quantities  produce  all  the  salient  features  in  the  cross  section 
behavior.  Very  recently,  Buckman  and  coworkers  [77]  have  measured  DCS  at  15  and  20  eV. 
Our  calculated  DCS  were  found  in  fair  agreement  with  the  measurements  of  Buckamn’s 
at  15  and  20  eV.  The  integral  cross  sections  were  reported  [10]  for  total  and  momentum 
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transfer  cross  sections.  For  the  total  cross  section,  we  found  a  fair  agreement  between  our 
results  and  the  measurements  of  Ref.  75.  The  first-Bom-approximation  (FBA)  results 
totally  fail  in  reproducing  the  shape  and  magnitude  of  the  cross-section  features.  This 
proves  the  simple  fact  that  first-order  theories  such  as  the  FBA  are  incapable  of  describing 
low  energy  electron  scattering  with  a  polar  molecule  such  as  the  H2S. 

In  Ref.  10,  we  have  given  values  the  DCS  and  integral  cross  sections  for  rotationally 
elastic,  inelastic  and  summed  processes.  We  have  given  rotationally  inelastic  results  only 
for  the  principal  quantum  numbers  (0  J')  (thus  summing  over  the  final  t'  quantum  numer). 
From  these  tables  on  the  DCS,  we  immediately  notice  that  the  0  — ♦  4  is  stronger  than  the 
0  — *  3  transition.  The  elastic  cross  sections  exhibit  a  dip  in  the  backward  direction  (around 
120°). 
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7.  ELECTRON  SCATTERING  FROM  C2H2  MOLECULES 


7.1  On  the  2IIg  Shape  Resonance 

The  electron  (e)  -C2H2  system  is  known  to  have  a  shape-resonance  in  its  2ng  symme¬ 
try  configuration)  in  the  2  eV  region.  There  is  considerable  discrepancy  about  the 

location  of  this  resonance  both  in  experimental  and  theoretical  studies  [81].  Total  scat¬ 
tering  [82-83]  ,  electron  transmission  [84]  and  vibrational  excitation  [85-86]  experiments 
have  observed  this  resonance  for  the  e-C2H2  scattering  around  2.4-2.6  eV.  Very  recently, 
Krumbach  et  al  [81]  have  employed  Feshbach  formalism  to  investigate  the  2IIg  state  of 
C2H^  and  found  the  position  of  this  resonance  around  2.9  eV.  Tossel  [87]  has  used  MS- 
Xa  method  to  investigate  the  e-C2H2  scattering.  No  other  theoretical  calculations  are 
available  in  the  literature  in  this  low  energy  (below  5  eV)  region  for  the  e-C2H2  system. 

In  a  shape-resonance  phenomenon,  the  projectile  electron,  with  non-zero  angular 
momentum,  is  trapped  by  its  own  centrifugal  barrier  to  form  the  temporary  molecular  ion 
(C2H^)  which  decays  via  the  quantum  mechanical  tunneling  of  the  extra  electron  through 
the  barrier.  The  position  (Er)  and  width  (T)  of  the  resonance  and  the  magnitude  of  the 
corresponding  cross  sections  are  highly  sensitive  to  the  treatment  of  non-local  short  range 
correlation  (electron  exchange  and  charge  distortion)  and  long  range  polarisation  effects. 

Two  well  known  local  approximations  for  exchange  effect  in  e-molecule  problems  are 
based  on  the  Hara  free-electron-gas-exchange  (FEGE)  [24]  and  semi-classical-exchange 
(SCE)  [88]  theories.  The  FEGE  model  has  also  been  employed  along  with  the  orthogo- 
nalization  procedure  (OFEGE,  see  Ref.  89)  where  the  continuum  function  is  forced  to  be 
orthogonal  with  the  bound  function  of  the  same  symmetry.  The  FEGE  approximation  gives 
incorrect  asymptotic  behaviour,  therefore,  its  asymptotically  adjusted  version  (AAFEGE) 
has  also  been  employed  in  the  literature  (see  Ref.  89).  Both  the  approximations  (FEGE 
and  SCE)  have  serious  drawbacks  when  applied  to  low  energy  e-molecule  case,  particularly 
when  the  scattering  event  is  dominated  by  resonance  effects.  The  success  of  these  local 
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potentials  depends  upon  how  well  they  reproduce  the  shape  of  the  differential  and  integral 
cross  sections  or  parameters  ( Er  and  T)  of  a  shape-resonance. 

Unfortunately,  the  FEGE,  AAFEGE  and  SCE  models  totally  fail  for  the  2ng  e- 
C2H2  scattering  (see  later  in  Fig.  7.1.1).  Here  we  employ  a  modified  version  of  SCE 
approximation  (to  be  denoted  as  MSCE),  recently  proposed  by  Gianturco  and  Scialla 
[90]  for  polyatomics.  We  also  time  the  FEGE  model  (to  be  denoted  as  TFEGE)  with 
respect  to  the  position  of  the  2ng  resonance  of  the  e-C2H2  scattering.  Since  the  electronic 
configuration  of  C2H2  molecule  does  not  have  a  Ttg  bound  orbital,  the  orthogonalisation 
techniques  is  not  required  for  the  present  case  of  2II9  symmetry. 

Full  details  of  the  present  method  can  be  found  else  where  [91].  In  brief,  we  solve 
the  following  set  of  coupled  differential  equations  (  after  carrying  out  all  the  single-center 
expansions  of  the  potential  and  scattered  function)  for  the  radial  components  of  the  con¬ 
tinuum  function  belonging  to  the  2IIg  symmetry, 
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where  the  potential  matrix  V^,*(r)  for  the  2ng  symmetry  is  given  by, 

^max  f 


A=0,2,.. 


where  C(£  1,^2, ^3 ',"*1,^2)  is  a  Clebsch-Gordan  coefficient.  The  potential  expansion  pa¬ 
rameter,  v\ (r)  is  a  sum  of  static  (vj),  exchange  (u”)  and  polarisation  (v£o1)  terms.  The 
calculation  of  target  charge  density  [p(r,  9)]  and  vsJ  is  done  ab  initio  [92-93]  by  employing 
molecular  wavefunctions  [94]  at  the  equilibrium  geometry  of  C2H2.  The  v^ol(r)  coefficients 
were  determined  in  the  correlation  polarisation  approximation  [95]. 

Both  the  FEGE  [24]  and  SCE  [88]  approximations  for  a  local  exchange  potential 
depend  on  p(r,  9)  and  the  incident  electron  energy  E.  The  SCE  potential  is  given  by, 


rSCE 


(r,0)  =  i[£  -  K.M)]  -  j[{£  -  V„(r,#)}2  +  4*p(r,«)]i,  (7.3) 
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The  main  assumption  in  the  SCE  approach  is  that  the  local  momentum  of  the  bound 
electrons  can  be  disregarded  with  respect  to  that  of  the  projectile  electron.  This  makes 
the  SCE  model  a  high  energy  approximation.  Very  recently,  Gianturco  and  Scialla  [90] 
have  introduced  corrections  in  the  SCE  theory  by  including  local  momentum  of  the  bound 
electrons,  which  can  play  a  significant  role  for  lower  collision  energies.  Thus  treating  the 
bound  electrons  within  the  isotropic  free-elect ron-gas  (FEG)  model,  the  MSCE  potential 
is  derived  to  be  [90], 

V«SCE(r,0)  =  \[ E  -  V„(r, 0)  +  a)]  -  |[{£  -  V.,(r, 0)  +  a)2  +  4xp(r,#)]i ,  (7.4) 

where  a  =  ■^j[3ir2p(r,  0)]^.  The  new  MSCE  potential  can  further  be  modified  by  realizing 
that  the  electronic  density  of  the  target  is  distorted  by  the  impinging  electron,  i.e.,  when 
polarisation  of  the  target  wavefunction  is  taken  into  account.  We  would,  therefore,  call  this 
potential  a s  MSCEP  when  in  the  MSCE  potential  (Eq.  7.4)  V,t  is  replaced  by  Vst  +  Vpoi- 
We  will  present  results  in  both  the  models,  i.e.,  MSCE  and  MSCEP. 

The  FEGE  potential  is  given  by  [24], 

V™°B(r,0)  =  hr(r,0)  i  +  ^la|i±5|  ,  (7.5) 

where  r\  —  (k2  +  2Ip  +  kj?)* /kp,  kp  is  defined  as  kp  =  [37r2p(r,  0)]  3  and  Ip  is  the  ionization 
potential  of  the  target.  In  the  limit  r  — *■  00,  the  local  momentum  kp  of  the  scattering 
electron  does  not  approach  the  correct  value  of  k.  Consequently  FEGE  with  Ip  =  0  is  the 
AAFEGE  potential  as  mentioned  earlier. 

We  solved  the  scattering  problem  (Eq.  7.1)  for  the  2ng  state  with  proper  convergence 
tests  on  the  K-matrix  elements  (including  the  eigenphase  sums  and  the  partial  cross  sec¬ 
tions).  This  included  the  number  of  terms  (Amax  =  20)  in  the  single-center  expansion  of 
the  interaction  potential,  the  size  of  the  close-coupling  K-matrix  (£max  =  20).  In  order 
to  include  the  contribution  of  the  long  range  quadrupole  and  polarisation  fields  in  all  the 
K-matrix  elements,  we  integrated  Eq.  (7.1)  up  to  a  radius  rmax  ~  10 £/fc,  where  k  is  the 
wave  number  of  the  incident  electron. 
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First,  in  Fig.  7.1.1,  we  have  shown  the  2IIg  symmetry  eigenphase  sums  (£a)  in  the 
SCE,  FEGE  and  AAFEGE  models  including  polarisation  effect.  The  SCE  and  AAFEGE 
models  do  not  exihibit  any  resonance  phenomenon,  while  the  FEGE  potential  (curve  C) 
gives  a  weak  structure  around  4  eV.  On  the  other  hand,  as  mentioned  earlier,  the  position  of 
the  e-C2H2  2ng  resonance  has  been  located  in  various  experimental  techniques  somewhere 
around  2.5  eV  [81]. 

In  Fig.  7.1.2,  we  have  shown  the  6a  quantity  in  MSCEP,  MSCE  and  TFEGE  (with  Ip  = 
5.25  eV)  models  including  polarisation  interaction.  We  see  that  the  MSCE  and  MSCEP 
models  reproduce  the  2IIJ  shape  resonance  in  the  vicinity  of  2  eV  in  good  agreement  with 
various  experimental  findings.  It  is  to  be  noted  that  there  is  no  tuning  involved  in  the 
MSCE  or  MSCEP  results  of  Fig.  7.1.2.  On  the  other  hand,  the  tuned  results  of  the 
TFEGE  approximation  (in  fact,  we  arbitrarily  chose  three  values,  5.0,  5.25  and  5.50  eV, 
of  Ip  and  found  that  the  I  =  5.25  eV  is  closest  to  the  2IIg  parameters)  also  predict  the 
features  of  the  resonance;  however,  a  tuning  procedure  is  always  inconvenient  and  requires 
prior  knowledge  of  cross  section  behaviour.  The  values  of  (Er,  T)  quantities  in  MSCEP, 
MSCE  and  TFEGE  models  are  respectively,  (2.20,  1.21),  (2.08,  1.06)  and  (2.48,  1.32). 
Thus,  we  conclude  that  the  MSCE  (MSCEP)  theory  of  Gianturco  and  Scialla  [90]  is  a  very 
good  local  approximation  in  order  to  represent  the  non-local  exchange  correlation  in  low 
energy  electron-C2H2  scattering.  We  would  emphasize  that  the  MSCEP  model  is  better 
than  the  MSCE  one.  We  also  found  that  it  is  the  d- wave  which  is  responsible  for  the  2IIg 
resonance  behaviour. 

Finally,  in  Fig.  7.1.3,  we  present  our  results  in  the  MSCEP,  MSCE  and  TFEGE 
models  without  including  polarisation  effects,  i.e.,  at  the  static-exchange  (SE)  level.  In 
general,  the  2  eV  shape-resonance  structure  in  other  isoelectronic  molecules  (CO  and 
HCN,  for  example)  occurs  at  a  higher  energy  if  polarisation  potential  is  neglected  in  the 
SE  calculation.  As  expected,  the  2  eV  resonance  (Fig.  7.1.2)  is  shifted  at  higher  energies 
in  all  the  three  MSCEP  (~  5  eV),  MSCE  (~  4.5  eV)  and  TFEGE  (~  5.5  eV)  models.  No 
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previous  static-exchange  level  results  are  available  for  the  present  C2H2  case.  Our  present 
MSCEP  and  MSCE  results  without  polarisation  effects,  shown  in  Fig.  7.1.3,  may  be  useful 
for  future  ab  initio  calculations  on  the  C2H2  molecule. 

7.2  Low  Energy  (0.01-20  eV)  e— C2H2  Crojs  Sections 

The  electron  (e“)  scattering  from  C2H2  molecules  plays  important  role  in  many  ap¬ 
plied  sciences  [96].  However,  there  is  a  paucity  of  theoretical  and  experimental  electron 
impact  cross  section  data  for  this  molecule  [97].  Earliest  measurements  on  the  relative 
differential  cross  sections  (DCS)  were  made  by  Hughes  and  McMillen  [98]  at  10-100  eV  in 
the  10°-150°  angular  range.  Recent  measurements  of  Kochem  et  al  [99]  on  the  elastic  DCS 
covered  a  low  energy  (0. 5-3.6  eV)  regime  from  5°  to  110°.  Very  recently  absolute  DCS  for 
the  electron  C2H2  system  have  been  measured  by  Khakoo  et  al  [96]  from  5  to  100  eV  in 
an  angular  range  of  10°-125°.  There  is  no  low  energy  theoretical  work  on  the  e"-C2H2 
system  except  the  intermediate  energy  (10-100  eV)  calculations  by  Mu-Tao  et  al  [100] 
and  of  Thirumalai  et  al  [101]  (only  at  10  eV).  Thus,  the  present  theoretical  study  is  the 
first  one  in  order  to  report  low  energy  vibrationally  and  electronically  elastic  (rotationally 
summed)  differential,  intergral  and  momentum  transfer  cross  sections  at  0.01-20  eV. 

The  total  cross  section  (<rt)  measurements  on  the  e--C2H2  system  are  available  due 
to  Briiche  [102]  (1-33  eV),  Hai-Xing  and  Moore  [83],  and  Sueoka  and  Mori  [82]  (1-400 
eV). 

In  our  (see  Section  7.1)  investigation  of  the  2ng  shape-resonance  phenomenon  in 
C2H2  due  to  low  energy  electron  impact,  we  found  that  this  collision  is  very  sensitive 
to  the  treatment  of  short-range  correlation  (electron  exchange  and  relaxation  of  bound 
orbitals  due  to  impinging  electron)  and  long-range  polarisation  effects.  In  the  cases  of 
several  features  in  the  energy  and  angular  dependence  of  the  cross  sections,  the  sensitivity 
of  the  cross  sections  with  respect  to  exchange  and  polarisation  forces  is  even  more  critical. 
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In  a  model  potential  approach,  where  non-local  nature  of  exchange  and  polarisation  effects 
is  approximated  by  local  forms,  it  is  very  difficult  to  describe  all  the  observed  features  in 
the  cross  sections  without  any  fitting  procedure. 

The  MSCEP  potential  alongwith  the  correlation  polarisation  term  was  found  to  be 
more  successful  than  the  other  models  (  see  Section  7.1).  Therefore,  we  will  present  our 
final  results  in  the  MSCEP  model  only. 

We  solved  the  scattering  problem  (Eq.  7.1)  for  various  symmetries  as  follows: 


Eg  :  A  =  0, 

*  =  0,2, 4,. ..20, 

Eu  :  A  =  0, 

*  =  1,3, 5,...  20, 

IIg  :  A  =  1, 

*  =  2,4, 6,  ...20, 

n„  :  A  —  1, 

*  =  1,3, 5, ...20, 

Ag:A  =  2, 

*  =  2,4, 6,... 20, 

Au  :  A  =  2, 

*  =  3, 5, 7,. ..20, 

*g:A  =  3, 

*  =  4,6, 8,...,  20 

$u  :  A  =  3, 

*  =  3,5, 7.  ..,20 

ConLribution  from  the  higher  order  terms  ( A  >  3)  is  very  small  and  therefore  neglected 
here.  We  carried  out  proper  convergence  tests  on  the  K-matrix  elements  (including  the 
eigenphase  sums  and  the  partial  cross  sections).  This  included  the  number  of  terms  (Amax  = 
20)  in  the  single-center  expansion  of  the  interaction  potential,  the  size  of  the  close-coupling 
K-matrix  (*max  =  20).  In  order  to  include  the  contribution  of  the  long  range  quadrupole 
and  polarisation  fields  in  all  the  K-matrix  elements,  we  integrated  Eq.  (7.1)  up  to  a  radius 
rmax  ~  10 */fc,  where  k  is  the  wave  number  of  the  incident  electron. 

In  terms  of  the  calculated  body-frame  K-matrix  for  a  particular  symmetry,  the  S  and 
T  matrices  are  defined  as 
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(7.6) 


SA  =(l+;KA)(l-iKA)-'  ;  Ta  =  1-Sa, 

The  integral  cross  sections  (<rt)  for  elastic  (rotationally  summed)  scattering  at  energy 
k2  is  given  by, 


(7.7) 


U>  A 


The  expression  for  the  DCS  is  obtained  in  a  standard  way  by  expanding  it  in  terms 
of  Legendre  polynomials, 


^(k\8)  =  Y.Al.PL(  cosO). 

Tj 

The  expansion  Al  coefficients  are  given  as  [89), 

V  “'Ivor  i  i  i\(o7j  i  i \L 


(7.8) 


Al  = 


4  fc2 


]T[(2Z  4- 1)(2/'  +  1)(2Z  +  l)(2f  +  l)]n,-i,(-*)'"''(-l)i 


iviv 


(o  0  o)(o  0  o)  £  (2i  +  l)(-l yW(ll'll'JL)Mj„Mli;  (7.9) 


where  the  M-matrix  is  constructed  as 


-£(-»>*  (I  -A  o)2*  ■ 


(7.10) 


From  equation  (7.8),  it  is  easy  to  derive  simple  form  for  the  momentum  transfer  (<rm) 
cross  sections  in  terms  of  Ai  coefficients,  namely, 


<rm  =  4tt[Aq  -  — -4i]. 


(7.11) 


Figs.  7.2.1-7.2.4  show  our  DCS  at  1,  2.0,  2.5  and  3  eV  respectively  along  with  the 
measurements  of  Kochem  et  al  [99] .  The  present  DCS  at  40°,  60°,  90°,  120°  and  150°  angles 
as  function  of  energy  (0-4  eV)  are  illustrated  in  Figs.  7.2.5-7.2.8  along  with  experimental 
data  of  Kochem  et  al  [99].  For  a  better  comparison  between  theory  and  experiment,  all 
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the  numbers  of  Kochem  et  al  [99],  shown  in  Figs.  7.2. 1-7.2. 8  are  multiplied  by  a  factor 
of  two.  First  of  all,  we  find  a  very  good  qualitative  agreemnet  between  our  theory  and 
measurements  of  Kochem  et  al  [99].  For  example,  the  structure  in  the  DCS  in  resonant 
region,  observed  in  the  measuremnets,  is  very  well  reproduced  in  the  present  parameter- 
free  model.  At  1  eV,  a  dip  around  40°  and  a  peak  around  90°  in  the  theoretical  DCS  is  also 
present  in  the  measured  data  of  Kochem  et  al  [99]  (Fig.  7.2.1).  Similaryly  at  2  eV,  there  is 
a  qualitative  agreement  between  theory  and  experiment.  At  2.5  and  3  eV  (Figs.  7.2.3  and 
7.2.4),  we  can  see  even  a  much  better  agreement  between  our  DCS  and  the  experimental 
results.  We,  however,  notice  that  for  this  comparison  between  theory  and  experiment,  the 
measured  cross  sections  are  multiplied  by  a  factor  of  two. 

We  make  a  further  comparison  between  our  DCS  and  experimental  DCS  of  Kochem 
et  al  as  a  function  of  energy,  where  measurements  are  available  at  few  angles  in  the  energy 
range  of  0.2-3.6  eV.  In  Figs.  7.2.5-7.2.8,  these  DCS  are  shown  at  40°,  60°,  90°,  120°  and 
150°  angles  in  the  0-4  eV  range.  Again,  we  see  a  good  qualitative  agreement  between 
present  calculations  and  observed  data.  As  expected,  the  DCS  around  2  eV  region  are 
significantly  influenced  by  the  2ng  shape  resonance  phenomenon  in  both  the  forward  and 
backward  directions.  A  low  energy  (around  0.2  eV)  minimum  (known  as  the  Ramsauer- 
Townsend  (RT)  effect)  is  also  present,  which  seems  to  be  supported  by  the  trend  in  the 
low  energy  experimental  points  (see  Figs.  7.2. 5-7. 2. 7).  Additional  experimental  studies 
below  1  eV  are  needed  at  this  time  in  order  to  confirm  the  existence  of  RT  effect  for  this 
molecule. 

We  now  compare  our  DCS  at  higher  energies  (5-20  eV)  with  a  very  recent  absolute 
experiment  by  Khakoo  et  al  [96],  In  fact,  we  provided  our  present  calculated  cross  sections 
to  this  experimental  group  (Khakoo  et  al)  before  they  obtained  their  final  results.  In  Figs. 
7.2.9-7.2.12,  we  have  plotted  our  DCS  at  5,  10,  15,  and  20  eV  respectively  along  with 
experimental  points  of  Khakoo  et  al.  At  10  eV,  we  have  also  included  theoretical  results  of 
Thirumalai  et  al  [101]  and  experimental  relative  DCS  of  Hughes  and  McMillen  [98].  From 
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Figs.  7.2.9-7.2.12,  we  conclude  that  our  model  potential  calculations  for  the  e“-C2H2 
system  describes  the  elastic  scattering  process  very  well  as  compared  to  a  most  recent 
experimental  findings.  The  structure  observed  in  the  DCS  at  all  these  energies  (5-20  eV) 
by  Khakoo  et  al  is  very  faithfully  reproduced  by  our  model.  For  example,  a  dip  structure 
around  90°  (5  eV),  80°  (10  eV)  and  75°  (15  and  20  eV)  is  present  in  both  the  theoretical 
and  experimental  investigations.  From  magnitude  point  of  view,  theoretical  results  are  a 
higher  (within  20%)  than  the  measured  values.  Other  previous  calculations  by  Mu-Tao 
et  al  [100]  are  very  poor  in  the  present  energy  region  and  therefore  not  shown  in  Figs. 
7.2.10-7.2.12. 

Fig.  7.2.13  displays  our  total  (vibrationally  and  electronically  elastic  and  rotationally 
summed)  cross  sections  in  a  wide  energy  range  0.01-20  eV.  There  is  low  energy  minimum 
around  0.2  eV  which  may  be  recognized  as  RT  effect.  The  shape  resonance  feature  around 
2  eV  and  a  broad  hump  around  7.0  eV,  observed  in  the  measurements  of  Sueoka  and  Mori 
[82]  and  also  in  the  old  experiment  of  Briiche  [102],  is  clearly  reproduced  in  our  calculations. 
In  fact  the  position  and  width  of  the  2ng  resonance,  calculated  in  the  present  non-empirical 
model,  agrees  very  well  with  previous  theoretical  and  experimental  results  (see  Sec.  7.1). 
The  disagreement  between  our  results  and  the  experimental  data  (Fig.  7.2.13)  with  respect 
to  the  magnitude  of  the  cross  section  around  2  eV  (our  calculated  results  are  higher  by 
about  50  %  than  the  experiment)  is  consistent  with  the  similar  situation  found  in  the 
e~-CO  scattering  [103].  One  reason  of  this  discrepancy  between  theory  and  experiment 
has  been  attributed  to  the  neglect  of  nuclear  dynamics  in  the  present  theory  (vibrational 
channel),  which  may  influence  elastic  scattering  significantly  in  this  resonant  region.  It 
is  also  to  be  noted  that  the  measurements  of  Sueoka  and  Mori  [82]  are  relative  and  do 
not  consider  corrections  due  to  forward  scattering.  Another  set  of  experimental  data  may 
further  resolve  this  discrepancy.  Unfortunately,  we  do  not  have  any  other  measuremnts  on 
the  abosolute  value  of  the  e~-C2H2  total  cross  section.  At  10  eV,  the  theoretical  crt  value 
of  26.85  x  10-16  cm2  of  Thirumalai  et  al  [101]  is  in  fair  agreement  with  the  present  value 
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of  24.83  x  10-16  cm2.  The  corresponding  theoretical  value  of  63.6  of  Lee  Mu-Tao  et  al 
[100]  does  not  seem  to  be  realistic. 

Finally,  in  Fig.  7.2.14,  we  have  plotted  our  am  cross  sections  in  the  same  energy  regime 
(0.01-20  eV),  where  the  shape  of  the  cross  section  is  very  similar  to  the  at  curve  of  Fig. 
7.2.13.  To  the  best  of  our  knowledge,  no  theoretical  or  experimental  data  is  available  in 
the  literature  for  comparison.  At  10  eV,  our  value  of  17.1  x  10-16cm2  compares  well  with 
the  18.9  x  10“16cm2  value  of  Thirumalai  et  al.  In  brief,  the  cm  curve  is  also  charaterized 
by  low  energy  (0.2  eV)  RT  minimum,  a  shape  resonance  around  2  eV,  and  a  broad  hump 
around  7  eV. 

We  presented  differential,  integral  and  momentum  transfer  cross  sections  for  the  elastic 
scattering  of  low  energy  electrons  from  C2H2.  No  fitting  procedure  is  involved  in  the 
present  model  potential  approach,  where  exchange  and  polarisation  effects  are  treated 
approximately.  The  numerical  method  is  based  on  the  fixed-nuclei,  single-center  and 
closed-coupling  formalism.  All  the  salient  features  in  the  cross  sections,  observed  in  various 
measurements  for  both  the  differential  and  total  cross  sections,  are  reproduced  very  well 
by  the  present  theoretical  model.  This  study  also  suggest  that  more  experimental  studies 
are  required  in  order  to  understand  some  discrepancy  between  theory  and  experiment.  In 
order  to  make  a  best  comparison  of  our  DCS  with  measurements  of  the  German  group 
(Kochem  et  al)  we  multiplied  their  numbers  by  a  factor  of  two;  no  simple  explanation  is 
available  for  this  discrepancy  in  the  magnitude  of  the  DCS.  their  values 
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Energy  (eV) 


!lUre  711  Eigenphase  sums  ftr  the  ng  e-C2H2  scattering  in 
(curve  A),  AAFEGE  (curve  (B)  and  FEGE  (curve  C)  models. 
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Figure  7.1.2  Eigenphase  sums  for  the  Ilg  e-C2H2  scattering  in 
MSCEP  (curve  A),  MSCE  (curve  (B)  and  TFEGE  (curve  C)  models. 
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Eigenphase  Sums  (Radians) 
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Figure  7.1.3  Eigenphase  sums  for  the  Tig  e-CgK^  scattering  in 

MSCEP  (curve  A),  MSCE  (curve  (B)  and  TFEGE  (curve  C)  models. 
Results  are  at  the  static-exchange  (without  polarization)  level. 
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Figure  7.2.1  DCS  for  the  e"-2h2  scattering  at  leV.  Solid  line, 

present  calculations;  Crosses  and  open  boxes  (multiplied  by  a 
factor  of  two),  experimental  data  of  Kochem  et  al  (1985). 
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Figure  7.2.2  DCS  for  the  e‘-2h2  scattering  at  2eV.  Solid  line, 

present  calculations;  Open  boxes  (multiplied  by  a  factor  of  two), 
experimental  data  of  Kochem  et  al  (1985). 
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Figure  7.2.3  DCS  for  the  e"-2h2  scattering  at  2.5  eV.  Solid  line, 

present  calculations;  Open  boxes  (multiplied  by  a  factor  of  two), 
experimental  data  of  Kochem  et  al  (1985). 
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Figure  7.2.4  DCS  for  the  e'~2h2  scattering  at  3  eV.  Solid  line, 

present  calculations;  Open  boxes  (multiplied  by  a  factor  of  two), 
experimental  data  of  Kochem  et  al  (1985). 


da/dQ  (10"16  cm2) 


Energy  (eV) 


Figure  7.2.5  Present  DCS  (solid  line)  as  a  function  of  energy  at  40° 
angle  are  compared  with  experimental  data  (multiplied  by  a  factor 
of  two)  of  Kochem  et  al  (1985). 
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Figure  7.2.6  Present  DCS  (solid  line)  as  a  function  of  energy  at  60° 
angle  are  compared  with  experimental  data  (multiplied  by  a  factor 
of  two)  of  Kochem  et  al  (1985). 
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Figure  7.2.7  Present  DCS  (solid  line)  as  a  function  of  energy  ji  90° 
angle  are  compared  with  experimental  data  (multiplied  by  a  fa' cor 
of  two)  of  Kochem  et  al  (1985). 
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Figure  7.2.8  Present  DCS  as  a  function  of  energy  at  angles  of  120° 
and  150°. 
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Figure  7.2.9  DCS  for  the  e-2h2  scattering  at  5  eV.  Solid  line, 
present  calculations;  Crosses,  experimental  data  of  Khakoo  et  al 
(1993) . 
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Figure  7.2.10  DCS  for  the  e-2h2  scattering  at  10  eV.  Solid  line, 
present  calculations;  Crosses,  experimental  data  of  Khakoo  et  al 
(1993) .  The  dash  curve  is  the  calculations  of  Thirumalai  et  al 
(1981),  while  open  circles  are  relative  measurements  of  Hughes 
and  McMillen  (1933). 
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Figure  7.2.13  Integral  cross  sections  for  the  e-2h2  collisions  in 
the  0.01-20  eV  range.  Present  calculations,  solid  line;  x, 
experimental  data  of  Sueoka  and  Mori  (1989);  o,  experimental 
data  of  Briiche  (1929). 
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Figure  7.2.14  Present  results  on  the  momentum  transfer  cross 
sections  for  the  e"-2h2  scattering  in  the  0.01-20  eV  energy  range. 


8.  EFFECT  OF  GAS  TEMPERATURE  ON  THE  CROSS  SECTIONS 

We  made  an  extensive  study  of  the  effect  of  gas  temperature  on  various  elastic  and  in¬ 
elastic  (particularly  the  rotational  excitation)  cross  sections  [9,  18].  At  a  given  temperature 
T,  a  typical  polyatomic  molecule  has  many  rotational  states  populated.  From  theoretical 
point  of  view,  it  is  easier  to  calculate  rotational  excitation  cross  sections  from  the  ground 
state  ( 7  =  0)  to  any  final  7  value  of  a  given  molecular  system.  In  many  practical  cases, 
the  transfer  of  rotational  angular  momentum  Jt  (j  J'  —  J|  <  Jt  <  J  +  7')  in  any  scattering 
event  involves  very  high  values  of  7  and  J'.  For  spherical  top  molecules,  the  differential 
cross  sections  (DCS)  for  rotational  excitation  from  any  initial  state  7  to  a  final  state  7' 
can  be  evaluated  as  follows  [18], 

(T  j> - 1-2  m  2J^J  (27+1+2A7)  k'  da  2 

-(J  J,k,6)~  ^  (27  -f  l)(2Jt  + 1)  kt  dD  °  J"k  t8-1) 


da 


where  k'  and  kt  are  kinematic  factors  of  the  wavevector  of  projectile  in  the  7  — >  J'  and 
0  — ►  Jt  cases  respectively  and  7'  =  7  +  A7.  In  the  summation  over  7t  in  Eq.  (8.1), 
usually  only  a  few  terms  are  needed  with  small  7t  values.  Thus,  given  a  set  of  cross 
sections  for  a  few  0  — *  Jt  transitions  only,  Eq.  (8.1)  can  easily  be  used  to  determine 
rotational  excitation  quantities  from  any  initial  7  value  including  very  high  ones.  This  is 
essential  when  one  needs  to  average  the  cross  section  over  the  distribution  of  rotational 
states  at  a  given  temperature.  Eq.  (8.1)  is  valid  under  the  adiabatic-nuclear-rotation 
(ANR)  approximation,  i.e.,  the  incident  and  scattered  electron  is  moving  much  faster  than 
the  rotational  motion  of  the  molecule.  In  the  present  case  of  spherical  top  molecules  and 
the  energy  range  of  0.5-20  eV,  the  ANR  approximation  is  valid  and  we  used  Eq.  (8.1)  to 
generate  various  cross  sections  with  different  initial  7  values  for  all  possible  final  7'  values 
allowed  by  the  Jt  selection  rule. 

If  one  neglects  the  kinematic  factors  (k'  and  kt)  from  Eq.  (8.1),  the  cross  sections 
between  the  rotational  states  7  and  J'  =  J  dh  A7  can  further  be  written  as  [14-15], 

=  (8.2) 
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with 


2 J±AJ  , 

-4(|AJ|)=  £  35(0  -  J,)/(2J,  +  1).  (8.3) 

J«=|AJ| 

where  we  assume  that  2J  —  |  A  J|  is  larger  than  any  that  contributes  appreciably  to 
the  sum  in  Eq.  (8.3).  In  other  words,  we  assume  that  dcr/d£l(0  —>  Jt)  is  almost  zero  for 
Jt  >  2J  —  |AJ|,  so  that  the  upper  limit  of  the  summation  in  Eq.  (8.3)  is  equivalent  to 
infinity.  Thus,  we  can  easily  write  |AJ|  in  the  argument  of  A  and  not  J  or  A  J. 

The  d<x/dQ( 0  — »  Jt)  can  be  expanded  as, 

(8.4) 

L 

where  Pl(cos6)  are  Legendre  polynomials  and  the  expansion  Al  coefficients  are  given  by 
[  see  Section  4.4], 


Al(0  -  Jt)  =  Y,  BC(^Li-  QQQ)C(£'L£'\ 000)(-l)J‘  W Jt L ) m M ft™  * . 

tv  it' 

(8.5) 

The  M  matrix  is  defined  as, 


mm'  hh'pn 


k’’ 


(8.6) 


and  B  =  [(2£  + 1)(2 £'  +  1)(2£  + 1)(2£/  +  1)]^^~^  (— i)*-*’ .  The  C-coefficients  are  the  usual 
Clebsch-Gordan  coefiecients  and  the  T-matrix  is  defined  in  terms  of  scattering  S  matrix 
as  T  =  (S  —  1)  where  S  =  (1  +  iK)/(l  —  iK)-1  is  written  in  terms  of  reaction  K  matrix  for 
each  symmetry  (p/z).  In  Eqs.  (8.5)-(8.6),  W(abcd;ef)  is  the  well  known  Racah  coefficient 
and  the  coefficients  are  expansion  terms  in  the  definition  of  symmetry  adapted  basis 
functions  in  terms  of  real  spherical  harmonics  [1].  For  the  case  of  total  and  momentum 
transfer  cross  sections,  Eq.  (8.5)  can  be  greatly  simplified  to  give  simple  expressions  [see 
Eqs.  4.4.9-4.4.10]. 

Here  we  have  investigated  the  behavior  of  various  scattering  parameters  on  the  initial 
rotational  state  J  and  temperature  of  the  gas  (CH4  and  SiH4).  This  is  important  because 
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in  rotational  equilibrium  at  room  temperature  or  higher,  the  gas  molecules  occupy  many 
different  rotational  states.  In  addition,  we  also  studied  the  dependence  of  the  rotationally 
elastic  and  inelastic  processes  on  the  initial  state  J. 

The  basic  scattering  parameter  is  the  differential  cross  section  as  a  function  of  scatter¬ 
ing  angle  (6)  and  energy  (Ar2)  for  any  rotational  transition  J  —>  J  +  AJ.  A  large  set  of  cross 
sections  can  be  obtained  for  state-to-state  rotational  transitions  and  rotationally  summed 
quantities;  however,  it  would  be  meaningful  to  derive  compact  information  from  this  huge 
amount  of  cross  section  data.  For  example,  the  differential  moments  (the  dependence  on 
6  and  k2  is  suppressed), 

A  <8-7> 

and  the  integral  moments  S(p;  J)  characterize  the  behavior  of  the  electron-molecule  sys¬ 
tem.  The  zeroth-order  integral  moment  with  /x  =  0,  S(0;  J),  is  the  usual  integral  (or  total) 
cross  section  (summed  at(  J)  or  state-to-state  crt(J  — ►  J')),  while  first-order  integral  mo¬ 
ment,  5(1;  J),  is  the  energy-loss  or  stopping  cross  section.  Another  useful  parameter  is  the 
momentum  transfer  cross  section  (summed  crm(J)  or  state-to-state  crm(J  — ►  J'))  which 
is  obtained  by  integrating  the  differential  cross  section  (Eq.  8.1)  with  the  weighted  factor 
of  (1  —  cosfl).  In  Eq.  (7.7),  Ej  =  BJ(J  +  1)  is  rotational  energy  of  the  molecule  in  Jth 
level  and  B  is  the  molecular  rotational  constant.  Higher  order  moments  (p  >  2)  in  Eq.  8.7 
provide  information  on  statistical  fluctuations  of  the  energy-loss  spectrum  from  its  mean 
value. 

In  practice,  all  observable  physical  quantities  are  connected  with  these  moments  av¬ 
eraged  over  the  rotational  distribution  function  N(J;T),  namely 

J 

where  N(  J;  T)  =  g(  J)e~Ej/*T  is  the  Maxwellian  distribution  function,  k  is  the  Boltzmann 
constant  and  g(J )  =  (2 J  +  l)2  is  the  statistical  weight  of  the  Jth  state  of  a  spherical  top. 
In  the  definition  of  N(  J;  T),  we  have  neglected  the  effect  of  nuclear  spins. 
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Another  interesting  aspect  of  the  physics  is  the  relationship  between  rotational  exci¬ 
tation  and  de-excitation  processes.  Shimamura  [14-15]  has  derived  partial  sum  rules  for 
the  J  — *  J  ±  A  J  rotational  transitions.  From  Eqs.  (8.2)-(8.3),  we  see  that  (denoting  the 
cross  section  quantity  in  general  by  a) 

o( J  — *  J  4-  AJ)  2J  ■+■  1  -f-  2AJ 

a(  J  — ►  J  —  A  J)  =  2  J  +  1  -  2A  J  ’  (8'9) 

and  the  sum  a(J  — ►  J  +  AJ)  +  a(J  —>  J  —  AJ)  is  a  constant  equal  to  2A(|AJ|).  How¬ 
ever,  if  we  include  the  effects  of  kinematic  factors,  particularly  at  the  lower  bound  of  the 
present  energy  region,  Eq.  (8.9)  may  be  true  only  approximately.  Shimamura  [14-15]  has 
also  discussed  the  relationship  between  weighted  rotational  excitation  and  de-excitation 
quantities,  i.e.,  AE+a( J  — *  J  +  AJ)  and  A E~a(J  — »  J  —  AJ).  Here  AE^  are  transition 
energies  given  by  jBAJ[±(2J  +  1)  +  AJ].  We  have  also  discussed  these  partial  stun  rules 
and  ratios  for  the  CH4  and  SiH4  cross  sections  [9,  18]. 

We  have  therefore  carried  out  a  systemetic  study  on  the  dependence  of  differential 
cross  section  (Eq.  8.1)  and  moments  (Eq.  8.7)  on  the  initial  rotational  state  J  and  of  the 
averaged  quantities  (using  Eq.  8.8)  on  the  gas  temperature  T.  We  compared  our  predictions 
with  the  theorems  as  discussed  by  Shimamura  [11-17]  regarding  the  J  or  T  dependence 
of  various  physical  quantities.  All  our  results  were  presented  only  for  the  spherical  top 
molecules  (CH4  and  SiH4).  The  basic  scattering  quantity,  i.e.,  da/dQ(0  — +  Jt;  k2,6),  is 
taken  from  our  recent  ab  initio  calculations  [see  Sections  6.1  and  6.2]  on  the  e-CH4  and 
SiH4  systems.  These  calculations  treat  electron-exchange  interaction  exactly  via  iterative 
procedure  and  the  charge  correlation  and  polarization  effects  are  included  approximately 
in  a  parameter-free  perturbative  approach  (see  Section  4.5). 

Our  previous  calculations  [4]  on  the  e-CH4  state-to-state  rotational  transitions  were 
in  qualitative  agreement  with  the  experimentally  extracted  values  of  Muller  et  al.  [53]. 
There  were,  however,  some  differences  between  theory  and  experiment  for  these  rotationally 
inelastic  channels  where  cross  sections  are  very  small  compared  to  the  pure  elastic  case. 
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For  example,  for  the  0  — *  3  transition,  calculated  values  show  zero  cross  section  in  the 
forward  direction,  while  experimental  points  exihibit  no  such  behaviour.  This  was  one  of 
the  purposes  of  this  study,  i.e.,  to  resolve  this  discrepancy  between  theory  and  experiment. 
For  the  CH4  and  SiEU  gases,  we  summarise  our  findings  [9,  18]  as  follows: 

1.  The  rotationally  summed  differential,  integral  and  momentum  transfer  cross  sections 
and  various  moments  are  independent  of  initial  rotational  state  J  and  hence  on  the 
gas  temperature  T.  The  higher  order  moments  depend  very  weakly  on  the  initial  J 
value. 

2.  we  found  that  the  partial-sum  rule  [14-15]  for  any  physical  quantity,  say  a , 

a(  J  -»  J  +  A  J)  4-  a(J  — »  J  —  A  J)  =  2A(|AJ|),  (8.10) 

is  approximately  true.  Here  A(|AJ|)  can  be  differential,  integral,  and  momentum 
transfer  cross  sections  or  the  moments  (Eq.  S.2-8.3). 

3.  There  is  a  significant  effect  of  rotational  averaging  on  the  rotational  excitation  cross 
sections.  The  pure  rotationally  elastic  channel  is  not  affected  with  rotational  avarag- 
ing.  In  general,  the  averaged  DCS  are  smaller  and  exihibit  less  structure  in  the  DCS 
as  compared  with  the  unaveraged  DCS.  (see  Ref.  9) 
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9.  ELECTRON-MOLECULE  SCATTERING  AT  HIGH  ENERGIES 

The  total  cross  sections  (crt)  (including  elastic  plus  all  energetically  possible  inelastic 
channels)  for  electron-molecule  systems  axe  important  in  many  applied  sciences  [104- 
107].  In  the  intermediate  and  high  energy  region  (roughly  from  ionization  threshold  up  to 
several  keV),  almost  all  inelastic  channels  (  rotational,  vibrational  and  electronic  excitation, 
ionization,  dissociative  processes  etc.  )  are  open.  In  this  energy  regime,  a  conventional 
close-coupling  theory  (see  Sec.  3)  of  electron-molecule  complex  is  an  arduous  task  and 
almost  impossible  to  carry  out  with  present  day  fast  supercomputers.  It  is,  therefore,  not 
surprising  that  almost  all  previous  calculations  on  the  crt  for  electron-molecule  systems 
have  been  restricted  to  low  energy  region  (typically  below  or  around  the  first  ionization 
threshold  )  only. 

Although  it  is  the  low-energy  electron  molecule  scattering  where  the  cross  sections 
exhibit  a  rich  structure  and  several  complex  phenomenon,  however,  a  knowledge  of  high- 
energy  o t  values  for  several  molecular  gases  is  required  [104-107].  The  importance  of  at 
in  radiation  physics  and  chemistry  has  been  discussed  in  the  literature  [108].  Such  cross 
sections  are  available  mostly  through  experimental  studies  and  theoretical  ax  values  in  the 
literature  are  scarce  (see  Table  7.1  of  Ref.  10),  even  for  the  simplest  H2  molecule. 

The  goal  of  this  calculation  was  the  following.  First,  to  suggest  a  general  and  quite 
simple  theoretical  method  which  can  predict  at  values  in  a  wide  energy  range  for  any  molec¬ 
ular  system  for  which  Hartree-Fock  level  wave  functions  are  available.  The  other  molecular 
quantities  needed  for  the  calculation  are  the  polarizability,  ionization  potential  and  various 
permanent  multipole  moments  (dipole,  quadrupole  etc.)  of  the  isolated  molecule;  second, 
to  provide  a  simple  fitted  formula  (see  Eq.  9.18)  from  our  calculated  data  for  the  total 
cross  section;  third,  to  provide  Bom-Bethe  parameters  [109-113]  for  all  the  molecules  listed 
above;  fourth,  to  emphasize  the  fact  that  individual  elastic  (  <xel)  and  inelastic  (crabs)  (or 
absorption)  cross  sections  (giving  rise  to  the  final  <7t  quantity)  should  compare  quite  well 
with  the  corresponding  experimental,  or  more  accurate  theoretical  results.  For  example,  in 
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the  present  case,  our  inelastic  (or  absorption)  cross  sections  should  be  an  upper  bound  to 
the  experimental  total  ionization,  electronic  excitation  and  dissociation  channels  together. 
The  u t  parameter  has  been  measured  for  a  large  number  of  molecules  from  low  to  high 
energies  [10]. 

Finally,  we  also  examined  correlation  between  molecular  and  scattering  parameters. 
A  quantitative  as  well  as  qualitative  correlation  picture  was  presented  [10]  which  depicts 
several  interesting  aspects  of  the  electron-molecule  scattering  at  high  energies.  In  partic¬ 
ular,  the  nuclear  charge  Z  or  the  occupation  number  of  the  target  plays  important  role 
in  the  values  of  ert  for  all  the  targets.  Such  correlations  would  allow  us  better  to  under¬ 
stand  the  scattering  mechanism  and  to  estimate  cross  sections  for  those  molecules  where 
experimental  or  theoretical  studies  are  difficult  to  perform. 

The  basic  philosophy  of  the  present  method  is  based  on  the  assumption  that  the 
non-spherical  nature  (providing  torque  to  the  molecule  for  rotational  excitation)  tf  the 
molecular  system  does  not  play  significant  role  in  shaping  up  the  total  cross  section  of  the 
high-energy  electron- molecule  collisions.  The  collision  time  is  too  short  and  rotational 
excitation  cross  sections  axe  insignificant  relative  to  elastic,  ionization,  etc.  processes.  In 
addition,  the  contribution  from  vibrational  excitation  process  is  assumed  to  be  negligble. 

We  first  assume  that  the  fixed-nuclei  approximation  is  valid  in  this  energy  region 
and  the  interaction  of  the  electron  molecule  system  can  be  represented  by  a  local  complex 
optical  potential,  namely: 

Vopt(r)  =  yR(r>+iyabs(r),  (9.1) 

where  the  real  part  is  a  sum  of  three  terms  : 

Vk(r)  =  yst(r)  +  lex(r)  +  Vpoir)-  (9.2) 

The  static  potential,  Vst(r),  is  calculated  from  the  unperturbed  target  wave  function  'I'o  at 
the  Hartree-Fock  level.  The  Vpo\  ( r )  represents  approximately  the  short-range  correlation 
and  long-range  polarization  effects,  while  the  Vex(r)  term  accounts  for  electron  exchange 
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interaction.  In  this  energy  region,  a  local  and  real  potential  model  for  exchange  and 
polarization  effects  is  adequate.  The  V^bs  in  Eq.  (9.1)  is  the  absorption  potential.  Due 
to  the  non-spherical  nature  of  a  molecule,  the  optical  potential  (Eq.  9.1)  is  not  isotropic. 
A  general  expression  for  Vropt(r)  for  any  target  can  be  written  in  terms  of  the  following 
multipole  expansion  around  the  center-of-mass  (COM)  of  the  molecule  [1]: 

^M  =  £v«(r)*<r>(r),  (9.3) 

th 

where  (pp)  denotes  the  ground-state  symmetry  of  molecule  and  the  symmetry  adapted  X 
functions  are  defined  in  terms  of  real  spherical  harmonics  33  W: 

= E  &!«?„,(*)•  (9.4) 

ru—O 

For  closed-shell  systems,  the  (pp)  is  the  totally  symmetric  1Ai  (non-linear  molecules) 
or  *£+  (linear  molecules)  irreducible  representation.  The  anisotropic  terms,  £  =  1,2, . .., 
in  expansion  (9.3),  provide  torque  to  excite  rotational  levels  in  the  molecule.  As  mentioned 
earlier,  our  main  assumption  in  this  work  is  that  such  higher-order  terms  are  weak  and 
can  be  treated  separately  in  the  first  Born  theory  and  added  incoherently  to  the  elastic 
part. 

First  we  determine  the  target  charge  density  p(r)  of  a  given  molecule: 

p(r)=  f  \V0\2dr1dr2...drz  =  2  ]T  |<£a(r)j2,  (9.5) 

J  a 

where  Z  is  the  number  of  electrons  in  the  target,  <£,  is  the  ith  molecular  orbital  and  a  factor 
of  2  appears  due  to  spin  integration  and  a  sum  being  over  each  doubly  occupied  orbital. 
It  can  be  shewn  that  for  closed-shell  molecules,  p(r)  belongs  to  totally  symmetric  one¬ 
dimensional  irreducible  representation  (*Aj  ,  or  1£+)  of  the  molecular  point  group 
(1).  All  the  four  potential  terms  (Fst,  Vex,  Vpoi  and  Vabs)  are  functions  of  p(r).  For  example: 

r  M 

Kt(r)  =  J  p(fi)|r  —  Jdri  -^Z.jr-R.l  *.  (9.6) 
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The  Vex  is  taken  in  the  Hara  free-electron-gas-exchange  (HFEGE)  model  [24]  and  Vpo\  is 
calculated  in  the  correlation-polarization  (COP)  approximation  [114-116]. 

The  imaginary  part  of  the  optical  potential,  Vabs(r),  is  the  absorption  potential  which 
represents  approximately  the  combined  effect  of  all  the  inelastic  channels.  An  ab  initio 
calculation  of  absorption  potential  is  still  an  open  problem.  Here  we  employ  a  semi- 
empirical  absorption  potential  as  discussed  by  Truhlar  and  coworkers  [117].  The  V^s  is 
a  function  of  molecular  charge  density,  incident  electron  energy  and  the  mean  excitation 
energy  ,  A,  of  the  target  : 

Kbs(r)  =  -p(r)(^0c/2)1/2(87r/5PfcJ)/f(fc2  -  k)  -  2A)(Aj  +  A2  +  A3),  (9.7) 


where: 


Vioc(r)  =  k2-  Vst(r)  -  Kx(r)  -  Vpol 

(9.8a) 

A\  =5ifcJ/2A, 

(9.8  b) 

A2  =  -kj(5k2  -  3 k})/(k2  -  kj)2, 

(9.8 c) 

,  ,  (2k2f  +  2A-  k2f!2 

^3  =  2H(2k]  +  2A  -  t2)  V-ty*  • 

(9.8<f) 

where  \k2  is  the  energy  of  the  incident  electron  in  Hartrees.  Here  H{x)  is  a  Heaviside 
function  defined  by,  H(x)  =  1,  for  x  >  0  and  H(x)  =  0  for  x  <  0.  Here  also,  we  consider 
only  the  spherical  term  in  the  single-center  expansion  of  the  absorption  potential. 

After  generating  the  full  optical  potential  (Eq.  9.1)  of  a  given  electron  molecule 
system,  we  treat  it  exactly  in  a  partial-wave  analysis  by  solving  the  following  set  of  first- 
order  coupled  differential  equations  for  the  real  (x<)  and  imaginary^/)  parts  of  the  complex 
phase-shift  function  under  the  variable-phase- approach  (VPA)  [118]: 

X't(kr)  =  -|[2VR(r)(A2  -  B2)  +  2V»ta(r)AB],  (9.9) 

X't{kr)  =  -|[2Fr (r)AB  -  2Vabs(r)(A2  -  B2)}.  (9.10) 
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where: 


A  =  coshxe(kr)[cosxt(kr)je(kr)  -  sin  Xt{kr)r)t(kr)],  (9.11a) 

B  =  -sinhx*(A:r)[sinx*(/i:r)j*(A:r)  -  cosxe(kr)rjt(kr)],  (9.116) 

and  je(kr)  and  T]t(kr )  are  the  usual  Riccati-Bessel  functions.  The  Eqs.  (9.9)  and  (9.10) 
are  integrated  up  to  a  sufficiently  laxge  r  different  for  different  £  and  k  values.  Thus  the 
final  5-matrix  is  written  as: 

Se(k)  =  exp(—2xt)  exp(i2xt),  (9.12) 

and  the  corresponding  DCS  axe  defined  as: 

tma*  2 

%  =  7H  E  (“  +  1)[«K*)  -  l]P<(cos  0)  .  (9.13) 

1=0 

where  P/(co s  0)  is  a  Legendre  polynomial  of  order  £.  The  integrated  elastic  (aej ),  absorption 
(crabs)  and  <rt  cross  sections  are  described  in  terms  of  the  5-matrix  as  follows: 


<4  =  +  1)|1  -  St(k) I2; 

lnux 

<?d  =  , 

(9.14) 

/==  0 

^b,  =  p(«  +  l)[l-lS,(»:)|2); 

^max 

^abs  =  Yj  °absi 

(9.15) 

1=  0 

<r‘  =  ^(2«+l)[l-ReS,(t)); 

^max 

=  Y 

1=0 

(916) 

We  note  that  of  =  +  aabs  is  the  contribution  from  the  spherical  term  only.  In  the  above 

analysis  the  inelasticity  or  the  absorption  factor  is  defined  by  |5*(Jfc)|  =  exp(—2xt)- 

We  presented  [10]  total  cross  sections  for  a  large  number  of  molecules  divided  in  terms 
of  isoelectronic  sequence  such  as,  2-  (H2),  6-  (Li2),  10-  (HF,  CH4),  14-  (CO,  N2,  C2H2, 
HCN),  16-  (02),  18-  (HC1,  H2S,  PH3,  SiH4)  and  22-  (C02)  electron  systems.  All  the 
calculated  results  *  ,re  compared  with  measured  data  and  we  found  a  very  good  agreement 
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between  theory  and  experiment.  In  particular,  the  present  theory  was  successful  above  100 
eV  for  almost  all  the  above  species. 

We  also  emplyed  our  high  energy  cross  sections  (ct^,  <?abs  and  <7t)  to  determine  Bom- 
Bethe  parameters  [109-113].  The  Bethe  theory  defines  the  total  inelastic  (<rabs)  cross 
section,  while  the  Bom  approximation  gives  the  total  elastic  (<rei)  cross  section.  Thus,  a 
combined  Bom-Bethe  theory  expresses  the  total  cross  section  in  terms  of  the  following 
analytic  formula: 


E  at 
R  nal 


+ 


4Mt2otln[4Ctot-l+... 


R 


(9.17) 


where  E  is  the  incident  energy  in  eV,  R  is  the  Rydberg  energy,  a0  is  the  Bohr  radius.  The 
constants  Ae\ ,  Bei ,  Cei ,  Mtot  and  Ctot  have  physical  significance.  The  Mt2ot  is  related  with 
the  optical  oscillator  strength  of  the  molecule.  To  our  knowledge,  there  were  hardly  any 
results  in  the  literature  on  the  Bom-Bethe  parameters  (Eq.  9.17)  for  most  of  the  molecules 
studied  here.  In  addition,  we  fit  a  general  formula  for  the  crt  from  our  calculated  values  in 
the  following  form: 

R  R  .  /? 

(9.18) 


E  at  ,  E  ifR  t 
—  — ^  =  a'ln  —  +b'—+c, 
R  47ra2  R  E 


and  provided  the  values  of  constants  a',  b'  and  c'  for  the  present  set  of  molecules  [10]. 

In  order  to  further  analyse  our  total  cross  section  results,  we  have  made  Bethe  plots 
(<7absI2/(47rIlao)  vs.  In (E/R))  for  all  the  molecules  [10].  From  these  plots  we  found  the 
Bethe  theory  valid,  which  says  that  the  Bethe  plot  should  be  a  straight  line  with  a  gradient 


of  AC- 


We  further  put  together  the  above  description  of  high  energy  electron  scattering  with 
a  variety  of  diatomic  and  polyatomic  molecules.  It  is  always  useful  to  know  the  variation 
of  cross  sections  in  different  targets.  We  have  already  seen  that  isoelectronic  collisional 
systems  possess  similar  cross  sections  both  in  quantity  and  quality.  Except  very  highly 
polarizable  Li2  or  the  simple  H2  molecules,  we  might  expect  some  kind  of  correlation 
between  at  and  any  molecular  property  (number  of  electrons  Z,  polarizability,  multipole 
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moments,  molecular  size  etc.).  A  general  feature  above  100  eV,  common  for  all  the  above 
targets,  is  the  decrease  of  <rt  with  increase  in  impact  energy.  We  clearly  found  [10]  that 
the  <rt  values  increase  with  an  increase  in  Z,  i.e.,  there  is  strong  correlation  between  the 
total  cross  section  and  the  size  of  the  target.  On  the  other  hand,  we  did  not  see  such  a 
correlation  with  respect  to  target  polarizability.  However,  at  lower  energies  (around  100 
eV),  the  at  seems  to  have  a  strong  dependence  on  ao-  From  our  <rt-Z  correlation  diagram 
[10]  a  rough  estimate  of  the  total  cross  section  can  be  made  for  any  molecular  system  with 
Z  smaller  or  greater  than  22. 
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11.  CONCLUSIONS 


During  a  period  of  three  years,  several  important  developments  and  calculations  were 
made  on  the  electron-molecule  scattering.  We  presented  rotationally  elastic,  inelastic  and 
summed  cross  sections  for  electron  scattering  with  several  polyatomic  molecules  (CH4, 
SiKU,  GeH4,  H2O,  H2S  and  NH3)  in  a  highly  sophisticated  close-coupling  non-empirical 
theory.  Exchange  effects  were  included  exactly,  while  polarization  corrections  were  consid¬ 
ered  approximately  but  without  involving  any  fitting  parameter.  Results  were  compared 
with  measurements  where  such  data  were  available.  No  previous  theoretical  investigation 
is  available  on  these  molecules  in  such  detail  and  energy  range.  The  present  study  has 
drawn  several  important  conclusions:  (1)  exchange  effects  at  the  exact  level  and  model 
polarization  potentials  can  describe  low  energy  e-molecule  scattering  quite  accurately,  (2) 
polarization  effects  should  be  included  in  a  theory  where  target  orbitals  are  relaxed  in  the 
presence  of  incoming  particle,  (3)  the  effect  of  gas  temperature  on  the  cross  sections  is  cru¬ 
cial  for  rotationally  inelastic  channels,  (4)  a  single-center  approach  can  be  useful  for  many 
polyatomics  without  any  convergence  problem,  (5)  the  iterative  scheme  is  quite  promising 
tool  in  future  if  employed  in  an  optimized  way  as  suggested  in  this  study. 

For  the  first  time,  a  large  variety  of  molecules  are  studied  at  intermediate  and  high 
energies.  The  total  cross  section  (sum  of  elastic  plus  all  possible  inelastic  channels)  for 
several  selected  gases  were  calculated  and  presented  in  tabular  form  in  the  energy  range  of 
10-5000  eV.  Such  theoretical  data  are  not  found  in  the  literature  for  most  of  the  targets 
studied  here.  In  addition  we  also  reported  Bom-Bethe  parameters  for  all  these  molecules. 
An  attempt  was  made  to  find  a  correlation  between  molecular  properties  and  the  corre¬ 
sponding  total  cross  section  quantity. 

Several  computers  codes  were  modified  and  developed  in  order  to  deal  with  big 
molecules  and  also  polar  ones.  Our  Single-center-expansion  scheme  was  made  possible 
to  work  for  much  heavier  systems  (e.g.  GeHU,  SF6,  CF4,  etc.)  than  the  typical  CH4,  etc. 
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ones.  An  interface  program  was  developed  which  generates  SCE  quantities  from  molecular 
wave  functions.  This  progress  was  realized  in  terms  of  some  actual  calculations  on  some 
molecules  of  larger  size. 

This  work  is  still  in  progress  in  two  ways:  (1),  to  investigate  electron  interactions  in 
large  systems,  and  (2)  to  employ  the  present  set  of  programs  for  vibrational  and  electronic 
excitation  of  polyatomic  molecules.  Some  preliminary  progress  has  been  made  for  the 
CH4  and  H2O  molecules.  However,  this  work  requires  a  huge  amount  of  effort  in  terms  of 
manual  and  computational  facilities. 
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1.  ABSTRACT 


This  is  to  be  considered  as  a  final  report  on  a  subsection  of  the  Air  Force  contract 
27256-PH-SAH  (Dr.  Ashok  Jain,  Principal  Investigator,  Dr.  Charles  A.  Weatherford, 
Co-Investigator) . 

The  present  portion  of  the  research  involves  an  attempt  to  go  beyond  the  standard 
approach  to  electron-polyatomic  scattering,  wherein  the  use  of  the  partial  differential  equa¬ 
tion  (PDE)  approach  to  electron  or  positron  molecule  scattering  in  the  low  to  intermediate 
energy  regime  is  utilized.  The  research  attempts  to  avoid  single-center  expansions  of  bound 
and  continuum  wavefunctions. 

The  body  of  the  report  consists  of  a  narrative  describing  the  rationale  for  three  papers 
which  constitute  the  appendix  (two  are  unpublished,  and  one  is  accepted  for  publication 
in  the  Physical  Review  A),  as  well  as  a  description  of  the  new  application  of  finite  elements 
to  the  PDE  approach  to  electron-polyatomic  scattering. 
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2.  INTRODUCTION 


The  present  document  c  utes  the  final  report  on  a  subsection  of  the  contract 
27256-PH-SAH  (Dr.  Ashok  Jain,  Principal  Investigator,  Dr.  Charles  A.  Weatherford,  Co- 
Investigator)  for  the  benefit  of  the  United  States  Air  Force.  The  present  project  is  entitled 
“Partial  Differential  Equation  Approach  to  Electron,  Positron-Molecule  Scattering.”  The 
research  in  scattering  theory  involved  the  calculation  of  electron  collision  cross  sections  for 
low  to  intermediate  energy  interactions  with  various  diatomic  and  polyatomic  molecules. 
This  subsection  involves  the  application  and  extension  of  the  Partial  Differential  Equation 
(PDE)  method1-10  to  electron-molecule  sc?  i~  In  doing  so,  single-center  expansions 
are  avoided  and  thus  new  methods  of  calculating  he  required  multicenter  integrals  have 
been  developed. 

Among  the  computational  techniques  that  have  been  used  to  compute  the  scattering 
cross  sections,  one  of  the  most  common  elements  is  the  expansion  of  the  continuum  wave 
function  in  partial  waves  -  that  is,  in  spherical  harmonics.  The  resulting  variational  treat¬ 
ment  leads  to  sets  of  coupled  integro-ordinary  differential  equations.11-17  Over  the  past 
number  of  years,  the  direct  solution  of  the  partial  differential  equations  (PDEs)  has,  been 
accomplished1-10  with  great  success.  The  code  is  still  in  the  development  stage,  how¬ 
ever.  The  code  is  being  generalized  to  use  a  finite  element  representation  of  the  Laplacian. 
Also,  no  attempt  has  as  yet  been  made  to  exploit  the  multiprocessor  environment  of  the 
newest  supercomputers.  In  the  present  research,  we  have  adapted  the  finite  element  PDE 
(FEPDE)  code  to  a  mutiprocessor  environment.18  Specifically,  we  are  using  the  massively 
parallel  environment  (Connection  Machine,  CM-5).  In  doing  so,  we  have  begun  to  use 
multigrid19-21  technology. 

In  the  remainder  of  this  report,  we  will  present:  section  3  briefly  explains  the  rationale 
for  the  three  papers  included  in  the  present  report  as  appendices;  section  4. A  describes  the 
basics  of  the  Partial  Differential  Equation  Theory;  section  4.B  describes  the  generalization 
to  a  finite  element  discretization;  section  5  contains  the  references;  the  Appendices  are 
given  as  section  6. 
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3.  RATIONALE  FOR  APPENDICES 


The  present  research  constitutes  an  attempt  to  establish  an  ab  initio  computational 
approach  to  electron  and  positron- polyatomic  scattering.  As  such,  it  avoids  single-center 
expansions  of  the  bound  and  continuum  wavefunctions.  The  longer-range  objective  is  to 
also  go  beyond  phenomenological  potentials  and  construct  a  method  with  an  ab  initio 
potential  including  polarization  effects. 

Appendix  6. A  (“Addition  Theorems  for  Coulomb  Sturmians  in  Coordinate  and  Mo¬ 
mentum  Space”  )  consists  of  an  as  yet,  unsubmitted  paper  describing  a  new  approach 
(developed  by  the  author)  to  calculate  all  required  bound-state  molecular  multicenter  in¬ 
tegrals  using  exponential  basis  functions.  The  method  is  still  being  developed,  but  show 
great  promise. 

Appendix  6.B  (“Schrodinger  Equation  Mesh  Requirements  in  the  Finite  Difference 
Discretization  for  Non-Spherical  Potentials”  )  describes  a  mesh  convergence  study  using 
the  finite  difference  discretization  of  the  continuum  Schrodinger  equation.  The  convergence 
turned  out  to  be  rather  disappointing.  This  along  with  ease  of  coding  considerations,  has 
led  the  author  to  adopt  the  finite  element  discretization  descrived  in  Section  V.B  of  this 
report. 

Finally,  Appendix  6.C  (“  Completion  of  hybrid  theory  calculation  of  the  IIS  resonance 
in  electron-N2  scattering  describes  the  latest  calculation  which  has  been  completed  using 
the  finite  difference  discretization  of  the  continuum  Schrodinger  Equation.  This  paper  is 
to  be  published  in  the  April  1994  issue  of  the  Physical  Review  A. 
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4.  PDE  APPROACH  TO  ELECTRON-MOLECULE  SCATTERING 


Narrative  A: 

Partial  Differential  Equation  Theory 

This  proposal  concerns  the  Quantum  Scattering  description  of  two  colliding  systems. 
It  considers  the  solution  of  the  time-independent  Schrodinger  equation  with  two  bodies  in 
both  the  initial  and  final  channels.  The  non-relativistic  velocity  regime  is  considered.  The 
collision  energy  is  low  enough  so  that  the  Bom  approximation  is  not  valid.  Although  the 
basic  equations  to  be  solved  axe  prototypes  for  a  large  number  of  collision  systems  in  atomic 
and  molecular  physics  and  quantum  chemistry,  the  theory  herein  presented,  specializes  to 
electron-molecule  scattering. 

The  central  idea  behind  the  theory  is  the  desire  to  avoid  partial-wave  expansions  of 
the  scattering  wavefunction.  This  is  particularly  desirable  when  the  interaction  potentials 
are  highly  non-spherical.  The  Partial  Differential  Equation  (PDE)  method  has  proven  to 
be  an  extremely  stable  and  efficient  solution  technique  for  the  close-coupling  equations 
that  occur  in  electron-molecule  scattering1-10.  This  has  proven  to  be  especially  true  on 
vector  computers.  It  is  anticipated  that  the  use  of  parallel  processors  will  also  allow  for  an 
even  greater  efficiency  of  computation. 

The  basic  equations  to  be  solved  are  illustrated  by  the  static  exchange  equation  3 

[V2  +  k2]F(f)  =2  V(r)F(f) 

Ne 

i=l 

This  is  the  most  simple  but  reasonably  realistic  equation  that  is  encountered  in  electron- 
molecule  scattering.  F  represents  the  Continuum  orbital;  r  represents  the  continuum 
electron  position  vector;  W,  the  exchange  kernel,  is  given  by 


$  is  a  bound  orbital  of  the  target  molecule.  Realize  that 


V2(|r  —  x|  — 47r6(f'— x).  (3) 

Then,  two  coupled  equations,  which  must  be  solved  simultaneously,  are  then  obtained 

[V2  +  k2]F(r)  =2  V(r)F(r) 

Ne 

-2£*„,(rW,,(r-) 

!=1 

V2^ai(r)  =  —  4n**a{(r)F(?)  (5) 

W  plays  the  part  of  a  pseudo-continuum  orbital.  It  has  the  value  of  zero  on  the  large 
r  boundary.  There  are  as  many  inhomogeneous  equations  like  (5)  as  there  are  bound 
molecular  orbitals. 

Eqs.  (4)  and  (5)  are  elliptic  PDEs.  They  are  solved  by  applying  a  finite  element 
approximation  to  V2  to  produce 


(4) 


AF  -  B  (6) 

A  is  the  block  tridiagonal  coefficient  matrix,  F_  is  the  solution  vector,  and  B_  is  the  boundary 
vector. 

At  the  r  =  p  boundary,  the  partial-wave  solutions  are  extracted  by  expanding  the 
PDE  solution,  for  each  symmetry  m,  in  partial  waves. 


*"•(?> = E  E  *)  (?) 

u  h 

The  a's  are  arbitrary  amplitudes.  The  f's  have  the  form 

fkij  (r)  =  sii(r)piiij  (r)  +  cii(r)Qiiij  ( r )  (8) 

The  S  and  C  are  spherical  Bessel  functions  if  the  target  is  neutral  and  Coulomb  functions 
if  the  target  is  charged.  The  value  r  =  p  is  determined  by  the  range  of  the  exchange  con¬ 
tribution.  The  P  and  Q  can  be  combined  into  a  multidimensional  vector  I  and  propagated 
to  an  approximate  oo. 
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I(r/)  =  T(r/)I(rt) 


(9) 


T  is  the  propagator.  The  reaction  matrix  is  then  computed  by 

K  =  QP -1  (10) 

We  have  generalized  the  PDE  method  to  treat  charged  targets  and  multichannels  in 
the  fixed  nuclei  approximation  (FNA).  This  means  that  the  S  and  C  are  Coulomb  funtions 
for  charged  targets.  To  treat  general  open  shell  targets,  we  must  expand 

TIT 

=A  ^  *i(l,  •  •  • ,  N)F-j(N  +  1) 

n!=1  (li) 

+  f>(l . N  +  l  )C.J 

a=l 

Here,  the  /3’s  are  correlation  terms.  This  results  in  the  equation 


[V2  +  *?]*!?(?)  =2f;[V.,(r)  + 

i= i 

+2^!7i„(f)Ca4 

a=l 

Note  that  the  Lagrange  multipliers  are  not  needed  since  the  orthogonality  is  enforced  by 
the  specification  of  the  value  of  0  on  the  r  =  p  boundary  in  the  B  vector.  They  can  be 
put  in  as  an  internal  check  (which  we  do).  The  C  in  (12)  satisfies 


£nalcM  +  £<tfi,,|F£)  =  o 


(13) 


6=1  t'=l 

where  a  =  (1, . . . ,  nc)  and  j  =  (1, . . . ,  nt),  and  j  is  fixed.  U  and  Q,  axe  defined  by 


ne 


a,=l 


Qab  =  (/?a|-HiV+l  ~  E\Pb) 


(14) 


(15) 
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We  might  note  that  the  correlation  terms  are  of  two  types:  (1)  those  that  axe  required  to 
relax  orthogonality  and  ;  (2)  those  that  are  required  to  produce  projectile-target  correlation 
(both  short  and  long  range). 

Eqs.  (12)  are  solved  in  the  standard  manner.  See  for  example,  Ref.  22.  In  this 
procedure,  F  is  expanded  as  (drop  the  m  symmetry  notation) 

=  4°V">  +  E  *ia\r)C.k  (16) 

a=l 

Then  the  actual  equations  to  be  solved  are 


[V2  +  i?]40,(f)  =2£[Vy(r)  +  (17) 

i= i 

[V2  +  k^\r)  =2  £>y(r)  +  W*°(r>]ff ’(?) 

i= l  (18) 

+2  Uia(f) 

The  W’s  are  here  defined  by 


V2wf  >(f)  =  — 4ir  £  5*,(r)  Fg'ff)  (19) 

k= 1 

V2I^;>(r)  =  -4t  £  «**(r|/f  >(r)  (20) 

fc=  1 

Thus  the  formalism  is  similar  to  that  of  Ref.  22.  The  primary  difference  is  that  we  are 
solving  PDE’s  instead  of  ordinary  differential  equations. 

Each  of  the  Eqs.  (17)  and  (18)  may  be  solved  independly  on  different  processors  and 
then  combined  after  solution.  For  those  correlation  terms  included  in  Eq.  (12)  that  are 
not  required  to  relax  orthogonality,  an  optical  potential  will  be  constructed  in  the  manner 
of  Schneider  and  Collins23. 
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Narrative  B: 


Finite  Element  Discretization 


A  short  outline  of  the  Finite  Element  method  is  presented  below: 

•  Time-independent  Schrodinger  equation  for  e- -molecule  scattering  in  the  Fixed  Nuclei 
Approximation  in  3D  using  spherical  coordinates: 


d 2  2d  Id2  cotO  d  1  d2 

. dr 2  r  dr  r2  d8 2  r2  d0  r2sin2B  d<f> 2 


+  fc2]*(r,M)  = 


2  V(r,M)*M,f> 


•  Solve  by  the  method  of  Finite  Elements: 


(21) 


1.  Divide  solution  space  into  non-overlapping  domains  called  “finite  elements”  ; 

2.  Define  local  basis  functions  within  each  element-in  the  present  case,  5th  order  Hermite 
polynomials; 

3.  Discretize  in  each  element  by  the  Galerkin  method  and  implicitly  form  global  element 
matrix; 

4.  Form  global  connected  matrix  by  enforcing  continuity  of  the  solution  and  its  deriva¬ 
tives  across  element  boundaries; 

5.  Form  boundary  reduced  connected  matrix  and  boundary  vector;  A;  B  3 


AX  =  B 

6.  Solve  simultaneous  equations  by  direct  LU 

X  =  A~1B 

•  Hermite  Interpolation 


-1 

0 

+1 

r(<) 

ro 

r(e)  +  !/,(*) 
r0  *  2  nr 

r$e)  +  ti 

(22) 


(23) 
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In  terms  of  local  variables  x,  :  {— 1,0,  +  l}i=1 


v>(e)(z)  =  XI  ^.-e)P‘(z) 

»=a 

where 


rtf'! 

'  rj>{e)(xi)  ' 

£V>(e)L 

^(e)(x2) 

i^(e)u. 

^(e)(x3) 

UH 

L£v>(e)lJ 

and  pi(x )  is  determined  from 

6 

Pi(x)  =  ^9ij 
j= i 


(24) 


(25) 


(26) 


Pi(xj)  =  6ij  for  i  =  1,3,5  and  je{x i,x2,x3} 


(27) 


^P*(x)Ui=^«i  /or  *  =  2,4,6  and  je{xi,x2,x3} 


(28) 


Pi(x)  =  9n  +  012Z  +  513X2  +  9ux3  +  gi5x 4  +  ^16x5  (29) 

^p2(x)  =  <722  2jf23Z  +  3^24®2  +  4025*3  +  &526X4  (30) 

p3(x)  =  <731  +  (J32X  +  gzzx2  +  g34x3  +  035x4  +  53ex5  (31) 

"j“P4(x)  =  <742  +  2<743X  +  3^44  x2  +  4g45X3  +  5^46  X4  (32) 

ax 

P5(x)  =  051  +  052*  +  9SZX2  +  g54X3  +  055X4  +  056X5  (33) 

^Pe(x)  =  562  +  2^63 X  4-  3<764*2  +  4565X3  +  5^66 x4  (34) 

Now  implement  conditions  (27)  and  (28) 
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-061 
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.-1 
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0 

1 

5. 

(35) 


The  expression  for  in  terms  of  global  variables  can  now  be  written  using 


r  =  r(0  +  k(«)(l+fi 

In  order  to  go  to  a  global  representation  of  (4),  tp(e\r),  it  is  not  necessary  to  let  x  — +  r 
in  Pi(x).  The  globed  operators  can  be  transformed  to  local  operators  which  operate  only 
on  the  piecewise  continuous  polynomials  within  each  element 


d  dx  d  2d 

dr  dr  dx  dx 


(37) 


So  we  can  simply  say 


f  2  ]2  cP 
dr2  hff)  dx2 


(38) 


6 

Vic)(r)  =  ^iJ;\e)pi(x)  (39) 

i=l 

However,  the  problem  with  (19)  is  that  rp\e\  for  i  =  2,4,6,  represents  derivatives  of 
\lde)(r)  with  respect  to  the  local  variable  x  and  we  need  it  to  represent  derivatives  with 
respect  to  the  global  variable  r. 


thus 


dxp('t\x)  _  dr  dip^^r)  __  1  ^(e)  dtp^\r) 
dx  dx  dr  2  r  dr 


^(£)(r)=  ^ie)pi(X)+  J2  Pi(X) 


<=1,3,5 


<=2,4,6 


(40) 


(41) 
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Let  us  now  write 


6 

tf(e)(r)  =  ^^e)P<e)(x)  (42) 

«=i 

such  that 


and 


for  i  =  1,3,5  1 
for  i  =  2,4,6  J 


(43) 


Pi(x) 
h[e)pi(x ) 


for  i  =  1,3,51 
for  i  =  2,4,6  J 


(44) 


In  Fig.  1,  a  short  description  is  presented  of  the  derivation  of  the  scattering  amplitude 
in  the  fixed-nuclei  approximation  for  electron-polyatomic  scattering. 

In  Fig.  2,  a  comparison  at  the  static  level  is  made  between  the  results  of  Jain,  using 
the  single-center  partial- wave  method,  and  the  present  finite  element  method  for  e-CH4,  in 
the  Ai  symmetry.  This  demonstrates  that  the  basic  discretization  is  working.  It  remains 
to  incorporate  the  exchange  effect  and  to  account  for  polarization. 
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FIXED-NUCLEI  AMPLITUDE  FOR  e-POLYATOMIC  SCATTERING 


Incoming  Wave  Transformed  To  Body  Frame 

r,“(n-)  =  ^r,-(n)D  <';r(A,) 

m 

Scattering  Of  Y,m(S7)  By  Polyatomic  In  Body  Frame 

srrtfi)  -  E  E 

h  m3 

Transform  Back  To  Lab  Frame 

m' 

In  Total 

General  Form  Of  Fixed  Nuclei  Scattering  Amplitude: 

/(&;«')  =  E 

Figure  1.  Short  derivation  of  the  form  of  the  scattering  amplitude  for  e-polyatomic  scattering 
in  the  fixed-nuclei  approximation. 
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ADDITION  THEOREMS  FOR  COULOMB  STURMIANS 
IN  COORDINATE  AND  MOMENTUM  SPACE 

by 

Charles  A.  Weatherford 
Physics  Department 
and 

Institute  For  Molecular  Computations 
Florida  A&M  University 
Tallahassee,  FL  32307,  USA 

1.  PRELIMINARIES 


The  problem  of  performing  multicenter  integrals  over  Slater-Type  Orbitals  is  one 
of  the  classic  problems  of  Quantum  Chemistry.1  Two  of  the  most  successful  techniques 
for  the  evaluation  of  such  integrals  are  the  Lowdin  a-function  method2-5  and  the  Fourier 
Transform  technique.6-8  An  interesting  variant  of  the  Fourier  Transform  technique  is  based 
on  the  Fock  hyperspherical  projection  of  the  momentum  space  hydrogenic  orbitals  onto  a 
four  dimensional  hypersphere.9’10  Some  important  applications  of  the  Fock  projection  to 
the  Slater  integral  problem  are  given  in  references  11  and  12. 

We  are  interested  in  performing  the  following  4  —  center  electron  repulsion  integral: 


*  S’  C,  D)  = 

x  J  *-Sprr^j +  c  +  s ) 

The  5— functions  are  Coulomb  Sturmians  (CSs).13  The  CSs  are  defined  by7 


(1) 
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S”,(o,f)  =  (-l)”+'-‘a 3>2N„I  e-"LZ*i\(2ar)yr (2a? )  (2) 

where  the  phase  factor  (— l)n+,_1  has  been  introduced  in  the  manner  of  Koga  and  Matsuhashi.14 
Then 


Nnt  =  23/2 


(n  —  /  —  1)! 


1/2 


(3) 


2n(n  +  Z)! 

is  the  normalization  constant,  L  is  a  Laguerre  polynomial  (same  definition  as  given  in 
Arfken15)  given  by 


£iw(*)«  ("^)  »«(-»! 0  +  i;»)  (4) 

where  i.Fi  is  a  confluent  hypergeometric  function16  and  y  is  a  solid  spherical  harmonic 
given  by 


JT(f)  =  >-,rl”(f)  (5) 

and  y  is  a  regular  spherical  harmonic.15  The  series  for  L  terminates  in  the  present  case 
and  gives 


n— l 


£i2iTi!(2ar)  =  2Br„,(2«r)r-1 


r=l 


(6) 


HQ^Hn  +  Q! 

rn/  (n  -  l  -  r)!(2/  +  r)!(r  -  1)! 

The  J5  integral  is  the  most  general  integral  which  is  encountered  when  non- correlated 
basis  functions  are  used  in  an  LCAO  procedure  for  Cl  calculations.  In  the  present  work, 
addition  theorems  (which  separate  integration  variables)  are  derived  for  Coulomb  Sturmian 
functions  and  modified  Coulomb  Sturmian  functions  (’S'), 


*3(«,f)s£s3(«,f)  (8) 

ar 

in  coordinate  and  momentum  space.  The  Coulomb  Sturmians  and  their  Fourier  Trans¬ 
forms,  are  complete  in  the  sense  of  distributions  in  Sobolev  space7  W^^R3),  which  is 
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a  proper  subset  of  Hilbert  space  L2(R3).  In  this  regard,  a  standard  Slater  Type  Orbital 
(STO)  defined  by 


XnlM  =  Nnan+1'2rn-1e-arYlm>(r ) 


Nn  = 


2n+a 

[(2n)!]i 


may  be  expanded  in  a  finite  linear  combination  of  Coulomb  Sturmians  and  vice  versa.17 


n/=l<  +  l 


L(nO.  _  /I  \n,+li  \2n  fih  +  *)("/  ~  U  ~  1)KW/  +  lj) [I  * 
n'  K  (2n*)! 

(~l)g(n,  +  lj  +  <?  -  1)! 

i rl  (nf  -  li  -  cr)\{2h  +  -  1)! 


and  vice  versa 


ST.  •„,(«, r“)  =  a-'/2  £  (13) 

n/=f,+l 

a(nt)<  =  (n»  ~  jj_ ~  l)-(n»  +  ^)-  2 

,  n<  (14) 

x _ l(n/)I)* _ 

(nj  -  n/)!(n/  +  li)!(nf  -  li  -  1)! 

It  is  useful  to  “equalize”  the  screening  parameters  in  one  electron,  single  center  densi¬ 
ties,  composed  of  CSs  basis  functions.  To  accomplish  this,  it  is  first  useful  to  equalize  the 
screening  parameters  for  densities  composed  of  STOs.  Some  trivial  algebra  produces 
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where 


0)  — 


Nni(a)Nna(0) 

NnMNnM 


where  7  =  and  T  is  the  STO  one  electron  density  defined  by 


=  X{Zl(0r  ) 

Then,  defining  the  one  electron,  single  center,  CS  density  by 


^o7(W“’M  s  sfcS&nsTM*’?)  (18) 

the  required  expression  (which  equalizes  the  CS  screening  parameters  for  the  density  Cl), 
is  obtained  by  first  using  Eq.  (13)  for  each  CS  in  Cl,  then  using  Eq.  (15)  to  equalize 
the  screening  parameters  in  T,  and  then  using  Eq.  (11)  to  recover  the  required  Cl  with 
equalized  screening  parameters.  The  result  is 

ni  713  712  71 4 

E  E  E  E 

ns=Jt+l  hj=/i+1  n*=/j+l  ng=lj+ 1 

(19) 

where 

=  2§7^7 (20) 
It  is  also  useful  to  linearize  the  product  of  two  CSs.  To  this  end,  we  linearize  the 
product  of  two  STOs  as  per  the  following  formula: 

li+h 

E  D!'m)'('m,’C+nTiu,(“  +  ^.0  (21) 

/3=|fi-/j|,2 

where 


JV„,(q)iV r„,Q?) 

'■'n\nt\0ii  P)  jit  (  ,  q\ 

Nni+n2-i{a  +  p) 
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and 


p(/m),(/m)2  _  FCg/l  +  l)(2k  +  l)(2/3  +  1)  1/2 

fh  h  I3  \  f  h  h  h  \ 

y  0  0  0  )  \7Tii  —m2  m2  —  mj  j 

Thus,  the  linearization  of  the  Cs  density  is  easily  shown  to  give 


nt+»»2— 1  li+h 

E  E  +  (24) 

n/=|/l— /j|+l  //  =Ml  — ^2l»2  _ 


where 


nf  n(ltnh(lmh 
[a(3(a  +  /3)]1/2 


D(lmh(lmh  £  fl(n0l  £  ajJ)i 


n3=l2+l  n<=ii+l 


It  is  important  to  note  that  all  of  the  sums  in  Eqs.  24  and  25  are  finite. 

The  CSs  satisfy  a  partial  differential  equation7  given  by 

[Vj -  a2]S3(a,r)  =  - 

This  can  be  contrasted  with  the  partial  differential  equation  satisfied  by  the  STOs 


[Vl-  a!kS(a,r)  =  Vnl(a,r)xZ(a,r) 


where 


Vni(oc,r)  =  - 


2am  [/(/  +  1)  —  n(n  —  1)] 


It  can  be  easily  demonstrated  that  the  CSs  satisfy  a  homogeneous  Fredholm  integral  equa¬ 
tion  of  the  second  kind15 
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(29) 


ST,(a,r)=~^  J  <SB(o,|r-2|)*3(a,i) 

where  G  is  the  Green’s  function  (same  as  for  the  modified  Helmholtz  equation)  and  is 
given  by 


<30> 

The  Green’s  function  satisfies  the  partial  differential  equation 

[V2?-  c2]G(a,  |f-  f|)  =  -6 (f-  x)  (31) 

The  translation  operator  Tj  =  (note  that  V  r  has  to  be  expressed  in  Cartesian 

coordinates  for  this  expression  to  be  valid)  is  defined  by 

TAF{?)  =  F(r  +  A)  =  F(rA)  (32) 

where  rA  =  r  +  A.  Now  an  important  characteristic  of  the  Laplacian,  which  is  well  known, 
is  its  invariance  with  respect  to  translations 

K+a  =  (33) 

Applying  the  translation  operator  to  both  sides  of  Eq.  (26)  and  using  Eq.  (33)  yields 

[Vi  -  o2]S3(a, f  +  A)  =  +  A)  (34) 

It  may  be  noted,  in  passing,  that  Eq.  (34)  can  be  used  to  generate  an  algebraic  eigenvalue 

— * 

equation  by  expanding  5^)(a,  r  +  A)  in  terms  of  5jjj(a,f)  and  letting  operate.  Then 
premultiplication  by  5*™(a,  r)  is  used  and  finally,  integration  over  dr.  The  result  is  an 
eigenvalue  equation  in  which  the  exact  eigenvalues  are  known.  However,  the  coefficients 
(eigenvectors)  in  the  aforementioned  expansion,  which  result  from  the  diagonalization  of 
the  coefficient  matrix,  sue  the  desired  quantities  and  unfortunately,  these  components  have 
first  order  errors  in  them,  as  opposed  to  second  order  errors  in  the  approximate  eigenvalues. 


Ill 


Thus  this  procedure  proves  too  slowly  convergent18  to  be  useful,  without  some  systematic 
way  of  correcting  the  eigenvectors. 

The  integral  equation  solution  to  Eq.  (34)  can,  nevertheless,  be  written 

Snl(a>r  +  A)  =  ^  J  dx  G(a,|r-xl)$;3(a,£+A)  (35) 

This  can  be  contrasted  with  the  result  of  applying  the  translation  operator  directly  to  both 
sides  of  Eq.  (29)  to  give 

■$£}(<*»  r  +  A)  =  ~  J  dx  G(a,\r  + A- x\)^(a,x)  (36) 

Finally,  the  CSs  satisfy  the  orthonormality  conditions7 

/*  ^(nOi  (**’ r  )  =  6(nlm)1(nlm)2  (37) 

and 

/  drS(nl) *  hnlmMnlmh  (38) 

and  the  closure  conditions19 

i  =  E  I  (39) 

nlm 

and 

i-El*3)<S3l  (40) 

nlm 

and  the  single-center  (equal  screening  constant),  coordinate-space  overlap  of  two  CSs  is 
given  by7,18 


r-jmi  m2  _ 

U(nlh(nl)2 


—  [Tjh/j  ^»i2,ni+l  +  ^n2,nx  +  7n2l2^n2,ni—  l] 


where 


(41) 
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(42) 


±  _  r(n  =F  0(n  ±  l  ±  l)ii/2 

lnl  L  4n(n±l)  J 
When  a  sum,  such  as  Ylnim  indicated,  it  means 


n  : 
/  : 


1  — ►  oo 

0  — ►  Tl  —  1 


^  m  :  —/,—/  + 1, ..., +/ 

or,  the  order  could  be  changed  so  that 


(  l  :  0  — *  oo 

m:  —/,—/  + 1, ..., +/ 
n  :  /  +  1  — >  oo 


The  particular  order  used  needs  to  be  numerically  tested.  Note  that  the  integrands  of 
Eqs.  (37),  (38)  and  (41)  are  composed  of  CSs  with  the  same  screening  parameters.  The 
expression  given  in  Eq.  (41)  is  derived  from  Eq.  (37)  using  the  following  recursion  relation: 


«rS3(<*.r)  =  S?+I/a,r) 

xnSS(a,r-)  +  S”,j(a,r) 

These  equations  will  prove  useful  in  the  discussions  given  below. 


(43) 


2.  FOCK  HYPERSPHERE  AND  ADDITION  THEOREMS 


The  derivation  of  computable  expressions  for  the  integral  in  Eq.  (1)  involves  the  use 
of  the  relationship  between  Coulomb  Sturmians  and  four  dimensional  spherical  harmonics 
(Fock  hyperspherical  projection)  and  proceeds  by  representing  the  translation  operator  in 
coordinate  space,  and  the  translation  operator  in  momentum  space,  on  a  basis  of  Coulomb 
Sturmians  and  their  Fourier  Transforms  respectively.  The  resulting  matrix  representations 
of  the  translation  operator  in  coordinate  and  momentum  space  are  respectively,  finite  linear 
combinations  of  Coulomb  Sturmians  and  their  Fourier  Transforms,  and  are  both  unitary. 
This  can  be  demonstrated  by  the  use  of  closure  relations  for  the  Coulomb  Sturmians  and 
their  Fourier  Transforms,  which  will  be  derived. 
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Filter  and  Steinborn20  have  given  an  addition  theorem  for  orthonormalized  Slaters 
(ONSs,  orthonormalized  functions  built  up  from  the  STOs  by  analytical  Gram-Schmidt 
orthonormalization)  by  using  the  translation  operator  explicitly  and  deriving  a  formula  for 
the  matrix  elements  of  the  translation  operator.  These  results  have  been  applied  to  the 
calculation  of  the  ground  state  wavefunction  of  H^“  with  somewhat  disappointing  results.21 

The  present  work  derives  an  addition  theorem  for  the  CSs,  similar  in  structure,  to  the 
addition  theorem  of  Filter  and  Steinborn  for  the  ONSs.  The  present  derivation,  however, 
is  quite  different  in  approach.  As  is  presented  below,  the  present  work  exploits  the  Fourier 
transform  of  a  CS7  and  the  relationship  of  the  Fourier  transform  to  four-dimensional 
spherical  harmonics22  using  the  Fock  hyperspherical  projection.9  The  method  is  similar 
in  many  respects  to  techiniques  of  Shibuya  and  Wulfman10  and  to  the  work  of  Alper11  and 
Novosadov.12  The  differences  between  these  techniques  and  the  present  technique  will  be 
implicit  in  the  derivation  presented,  and  in  the  approach  by  which  the  addition  theorem 
is  used  to  evaluate  the  multicenter  :ntegrals.  The  work  of  Alper11  is  perhaps  most  closely 
related  to  the  present  methods. 

The  symmetric  definition  of  the  three-dimensional  Fourier  transform  is  first  applied 
to  a  CS  so  that 


S3(«.f)  =  (  2*)~3/2j 

1  dp  etT' S„?(a,p) 

(44) 

S„?(a,p)=(  2tt)-3/2 

l  df  e-'?'  SZM 

(45) 

where  Snl  is  the  Fourier  transform  of  the  CS  (CSFT).  Thus,  it  can  be  shown7  that  the 
single-center  (equal  screening  constant),  momentum-space  overlap  of  two  CSs  is 


°r.i "m, = / (46) 

=  +  l  +  ^n2,nl  +  7n2f2 ^n2,m -1  ] 

The  Fourier  transform  of  can  also  be  taken.  Weniger7  demonstrated  the  relationship 
between  a  CSFT  and  a  four-dimensional  spherical  harmonic  (the  FDSH  is  the  one  used  by 
Judd,  Ref.  23)  given  by 
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y»?(r,M)  = 


(*2+p2)2o 


Snl(^P) 


(47) 


4a5/2  n' 

where  the  relationship  between  r  (hyperangle)  and  a  (screening  parameter)  is  defined  by 


e  = 

V  = 

C  = 

X  = 


2apx 
a2  4-  p2 
2<*Py 
a2  +  p2 
2  opz 
a2  4-  p2 


a2  4-  p2 


Thus  a  is  a  scaling  parameter  and 


sin(T)sin{9)cos(4> ) 
sin(T)sin{9)sin(4> ) 
sin(T)cos(0) 
cos(t ) 


(48a) 

(486) 

(48c) 

(48d) 


£2dV  +  C2+X2  =  l  (49) 

It  is  clear  then  that  a  point  p  of  three-dimensional  space  is  mapped  (one-to-one)  onto  a 
surface  of  a  four-dimensional  unit  sphere  described  by  the  angular  variables  f ,  rj,  (,  x-  The 
following  relationships  are  then  determined. 

sini(r)=  (50) 

dCla(p)  =  sin2(T)sin(0)dTd6d<j>  =  [  — — ■  2  ]  dp  (51) 

Do  not  confuse  the  four-dimensional  spherical  differential  surface  element  dCl  with  the 
single  center  charge  density  of  Eq.  (18).  The  orthonormality  conditions  for  the  CSFTs 
and  FDSHs  are  given  by 

r  dr  sin2(r)  f  d0  «n(«)  [  d*Y?J£(T,e,+)im 

Jo  Jo  Jo 

=  J  dP  S(n ih\ai,pj*{Z]2(a,p  )  (52) 

=  J  dP^(n/hl(a’P)^(nl)2(a’P) 

=  ^(n/m)i(n/m)  j 
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where 


TrTm/  -  \  (a2  +  J32)tt  to,  _ 

*nl(<*7P  )  =  - Snl  (<*’P  ) 


2a1/2 


(53) 


(a2  +  p2) 

The  Shibuya  and  Wulfman10  expansion  of  a  plane  wave  in  terms  of  CSs  and  FDSHs  is 
given  by 


=  [a2V]‘/J  £  _i__S„"Xa,  r)Y^t  m[(2Q(p)]  (54) 

nlm  '  ^  ' 

Eq.  (54)  may  be  compared  with  the  expansion  of  a  plane  wave  in  terms  of  CSs  and  CSFTs 
given  by 


or 


e*'  =  (2*)3/2  £  *„■,  ”W  )«•  (a,  f  ) 

nlm 


(55) 


e.r-r  =  (2t)3/2  £  S„',  ”(<.,?)*J(o,r-)  (56) 

nlm. 

For  later  use,  the  Fourier  transform  representation  of  a  three-dimensional  Dirac  delta 
function  and  the  Coulomb  potential  may  be  written  as 


*<*-*>= (2 


r(j>2-px) 


(57) 


[f*i  —  r2|  »I* 


0«p(ri-rj) 


|rx  -r2|  2tt2  J  p2  ^ 

The  derivation  of  the  addition  theorem  begins  by  applying  the  translation  operator 
defined  by  Eq.  (32)  to  both  sides  of  Eq.  (44).  The  result  is 


S3(«,?±  R  )  =  (2 tc)-3'2  I  dpe e**  '  S. ?(<*,?)  (59) 

The  Shibuya  and  Wulfman  plane  wave  expansion  of  Eq.  (55),  is  used  for  the  r-containing 
plane  wave  term  in  the  integrand  of  Eq.  (59).  The  result  is 
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(60) 


R  ^  ~  Q(nlTn)i(nlm)f{Qi  ±R)S(nl)  f(a’ 

(ntm)j 

where 

<?(.l«),(^-)/  (a,  ±S)  =  I  <&***  (61) 

Then  Eqs.  (51)  and  (53)  along  with 

sr^  =  2^2  [2K,0.[^(P)1  +  l?o[0»(p)l]  (62) 

is  used  to  transform  the  volume  integral  over  p  to  an  integral  over  the  surface  of  a  four¬ 
dimensional  unit  sphere  with  the  result 

Ql.taM.M,  (“.**)  =  /  <6.(fl«±iM  (63) 

Note  that  F^Q)  =  {tv21I2\~1I2  and  F2°0(ft)  =  [tt2  1/2]-1/2cos(r).  The  exphcit  evaluation 
of  the  integral  in  Eq.  (63)  proceeds  by  using  Eq.  (53),  (55),  and  (62)  along  with 

Y«mW)  =  (-l),-mFn?(n)  (64) 

and  a  very  important  linearization  property  of  FDSHs.23  This  is  a  property  of  the  FDSHs 
which  allows  the  product  of  two  FDSHs  on  the  same  hypersurface,  to  be  written  as  a 
finite  linear  combination  of  a  single  FDSH  on  the  same  surface,  and  is  central  to  the 
present  strategy.  This  expression  was  earlier  by  Shibuya  and  Wulfman10  and  used  to 
great  advantage  by  Alper.11  Actually,  the  result  is  probably  most  directly  attributable  to 
Biedenharn22.  It  may  be  expressed  as 

f»i+n2— 1  I1+I2 

x  5Z  C(n1l1m1,n2l2m2-,n3l3,m1  +  m2)  F™,l3+rri2[ftQ(p)]  (65) 

ns=|ni-n2|+l,2  (3=  j/i  — *2 1 ,2 

The  ‘,2’  on  the  sum  indicates  every  other  term  is  used.  The  C  is  defined  by 
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C(nilimi,n2l2m2’,n3l3m3)  = 


(2/i  +  l)(2/2  +  I)njn2n3 
2^ 


1/2 


x  (/im1,/2m2|/3."*i  +  m2) 


The  (•  •  •)  term  is  a  Clebsch  —  Gordon  coef  ficient,  and 


("1-1) 

("1—  1) 

Ll 

2 

2 

2 

("2-1) 

("2-1) 

L 

2 

2 

2 

("3  —  1) 

(n3-l) 

ii 

(66) 


I": 

is  a  Wigner  9 j  coefficient.24  The  efficient  evaluation  of  the  9 j  symbol  is  critical  to  the 
ultimate  computational  success  of  the  present  program.  However,  present  and  projected 
computer  technology  offers  great  encouragement  in  this  regard.  Even  now,  textbooks  have 
reasonably  adequate  computer  codes  included  in  the  appendices.25 
The  direct  evaluation  of  Q  is  now  straightforward,  resulting  in 


"i+"2  mtn(n— 1,/i+lj)  f'fiT'l 

x  a-3/2  V"'  ^(»Jm)i(«Jm)j  ^*771  l — 7712  (  t?'  ' 

/  v  /  gj  Til  Til 

n=mai(|tii— nj|,l)  lj|,2 


Then 


^(n/m) !(n/m) 2  _  2(_i)/j  m2 ir5/2B[(nl,mi  —  m2),(nlm)i,(nlm)2]  (68) 

and 


B[(nl,mi  —  m2),(nlm)i,(nlm)2]  =  2T[(nl,mi  —  m2),  (nlm)i,  (n2h,  — m2),  (100)] 

r[(ni,m!  -TO2),(nIm)1,(rtji2,-m2),(200)] 


(69) 


where  in  general 


r[(nZm)i(nJm)2(nlm)3(nfm)4]  2 

J  aa„(p) 


(70) 
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The  evaluation  of  T  proceeds  by  using  Eqs.  (52)  and  (65)  to  give 


r[(n/m)1(n/m)2(n/m)3(n/m)4]  = 
x  £  C(n3l3m3,nAUm4,n5l5,m3  +  m4) 

n8=|n3— n^l+1,2  ls  =  \l3  — 14\,2 
nS+>»2  — 1  h+l2 

x  E  E  C{ri2l2m2,nsh,rnz  +  m4,n5l5,m2  +  m3  +  m4) 

ns  =  |ns  — nj|+i,2  |/s  —  ijl.2 

Now,  an  addition  theorem  for  ’i>  can  be  derived,  in  an  analogous  fashion  to  the  one 
given  in  Eq.  (60)  for  S  with  the  result  of 

4^)l(tt,f±H)=  £  <3(„, (a, ±£)*"'); (a, r)  (72) 

(nlm)/ 

Note  that  the  coefficient  matrix  Q  is  the  same  as  in  Eq.  (60). 

In  a  very  similar  manner,  an  addition  theorem  can  be  given  for  S.  It  is 

5T4(a,p±p)=  £  ±^4,03,;!)  (73) 

with 

(.(.),  (a,  0,±P)  =  j  *(^7  (ft  ?  )S”j,  (a,  f )  (74) 

The  explicit  evaluation  of  f?  is  then  possible  and  results  in 


■R(n/m)j(nIm)i(Q(>  ^5  — 

nj+nj  min(n— 

x  E  E  G5'”),("'")i(o,/9)Cr'"’(«+«^ 

n=mai(|n1— n2|,l)  /=|/i  —/a  1,2 

(75) 

and 


G^lm)l{nlmh(a,0)  =  (27r)3/2n-2i-a--^)P^mh(nfm)2(a,;g) 


(76) 
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and  where  P  is  given  by  Eq.  (25). 

Finally,  an  addition  theorem  can  be  given  for  'if.  The  result  is 

^J,i)i(Q,P±p)=  ^2  nO/C&P)  (77) 

(nlm), 

where  it  can  be  easily  shown  that 

7(nlm)j(nim);(ai^iip)  ==  Q(nlm)  Tp)  (7®) 

The  R  and  Q  matrices  can  be  shown  to  be  unitary.  This  can  be  accomplished  using 
closure  relations.  The  closure  relations  for  the  FDSHs  were  given  by  Schwinger26 

S(  n  -  n ')  =  E  Y«(a  W  ( «)  (79) 

nlm 

In  this  regard,  an  important  relationship  is 

2  2 

« [n„(p)  -  n„ '(?')]  =  (2L2S£")s<(p-  p ')  (so) 

The  following  completeness  relations  can  then  be  derived  from  Eqs.  (39),  (40),  (79),  and 

(80) 

6(n  -  ft)  =  £  »3‘(a.irj)S3(«,fi)  =  £  S^(a,r2)4S(“.n)  (81) 

nlm  nlm 

f(fi  -*) = E^no.ftKto.p-i) = E^"(“’ft)*^(a’P>)  (S2> 

n/m  nlm 

The  unitarity  of  Q  and  i?  can  then  be  readily  demonstrated.  Therefore 

HI  ^(a’^)](„Jm)i(n/m)3^  (a’^l(ntm)3(nJm)2  =  ^(n,mh(nim)2  (83) 

(nlm)  3 

X/  (a>-S)](nlm)1(nlm)3^a,^(„lin)»(nlm)2  =  ^(nl”»)i(n|m)2  (8^) 

(nlm)s 

and  similarly  for  i? 
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^  y  [fi(®»  '^J)](n/m)i(n/m)3[—  (*■*’  @  i  ■^,)l(nim)3(nlm)j  ^(nlm)i(nlm)2  (88) 

(nim)s 

y!  LS  (a,^,i2)](nim)i(n/m)3LS(a,i5,i2)](n/m)3(„/m)s  =  ^nlm^nlmh  (86) 

(nlm)s 

A  number  of  important  relationships  can  be  proven  for  Q  and  R  including 


^(nlm)i(n/m)j(Q!’  P)  =  Q(nlm)t(nlm)i  (a>  ~R)  (87) 

P‘(nlm)l(nlm)i(ai  &■>  P)  =  R(nlm)i(nlm)i(0,a,—P)  (88) 

Q(n/m)i(n/m)j(^)  0)  =  ^(n/m)i(n/m)2  (89) 

^(nlm)i  (nlm)j  (®!  ®)  ^(n/m)i(nJra)2  (90) 

0<|<?(„/m)l(»,m)s(a,-R)|<l  (91) 

0  <  |-R(nIm)1(n/m)2(a^>-R)l  <  1  (92) 


A  number  of  other  relations  are  defined  by  Koga  et  al.14,27,28  for  the  Q-function  which 
they  call  the  S  function.  Corresponding  relationships  can  be  found  for  the  R  function. 

3.  MULTICENTER  INTEGRALS  OVER  COULOMB  STURMIANS 


The  addition  theorems  are  useful  for  application  to  problems  which  involve  the  eval¬ 
uation  of  multicenter  integrals  over  exponentially  decaying  orbitals.  Thus,  applications 
to  the  calulation  of  molecular  wavefunctions  for  discrete  and  continuum  processes  can  be 
considered.  In  particular,  the  general  4-center  electron  repulsion  integral,  defined  in  Eq. 
(1),  will  be  analytically  evaluated.  Now,  introducing  Eq.  (58)  into  Eq.  (1)  results  in 


a  b,  c,  d)  ^^jdpjz 


a+iPRcB 


p 2  ^  (93) 

*  P(nlm)a(nlm)i{®*'>  PABi  P  )  P(nlm)j(nlm)c(®di  ^ci  PDCi  P  ) 


where 
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(94) 


'  Rcb  =  C-B' 
i  Rab  =  A  —  P  ’ 

Rdc  =  D  —  C 

and  F  is  called  a  2-center  form  factor  and  is  given  by 

FinirnMnimhiatuatiiR;?  )  =  J  dr  S^h(a2, r)S^)h (<*1  ,r  +  P)e",pjr  (95) 
or,  by  introducing  the  appropriate  Fourier  transforms, 

i,,  («.,?)«-'*•'  (96) 

It  is  important  to  realize  that  £  in  Eq.  (93)  depends  only  on  the  differences  between 
intemuclear  vectors  (defined  in  Eq.  (94))  and  not  on  the  vectors  themselves.  Also,  the 
efficient  evaluation  of  the  2-center  form  factors  (F)  is  of  critical  importance. 

In  order  to  evaluate  the  form  factors,  the  2-center  overlap  of  two  CSs  in  coordinate 
and  momentum  space  needs  to  be  defined  and  evaluated.  Thus 


and 


(P,r)S^t).(a,r  +  R) 

(M^KP  +  P) 


It  should  be  clear  that 


(97) 


(98) 


(99) 


poS"U(“-a-5)  =  °SU 

where  ~  is  defined  in  Eq.  (41).  Also,  using  Eq.  (19) 


(100) 


^3  ri2  TI4 


E  E  E  E 

7»3=il+l  #5=ll+l  **4=f2+l  Hg  =  /2+l 


(101) 


TTl2 
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In  order  to  evaluate  and  p0^^Jn^  ,  the  translation  formulas  in  Eqs.  (60) 

and  (73),  along  with  Eq.  (101),  are  used  to  produce 


E  I, .,>(<«. rO”-"')(„;lj (a, ft 3)  (102) 

ni=max(l,n;—  1) 


rO(nl)i(nl)j(0liP'P)  ~  P(nlm)i(n1limj)(a»0>P)0(f/1|i)(n|)i 

ni=max(l,n,— 1) 

Then,  in  terms  of  Q,  the  form  factor  becomes 


(103) 


^P(n/m)i(n<m)2(f*li  a2i  P?  P)  —  ^  P^(n*)^(nl)2(al  ’  a2)  P)  (104) 

(n/m)  3 

where  a  pre-superscript  has  been  added  to  indicate  that  F  has  its  intemuclear  dependence 
expressed  in  terms  of  the  Q  function.  In  terms  of  R,  then 


P(f»Jm)i(n/m)2  (aa,a2  ;P;P)  =  ^  r0(^l,)^(3n{)a(ai,a1;^)fl(n/m)3(n/m)j(a1,a2,P) 

(nlm)  3 

(105) 

Using  Eqs.  (104)  and  (105)  in  (93)  results  in 


P  n—2  ^  >  Q(nlm)a(nlm)2(°la’  P-Ai?)  ^  ]  Q(nlm),j(n/m)i  (ad>  Rdc) 


(nlm)i 


(106) 


Q *(nlh(nl)!(nl)2(nl)b(ad'  a*’  ^Cb) 


*E  ~  9 Tt2  r^(nI)(,(n/)2(Q!“’ao>^fi)  ^  Pdc) 


(n/m)2 


(nim)i 


(107) 


Rhnl)iinl)2nlh(nl)h  (a<*’  aC’  aa’ab',  R-Cb) 


where 
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Q I(rd™(nl)!(nl)i(nl)t(ai'  a<: ’  a“ ’  a*’  )  ~ 


J  dp  po(-;,;“‘„,)t(aa,a4ip)  i8) 


,+iP^ci  Pn*m!' 


(108) 


fi^0iMMn02(n/)»(a<f,ac,Q!a,Q:6;  ^CB)  - 


J dP  p2e+  ~  CB  ■R(nlm)3(nlm)i(°‘a’  abi  P)R(nlm)i(nlm)c(ad>  aci  P) 


(109) 


It  is  then  easy  to  see  that 


Q  rmx mcm2 m*, 
J(n0i(n0c(n0j(n/)» 


'*1  I  * 

ER  imimtm}m( 
(nJ)i(»»3/c)(n 


n4=moi(l,nj  — 1)  nj=mai(l,ni  —  1) 
/imtmi  /\rncme 


(110) 


and  thus  only  RJ  needs  to  be  considered  at  this  point  in  the  development.  Thus,  using 
Eq.  (75)  in  Eq.  (109)  produces 

R^(nlh(nl)!(nl)i(nl)i (“<*’  a<=>  a“’  ^Cb)  = 

nj+ni  mtn(n4— 

E  E 

n4=maz(|n2-n»|,X)  i4=|{2-4|,2 

nt+rie  min(nj-l,li+lt) 

E  E 

n3=maz(|nl-nej,l)  /3=|h-/e|,2 

*  ^(^)3("),Km3  m,)(aJ  +  “c,«a  +  Kce) 

(111) 

where 


•7(”l)”n, ).(“»>“««)  5  J*Pjl'+iP  ii  )  (H2) 

In  order  to  perform  the  integral  in  Eq.  (112),  the  modified  CSFT  (\&)  must  be  ex¬ 
plicitly  given.  This  can  be  done  by  first  relating  the  CS  to  the  5-function  of  Steinbom8 
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(113) 


Sr,(a,r)  =  a3/2£<'B1?(a,r-) 


where 


nl  _  (~ l)n+Z  12/+5^2  r  n(n  + /)!  ii/2(— n  +  /  +  l)t-i(n  +  /  +  l)t-i 
fl<  “  (21  +  1)11  2(n  —  /  —  1)!  (t  —  1)!(/  +  3/2)t_i 

where  the  Pochmanner  symbol  is  defined  by 


(114) 


(115) 


Then  taking  the  Fourier  transform  of  Eq.  (113)  and  using  Eq,  (53),  the  Fourier  transform 


of  $  is  obtained  as 


.  o  1/2  ” _ j  pi 

*>•■?) = (-o'-  [a2 +P,r,yr(p) 


(116) 


Then,  J  can  be  evaluated  as 


*2=1  <1  =  1 

ll+l2 

x  X  ^i+i2+i3=e»en[(-l)'1iZl+'J+/3](^m2Kimi|Z3,m2  -  mi) 

i3=|h-*j| 


(117) 


where 


WkJ,n(<*U<*2>R)-  %JdP  Ja2  +  p2]-l[a2+P2]"  (U8) 

The  W  integral  may  be  evaluated  by  the  methods  developed  by  Harris  and  Michels.29 

Thus,  the  evaluation  of  E,  the  4-center  electron  repulsion  integral,  has  been  described. 
The  programming  of  these  results  is  underway. 
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ABSTRACT 

The  Schrodinger  equation  for  the  scattering  of  an  electron  by  a  hydrogen  molecule 
is  solved  by  converting  it  to  a  block-tri-diagonal  matrix  equation.  The  discrete  mesh  of 
points  in  the  angular  variable  must  be  made  more  dense  as  the  potential  energy  deviates 
from  spherical  symmetry  to  maintain  a  desirable  level  of  accuracy.  The  present  paper 
investigates  the  nature  of  this  dependence,  and  the  relative  effect  of  third-order  and  fifth- 
order  differentiation  rules.  It  is  found  that  the  existence  of  one  large  element  in  the 
K-matrix  sets  a  lower  bound  on  the  accuracy  of  all  elements.  This  bound,  and  therefore 
the  overall  accuracy  of  the  numerical  approximation,  can  be  substantially  improved  by  the 
use  of  the  higher-order  rules. 

*This  work  was  supported  by  NASA  contract  NAGW-2930  and  NCC  2-492,  and  an  Air 
Force  contract  27256-PH-SAH.  A  grant  of  computer  time  from  Florida  State  University  is 
gratefully  acknowledged. 
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1.  INTRODUCTION 


Traditional  electron  scattering  theory  is  concerned  with  the  results  of  collisions  be¬ 
tween  an  electron  and  a  target  under  the  assumption  that  their  interaction  depends  only 
on  their  separation,  and  is  independent  of  their  angular  orientation.  When  the  target  is 
an  atom,  this  asumption  is  justified,  but  when  the  target  atoms  are  bound  into  molecules, 
the  great  simplification  afforded  by  the  conservation  of  angular  momentum  are  no  longer 
available.  A  numerical  analysis  must  therefore  be  introduced  at  an  earlier  stage  in  the 
solution  of  the  problem. 

The  X- matrix  is  a  concise  and  computationally  simple  means  of  presenting  the  results 
of  a  scattering  calculation  using  the  method  of  partial  differential  equations.  This  work 
will  investigate  the  relationship  between  the  mesh  size  required  in  a  finite-difference  ap¬ 
proximation  to  the  -matrix  and  the  deviation  of  the  scattering  potential  from  spherical 
symmetry  for  various  levels  of  accuracy.  ‘Accuracy’  is  defined  as  the  difference  between 
the  numerical  approximation  and  a  hypothetical  continuum-limit  AT-matrix,  which  must 
be  inferred  from  the  behavior  of  the  approximation. 

The  angular  accuracy  is  roughly  independent  of  the  radial  properties  of  the  potential, 
so  the  numerical  requirements  of  the  polar-angle  dependence  can  be  investigated  without 
a  detailed  model  of  the  physical  potential.  The  angular  momentum  quantum  number  l 
provides  a  convenient  basis  for  calculating  the  AT-matrix,  but  one  which  is  not  a  basis  of 
eigenvectors.  In  an  eigenvector  basis,  the  AT-matrix  would  be  diagonal,  and  the  numerical 
convergence  of  the  various  elements  would  be  unrelated,  with  the  states  of  higher  rotational 
energy  requiring  a  finer  mesh.  Here,  by  contrast,  the  error  in  the  numerical  approximation 
is  spread  over  all  coupled  states.  In  particular,  if  one  element  dominates  the  others,  a 
small  relative  error  in  that  element  will  imply  a  large  relative  error  in  the  other  elements. 
The  convergence  of  the  largest  element  is  found  to  be  excellent,  with  accuracy  of  about 
10-3  using  29  points  in  the  angular  mesh.  The  convergence  of  the  smaller  elements  is 
limited  by  this  effectively  constant  bound.  Convergence  beyond  this  iimit  requires  vastly 
larger  meshes,  though  a  great  improvement  can  be  gained  by  the  use  of  higher-order 
differentiation  rules.  It  is  found  that  adding  rules  of  fifth  order  in  the  mesh  spacing 
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requires  little  additional  programming,  and  virtually  no  extra  computer  time  for  the  models 
considered. 


2.  THE  METHOD 

The  homonuclear  diatomic  molecule  is  the  simplest  system  for  which  non-spherical 
terms  contribute  to  the  scattering  potential,  so  it  will  serve  as  the  physical  anchor  for  an 
otherwise  computational  study.  When  the  energy  of  the  incident  electron  is  between  1  and 
10  electron  volts,  the  non-spherical  structure  of  the  molecule  is  particularly  important.  In 
such  a  situation,  we  use  the  axial  symmetry  of  the  molecule  to  reduce  the  Schrodinger 
equation  to  a  two-dimensional  partial  differential  equation.  We  assume  for  simplicity  that 
the  collision  does  not  excite  vibrational  states  of  the  molecule.  Rotational  excitations  may 
be  dealt  with  at  the  cross-section  stage  of  the  calculation  in  the  usual  adiabatic-nuclei 
manner.1 

The  target  molecule  may  be  assumed  to  be  oriented  in  a  particular  direction  in  space 
so  we  use  the  axis  of  symmetry  of  the  molecule  as  the  z-axis  of  our  coordinate  system  (Fig. 
1),  and  define  a  modified  wavefunction  u  which  is  related  to  the  true  wavefunction  for 
the  projectile  by 


*(r,«,*)  =  -rt<(r 
r 


(1) 


The  equation  for  u  in  atomic  units  is 


d2u 


+ 


1  f  ,, d2it  du 


}  +  {p2- 


2  V(r,6)- 


m 


'}“  = 


0 


(2) 


^2  ■  r2  '  dej  '  r2sin(0) 

where  V(r,  6)  is  the  scattering  potential,  and  p  is  the  momentum  of  the  projectile.  A 
finite  difference  approximation  to  this  equation  is  solved  on  the  domain  0  <  r  <  Rmax , 
0  <  0  £  f-  The  outer  limit  Rmax  is  chosen  to  be  sufficiently  large  that  the  scattered 
particle  is  effectively  free  at  that  distance.  The  range  of  6  has  been  cut  in  half  to  take 
advantage  of  the  symmetry  of  the  molecule  under  reflections  in  the  x  —  y  plane. 

The  boundary  condition  at  r  =  0  requires  that  u  =  0.  At  6  =  0,  there  are  two 
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possibilities:  if  m,  the  2-component  of  the  orbital  angular  momentum,  is  nonzero,  then 
tx  =  0;  otherwise  du/dd  =  0.  At  6  =  ir/2,  the  parity  of  the  system  dictates  either  that 
u  =  0  or  du/ d6  =  0.  At  r  =  Rmax,  we  impose  the  requirement  that  the  incoming  electron 
be  in  a  state  of  well-defined  orbital  angular  momentum,  so 


u(Rmaz,0)  =  Plm(cosd)  (3) 

The  parity  of  the  associated  Legendre  function,  (— l),+m,  determines  the  choice  of  bound¬ 
ary  condition  at  8  —  7r/2.  Though  this  parity  is  conserved  by  the  scattered  wave  function, 
the  value  of  l  for  the  incident  wave  is  not.  Outgoing  waves  therefore  contain  a  superposition 
of  l,  l  ±  2,  l  ±  4,  and  so  forth. 

Once  the  wave  function  has  been  found,  information  about  the  scattering  is  extracted 
by  matching  it  to  the  known  asymptotic  form 


*  =  £ 
v 


Awji'(pr)  +  Bu.yv(pr) 


PVm(cosO ), 


(4) 


where  l  and  m  are  the  index  and  order  of  the  Legendre  function  used  in  the  boundary 
condition  at  large  r,  and  j  and  y  axe  spherical  Bessel  functions.  The  coefficients  A  and  B 
form  matrices  on  the  space  of  allowed  values  of  orbited  angular  momentum. 

The  orthogonality  of  the  Legendre  functions  may  be  used  to  project  out  a  particular 
value  of  V .  When  this  is  done  for  two  different  large  values  of  r,  the  linear  independence  of 
the  Bessel  functions  allows  extraction  of  the  unknown  coefficients  A  and  B.  This  process 
must  be  done  for  each  allowed  value  of  then,  the  whole  solution  process  repeated, 
allowing  l  to  take  on  each  permitted  value.  The  result  is  two  matrices  describing  the 
transformation  of  a  angular  momentum  eigenstate  into  a  scattered  wave.  In  principle,  the 
dimension  of  the  matrices  is  infinite,  but  in  practice  the  finite  range  of  the  non-spherical 
potentials  provides  a  cutoff.  For  an  H2  molecule,  with  a  projectile  energy  below  10  eV, 
values  of  1  less  than  5  axe  all  that  are  required.  Hence,  the  matrices  are  at  most  3x3.  The 
K- matrix  is  defined  by  the  equation 


*  =  £ 

v 


jv(pr)  +  Kwyi'(pr) 


Pi>m(cos0), 


(5) 
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so  K  =  A-1B,  and  the  influence  of  the  artificial  imposition  of  an  initial  orbital  state  is 
removed.  The  diagonal  elements  of  K  are  the  tangents  of  the  elastic  phase  shifts,  and  the 
off-diagonal  elements  describe  the  coupling  of  states  of  different  orbital  angular  momentum 
by  the  potential. 


3.  THE  POTENTIALS 

The  potential  energy  function  V  in  Eq.  (2)  represents  the  electromagnetic  interaction 
of  the  projectile  with  the  target  molecule,  and  is  the  source  of  the  breaking  of  spherical 
symmetry.  In  principle,  this  interaction  is  highly  complicated.  A  molecule  is  by  no  means 
a  static,  rigid  body  in  the  presence  of  an  incoming  electron,  so  a  potential  model  is  not 
strictly  correct,  but  the  practical  usefulness  of  potentials  in  calculations  is  great  enough 
to  excuse  the  gross  simplification  of  the  scattering  problem  they  imply.  The  potentials 
involved  in  scattering  from  hydrogen  molecules  are  approximations  to  three  effects:  the 
static  electric  moments  of  the  molecule;  the  polarization  of  the  molecule  by  the  electric  field 
of  the  projectile;  and  the  substitution  of  an  electron  from  the  molecule  for  the  projectile 
in  the  outgoing  state,  known  as  ‘exchange.’  Toy  models  have  been  adopted  which  mimic 
to  some  extent  the  radial  structure  of  the  physical  potentials,  and  have  angular  properties 
which  are  expandable  in  the  low-order  Legendre  polynomials  of  even  parity. 

Because  the  hydrogen  molecule  is  composed  of  identical  atoms,  the  lowest-order  mul¬ 
tipole  field  seen  by  the  projectile  is  the  quadrupole  term 

VQ(r,0)  =  ^Mcose)  (6) 

where  Q  =  0.49  in  atomic  units.  This  form  is  valid  for  values  of  r  which  are  outside  the 
molecule  for  all  0.  At  shorter  ranges,  the  quadrupole  potential  is  cut  off  by  the  usual 
exponential  term  (1  —  exp{—(r/a)6])  as  a  multiplier.  The  adjustable  parameter  a  is  set  to 
1,  since  no  effort  is  made  to  fit  experimental  data.  The  static  potential  inside  the  molecule 
is  represented  by  a  Yukawa  function  of  unit  range  and  arbitrary  strength,  which  is  adjusted 
as  necessary  to  give  the  correct  order  of  magnitude  for  the  A-matrix. 
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Polarization  of  the  target  is  well  known  to  include  terms  of  both  spherical  and  P2 
character.1  The  form  used  here  is 

v>( r,0)  =  2^7  f1  _  e<'r)")  [°°  +  (7) 

with  a0  =5.5,  and  a2  =1.4. 

The  excha’  otential  is  a  crude  form  of  the  free  electron  gas  model.  The  basic  idea 
here  is  that  the  .ectron  in  the  target  may  be  approximated  by  a  degenerate  Fermi  gas 
whose  density  varies  with  position.1  The  exchange  potential  is  proportional  to  the  Fermi 
momentum  at  any  point  in  space,  giving 

Vx(r-)  =  -(|?)[3^(f)]  (8) 

The  electron  density  /?(r  )  is  approximated  by  a  simple  sum, 

/?(r,  9)  =  i-e-r  [1  +  1.3P2(cos0)]  ,  (9) 

which  has  contours  of  constant  density  which  resemble  those  of  the  sum  of  two  hydrogen 
atoms  separated  by  .7  angstroms.  The  resulting  potential  can  be  approximated  by  a  linear 
combination  of  Legendre  polynomials,  which  is 

Vx(r,e )  =  — 0.8e-r/3  [0.97  +  0A3P2(cos9 )  +  O.lP4(cos0)].  (10) 

The  ‘realistic’  potential  is  a  sum  of  these  three  terms. 

The  long-range  quadrupole  potential  is  the  most  important  term  in  real  scattering 
problems,  but  it  does  not  vanish  until  a  distance  of  14  Bohr  radii  from  the  molecule’s  center 
is  reached.  This  was  determined  a  posteriori ,  by  searching  for  a  region  in  which  the  K- 
matrix  does  depend  on  small  changes  in  Rma x.)  Such  a  long  range  requires  a  large  number 
of  radial  points,  and  consequently  a  long  running  time  for  the  computer  code.  Although 
the  convergence  properties  of  the  angular  approximation  are  roughly  independent  of  the 
quality  of  the  radial  approximation,  this  independence  is  only  valid  when  the  radial  mesh 
is  on  the  order  of  magnitude  of  the  mesh  which  gives  decent  radial  convergence.  Therefore, 
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the  quadrupole  term  is  generally  omitted  from  this  work,  in  the  interest  of  a  more  efficient 
use  of  computer  time.  The  P2  terms  in  the  polarization  and  exchange  potentials  are,  in 
any  case,  expected  to  yield  the  same  requirements  on  the  angular  mesh. 


4.  NUMERICAL  ANALYSIS 


The  usual  discrete  approximation  to  Eq.  (2)  is  obtained  by  replacing  the  radial  vari¬ 
able  with  a  set  of  Nr  points,  and  the  angle  with  a  set  of  N0  points.  The  wavefunction 
therefore  becomes  a  set  of  (real)  values  which  solve  the  equation2 


^2 


\\-r-cot(6j)(ui<j+1  -  uitj)  +  ~(ui)i+1  -  2uitj  +  «i,j-i)]  +  (11) 

/.■  fin  tin 


[P2  ~  2VS  j  - 


mM 


rfsin2(6j) 


Ki  =  0 


Here,  Uij  =  tt(rj,0j),  with  0  <  i  <  Nr  and  0  <  j  <  N0.  The  spacings  between  the 
radial  and  angular  points  are  hr  and  ha,  respectively.  (The  values  of  hr  and  h0  may  in 
principle  depend  on  i  and  j.)  Eq.  (11),  together  with  the  boundary  conditions,  gives  a 
system  of  (Nr  +  2 )(N„  +  2)  coupled  linear  equations  for  the  u^,  which  may  be  solved  by 
conventional  matrix  method.3 

The  choice  of  the  points  r*  and  6j  presents  an  unavoidable  conflict  of  interest.  There 
are  ways  to  improve  the  accuracy  of  a  numerical  solution  of  a  differential  equation  by 
judicious  choice  of  abscissae,  just  as  there  are  ways  to  improve  the  numerical  quadrature  by 
which  the  elements  of  the  if-matrix  Eire  projected  out  of  the  wave  function.  Unfortunately, 
the  ways  are  mutually  inconsistent.  Numerical  quadrature  is  optimized  by  choosing  points 
so  that  successive  errors  alternate  in  sign,  which  is  the  worst  possible  choice  of  points 
for  differentiation.  Therefore,  a  compromise  of  equally-spaced  abscissae  will  be  used,  and 
optimization  achieved  through  higher-order  integration  and  differentiation  rules. 

When  the  Schrodinger  equation  is  solved  on  a  two-dimensional  lattice,  the  resulting 
matrix  is  block  tridiagonEil.  This  property,  which  makes  the  matrix  easier  to  decompose 
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into  lower  and  upper  triangular  matrices,  depends  on  one  of  the  variables  having  a  simple 
three-point  derivative.  The  other  may  be  as  complicated  as  the  problem  demands.  In  this 
work,  the  radial  accuracy  is  of  secondary  interest,  so  the  ^-derivative  part  of  Eq.  (11), 
which  is  in  braces,  will  be  explored  with  rules  of  greater  accuracy  (Table  25.2  in  Ref.  4). 
For  a  point  in  the  interior  of  the  domain  of  approximation  of  0,  labelled  by  ‘j’  ,  rules  of 
order  h\  are 


%  =  5T*2*'2 " +  16/i+1 " 2*++  (12) 

W  =  Th(~  /j'2  +  “  30fi  + 16/2+1  ~ 2/j++  (13) 

where  D\  =  24h0  and  £>2  =  12 h2.  (To  simplify  the  formulae,  utiJ  has  been  replaced  by 
fj.)  For  points  next  to  the  boundary,  asymmetrical  rules  must  be  used: 

%  =  -6/i-i  - 20/y  +  36/,+,  -12/,+j  +  2/>+3},  (14) 

-Qfi  —  l  —  20/j  +  6fj+i  +4/>+2  -  /i+3 } ,  (15) 

These  rules  axe  used  for  points  next  to  the  boundary  at  6  —  0.  Next  to  8  =  7t/2,  the 
subscript  lj  +  3’  is  replaced  by  ‘j  —  3’  ,  lj  +  2’  by  ‘j  —  2’  and  so  forth,  and  there  is  an  overall 
minus  sign  in  the  first  derivative. 


At  the  boundary,  if  u  =  0  is  the  boundary  condition,  then  it  goes  directly  into  the 
matrix.  If  Neumann  boundary  conditions  apply,  the  process  is  more  involved.  First,  terms 
in  the  Schrodinger  equation  proportional  to  dr/ 69  are  dropped.  Then,  Eq.  (14)  is  used 
with  j  =  0  to  get  the  relation 


-6/_i  -  20/o  +  36A  -  12/2  +  2/3  =  0  (16) 

which  is  solved  for  f-\.  This  value  is  then  used  in  Eq.  (15)  (with  j  =  0  again),  to  derive 
an  expression  for  cPu/dQ2  at  0  =  0  which  is  accurate  to  fifth  order  in  hQ: 
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When  Eq.  (17)  is  used  in  place  of  the  second  derivative  term  in  Eq.  (11),  the  result 
is  an  approximation  which  significantly  improves  the  numerical  accuracy,  compared  to  the 
three-point  rules.  As  a  test,  when  a  spherical  potential  is  assumed,  and  the  7f -matrix 
calculated,  the  off-diagonal  matrix  elements  should  vanish.  In  fact,  for  a  given  number  of 
6  points,  they  are  at  most  one  fourth  of  those  calculated  with  three-point  rules,  and  of 
order  10-3  when  the  largest  diagonal  element  is  about  1. 

5.  RESULTS 


Convergence  of  the  If -matrix  is  determined  in  two  ways.  The  first  takes  advantage 
of  the  invariance  of  the  if-  matrix  under  time  reversal  to  require  that  the  If -matrix  be 
symmetric.  The  second  is  that  the  value  of  an  element  of  the  of  the  if -matrix  change  with 
changes  in  the  number  of  6  points  in  a  manner  which  points  clearly  to  an  asymptotic  limit. 

The  matrices  are  not  strictly  symmetrical.  When  the  scattering  potential  is  spherical, 
the  wave  function  should  have  the  angular  properties  of  the  large- R  boundary  condition, 
all  the  way  to  the  origin.  This  is  true  only  when  the  boundary  condition  is  for  l  =  0.  When 
higher  waves  are  used,  each  contains  a  contaminant  of  lower  1  values.  This  contamination 
decreases  with  each  step  in  r  away  from  the  origin,  reaching  an  eventual  limit  of  about 
.002  times  the  amplitude  of  the  correct  wave.  Therefore,  any  number  less  than  .002  in  the 
if -matrix  is  indistinguishable  from  zero.  This  seems  to  a  be  an  artifact  of  the  solution 
algorithm,  reducible  only  by  improving  the  J-derivative  rule.  (When  a  three-point  rule  is 
used,  the  corresponding  limit  is  .008.) 

The  second  criterion  makes  an  accurate  graphical  depiction  of  the  convergence  difficult. 
The  convergence  of  almost  all  elements  is  monotonic,  so  the  choice  of  a  specific  limit  is 
arbitrary,  yet  it  has  a  great  impact  on  the  perceived  quality  of  the  convergence.  For  this 
reason,  the  first  derivative  of  if//»  with  respect  to  Na  is  plotted  in  Figures  2-7. 

To  best  isolate  the  effects  of  various  non- spherical  potentials,  the  if-matrices  used  to 
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produce  Figures  2-7  are  derived  from  several  potentials.  First,  a  spherical  Yukawa  poten¬ 
tial  was  used.  Then  the  same  Yukawa  potential  multiplied  by  the  Legendre  polynomial 
P2(cos0).  Then  half  of  each  is  used.  The  length  scale  of  the  potential  was  3,  and  its 
strength  was  set  to  -0.4.  The  momentum  of  the  incoming  particle  was  0.8,  corresponding 
to  a  kinetic  energy  of  8.7  eV.  All  these  parameters  were  chosen  to  give  significant  K -matrix 
elements  up  to  1  =  2.  An  exception  is  the  P2  potential  in  Fig.  5.  The  1  =  2  element  was 
very  small,  so  its  relative  error  was  extremely  large,  and  not  indicative  of  the  numerical 
properties  of  the  algorithm.  The  fourth  curve  on  the  graphs  is  the  ‘realistic’  potential 
described  above,  with  Yukawa  of  unit  range  in  place  of  the  long-range  quadrupole. 

Figures  2  and  3  show  the  slope  versus  N0  of  the  if- matrix  element  for  transitions 
between  spherical  waves.  Fig.  2  is  for  a  three-point  rule,  and  Fig.  3  is  for  the  five-point 
rules.  The  advantage  of  the  five-point  rule  is  most  obvious  at  small  numbers  of  theta  points. 
As  the  mesh  becomes  finer,  the  relative  advantage  of  the  five-point  rule  diminishes,  but 
it  is  still  a  factor  of  two  better.  Figures  4  and  5  are  the  same,  for  the  d-wave  to  d-wave 
reaction. 

Figures  6  and  7  are  different.  They  show  the  difference  between  the  1  =  0  to  2  matrix 
element  and  the  2  to  0  element-this  provide  a  test  of  time-reversal.  The  difference  is 
normalized  to  the  value  of  K02 ,  usually  the  smaller  of  the  two.  In  judging  the  size  of  the 
relative  error,  it  should  be  borne  in  mind  that  there  is  always  a  diagonal  element  of  the 
if -matrix  which  dominates  this  term  by  an  order  of  magnitude.  The  uncertainty  in  this 
dominant  element  is  the  primary  source  of  this  error.  In  this  place,  the  five-point  rule 
shows  its  particular  virtue.  The  five-point  rule  shows  a  more  symmetrical  if -matrix  with 
11  points  in  the  mesh  than  the  three-point  rule  with  30. 

6.  CONCLUSIONS 

The  addition  of  five-point  approximations  to  the  angular  derivatives  requires  only  mi¬ 
nor  changes  in  the  computer  code  for  solving  partial  differential  equations.  It  preserves  the 
block-tridiagonal  character  of  the  equivalent  linear  system,  so  it  requires  no  extra  running 
time,  and  it  yields  significant  improvement  in  the  accuracy  of  the  numerical  solution. 
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The  chief  drawback  of  the  partial-differential-equation  method  for  non-spherical  scat¬ 
tering  problems  is  that  the  coupling  between  states  makes  the  relative  error  in  the  small 
elements  of  the  if -matrix  large,  since  it  is  a  small  fraction  of  a  large  number.  The  largest 
element  of  the  if -matrix  thus  creates  a  ‘floor’,  a  lower  bound  on  the  possible  error  in  the 
other  elements. 

The  contamination  of  the  wave  function  by  lower  partied  waves  them  the  input  bound¬ 
ary  condition  is  the  largest  single  source  of  error.  It  is  an  artifact  of  the  u  =  0  boundary 
condition  at  the  origin.  The  natural  tendency  of  this  algorithm  is  to  assume  the  smoothest 
possible  wave  function  near  r  =  0,  and  it  takes  several  steps  way  from  the  boundary  before 
the  large-r  boundary  condition  can  assert  its  dominance.  A  more  sophisticated  treatment 
of  the  radial  derivatives  might  well  be  rewarded  by  a  reduction  of  this  defect. 

Though  this  work  has  confined  itself  to  the  angular  behavior  of  the  wave  function,  the 
great  improvements  afforded  by  the  higher-order  approximations  to  the  derivatives  might 
have  counterparts  in  the  radial  direction.  This  would  destroy  the  tridiagonality  of  the 
matrix,  but  in  applications  where  the  size  of  the  computer’s  memory  is  less  of  a  constraint 
than  the  cost  of  processor  time,  it  may  be  worthwhile. 
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Figure  1.  The  coordinate  system  for  scattering  from  a  diatomic  molecule. 
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Figure  2.  The  relative  change  in  the  if-matrix  element  for  5-wave— *s-wave  transitions, 
using  three-point  derivative  rules:  □,  uses  spherical  potential,  with  a  spherical  angular 
term  Po;  •,  uses  only  term  with  P2,  the  second  Legendre  polynomial;  o,  uses  the  most 
realistic  potential  (see  text);  o,  uses  Pq  and  P 2  terms. 


Number  of  0  Points  (Ne) 


Figure  3.  The  relative  change  in  the  jFf-matrix  element  for  s-wave— ►s-wave  transitions, 
using  five-point  derivative  rules:  n,  uses  spherical  potential,  with  a  spherical  angular  term 
p0;  •,  uses  only  term  with  P2,  the  second  Legendre  polynomial;  o,  uses  the  most  realistic 
potential  (see  text);  o,  uses  P0  and  P2  terms. 
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Figure  4.  The  relative  change  in  the  if -matrix  element  for  d- wave— ><i- wave  transitions, 
using  three-point  derivative  rules:  □,  uses  spherical  potential,  with  a  spherical  angular  term 
Po;  *,  uses  only  term  with  P2,  the  second  Legendre  polynomial;  0,  uses  the  most  realistic 
potential  (see  text);  o,  uses  Pq  and  P2  terms. 
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Figure  5.  The  relative  change  in  the  if-matrix  element  for  d- wave  — > d- wave  transitions, 
using  five-point  derivative  rules:  o,  uses  spherical  potential,  with  a  spherical  angular  term 
Po;  o,  uses  the  most  realistic  potential  (see  text);  o,  uses  P0  and  P2  terms. 
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Figure  6.  The  relative  change  in  the  -matrix  element  for  s- wave— *d- wave  transitions, 
using  three-point  derivative  rules:  □,  uses  spherical  potential,  with  a  spherical  angular 
term  Po ;  •,  uses  only  term  with  P2,  the  second  Legendre  polynomial;  o,  uses  the  most 
realistic  potential  (see  text);  o,  uses  P0  and  P2  terms. 
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Figure  7.  The  relative  change  in  the  if -matrix  element  for  5-wave— xi-wave  transitions, 
using  five-point  derivative  rules:  □,  uses  spherical  potential,  with  a  spherical  angular  term 
Po;  •,  uses  only  term  with  P2,  the  second  Legendre  polynomial;  o,  uses  the  most  realistic 
potential  (see  text);  o,  uses  Pq  and  P2  terms. 
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Abstract 

A  calculation  of  e-N2  scattering  in  the  vicinity  of  the  2.4  eV  (II9)  resonance  has  been 
completed.  The  main  element  of  the  calculation  is  a  15  term  vibrational  close  coupling 
expansion,  reduced  to  coupled  two-dimensional  (2d)  partial  differential  equations  (pde’s), 
and  solved  using  the  noniterative  pde  technique.  The  potential  consists  of  static,  exchange, 
and  polarization  parts;  each  part  has  been  (previously)  derived  in  a  manner  appropriate  to 
its  importance  in  the  scattering  equation.  Results  for  the  absolute  total  cross  section,  both 
in  magnitude  and  shape  of  the  substructure  in  the  resonance  region  (1.5<  k2<3.0eV),  are 
in  excellent  accord  with  experiment.  Angular  distributions  are  also  calculated  and  found 
to  vary  significantly  in  shape  in  the  immediate  vicinity  of  the  center  of  the  resonance 
(2.05<k2<2.15eV),  indicating  the  need  for  differential  measurements  at  a  finer  grid  in 
energy,  and  therefore,  requiring  even  better  energy  resolution.  Comparison  with  other 
calculations  and  discussion  of  some  theoretical  aspects  are  also  included. 


PACS  34.80.Gs,  34.80.Bm 
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1.  INTRODUCTION  , 


The  purpose  of  this  paper  is  to  report  a  calculation  of  electron  scattering  from  molec¬ 
ular  nitrogen  (e-N2)  using  the  noniterative  partial  differential  equation  technique.1  This 
calculation  is  confined  to  the  energy  region  around  the  2.4  eV  resonance;  it  consists  of  a 
15-state  vibrational  close  coupling  expansion  for  the  resonant  (IIS)  partial  wave  combined 
with  fixed-(plus  adiabatic-)nuclei  amplitudes  for  the  other  (contributing)  partial  waves. 
The  underlying  method  is  the  hybrid  theory.2  We  shall  concentrate,  on  the  angular  distri¬ 
butions  for  the  lower  vibrational  states  (u  =  0, 1, 2),  showing  that  their  shapes  change  very 
rapidly  from  one  energy  to  the  next  as  one  traverses  the  energy  region  in  the  center  of 
the  resonance.  Thus,  despite  the  fact  that  several  differential  scattering  experiments  have 
been  reported  (to  be  discussed  below),  they  do  not  have  the  energy  resolution,  nor  have 
they  attempted  to  explore  the  variation  in  the  angular  distributions  in  this  narrow  energy 
region  in  detail.  It  is  one  of  our  aims  to  motivate  such  experimented  investigations. 

Several  important  developments  in  our  calculations  have  been  made  since  the  intro¬ 
duction  of  the  hybrid  theory2,  which  we  have  reported  piecemeal  since  then.  Briefly,  the 
most  salient  of  them  are:  the  noniterative  technique  itself1 ,  the  reduction  of  the  scattering 
equation  to  2d  form,  which  was  first  outlined  in  Ref.  3,  in  addition  to  which — it  turns  out — 
the  2d  technique  could  also  be  applied  to  the  derivation  of  the  polarization  potential4,  as 
well  as  to  a  method  for  exactly  including  exchange  in  the  static-exchange  approximation.5 
And  finally,  the  direct  (i.e.  non-exchange)  static  potential  could  be  calculated6  using  a 
much  better  (MCSCF)  approximation  of  the  N2  ground  state  than  the  SCF  approxima¬ 
tion,  used  originally.2-5  Some  interim  results  using  a  10  state  vibrational  expansion  were 
reported  in  Ref.  7. 


2.  THEORETICAL  ASPECTS 


We  start  with  an  antisymmetrized  ansatz  for  the  total  wavefunction  of  the  e-N2  system 

2,7. 


15 

*(m)  =  £(-!)"  R^n^-r)  (2.i) 

i=i 

Here  X{  are  the  coordinates  (space  and  spin)  of  the  ith  electron  and  is  the  collection  of 
coordinates  of  the  remaining  (fourteen)  electrons.  The  factor  (— l)Pi  is  the  parity  of  a  cyclic 
permutation  (P{)  of  the  sequence  1,2,...,  15,  thus  making  completely  antisymmetric. 
We  first  summarize  the  analysis,  mostly  of  Ref.  5,  whereby  the  basic  integro-pde  is 
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derived.  (Note,  however  that  we  are  here  using  Rydberg  units  as  opposed  to  atomic  units, 
used  in  Refs.  [5,6]).  Insertion  of  Eq.  2.1  into  the  Schrodinger  equation,  and  premulti- 
plication  by  the  target  ground  state,  $at2,  yields  the  static  exchange  approximation.  We 
first  consider  the  ground  state  (£+ )  to  be  represented  by  a  single  determinant  (i.e.  SCF 
approximation): 


=  det{la2g2a2gla2u2a2uZa2uU2xuU2yu\ £+ )  (2.2) 

Labelling  the  different  orbitals  <pa  [a  :  1<t?,  •  •  • ,  l7ryu],  we  recognize  that  each  <pQ  is  an 
explicit  function  of  the  coordinates  of  a  single  electron  and  an  implicit  function  of  the 
interauclear  separation  R: 


<f>a  =  ,  •  •  • ;  -R)Ca(spm)  (2.3) 

When  one  includes  the  R  dependence  of  4>Q  and  F^m\  one  derives7  a  3d  pde  for  F 

7 

[-V2  +  B.a(R)  +  V(r; R)  -  E«]  R)  =  2  £  W^\f-  «)*„(?;  R),  (2.4) 

Q=  1 

where  V(r;R )  is  (to  begin  with)  the  static  potential  between  the  scattered  electron  and 
the  target  (N2): 


Vstatic  (r;R)  —  ($iv2|V’e_moi!$N2)  (2.5) 

and  WQ  are  the  (static)  exchange  kernels 

‘(nR)  =  J  <jV«(f';B)p-H_^  *<"*>(?';*)  (2.6) 

The  energy  appearing  in  Eqs.  2.4  and  2.5  is 

Esc  =  +  eV5  (2-7) 

and  Hvh(R)  is  the  vibrational  part  of  the  target  Hamiltonian  from  whose  intemuclear 
potential,  vibrational  wave  functions  ( Xv )  and  energies  (ev)  are  evaluated.  (In  this  calcula¬ 
tion,  the  latter  sure  obtained  numerically  from  the  MCSCF  potential  energy  curve  reported 
in  Ref.  6. 

For  the  static  potential  (2.5),  we  use  an  MCSCF  approximation6  of  An  MCSCF 
wave  function  is  a  sum  of  determinants,  whereas  the  right  hand  side  (i.e.  the  exchange 
part)  of  the  scattering  Eq.  2.4  assumes  that  $jv2  is  a  single  determinant  (i.e.  SCF 
approximation).  Thus  (2.4)  constitutes  an  admittedly  unbalanced  approximation,  which 
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nevertheless  seemed  well  justified6,7  because  (a)  the  exchange  terms  usually  have  a  quanti¬ 
tatively  smaller  effect  than  the  direct  terms,  and  (b)  when  written  as  in  (2.4),  the  exchange 
terms  can  be  rigorously  reduced  to  a  coupled  set  of  ordinary  (i.e.  non- integral)  pde’s.6  This 
comes  about  by  expanding  and  W'i™'*  in  vibrational  states  Xv  of  iV2 

Nv 

F<n|(r;J!)  =  5]F<")(?)x.W 

(2'8> 

wir\p,  r) = E 

v=0 

and,  using  the  well  known  property  of  the  Coulomb  potential  as  the  Green’s  function  of 
the  kinetic  energy, 

V^pr^Tj)  =  -to«(f'  -  f),  (2.9) 

allows  the  scattering  equation  (2.4)  to  be  reduced  to  coupled,  but  nonintegral  pde’s7: 


[v2  +  ki]  Fim\r)  =  £  vvA^r\^-^<f>[aj(r)wira)(n  , 

V=0  L  oral 

Nv 

= -to  y, 


(2.10) 


The  double  indexed  quantities  in  (2.10)  axe  vibrational  matrix  elements  of  the  unsub- 
scripted  quantities.  Thus  the  matrix  elements  in  Eqs.  2.10  are  given  by  the  following 


expressions: 


V»y(f)=(x.(H)|V(f;Ji)|x„.(fi)), 

*{‘M  =  (x.(W«(r-;B)lxX*)>. 

n  •  ' 


(2.11) 


Eqs.  2.10  can  be  further  reduced  to  2d  pde’s  by  exploiting  the  cylindrical  symmetry  of  the 
various  functions.  As  derived  in  Ref.  5,  they  take  the  form 


[A(m)  +  k2v]  fim){z)  =  I  Vvv.{z  ,)/iri)(z)  -  i  ^2  4>{${z)w{™'0\z)  , 

v'=0  L  r  oral 

A(m  -  m„)«,<"--«>(i)  =  E  «©U)/‘,m)(z). 


(2.12) 
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here  z  —  (r,  9)  and  A (m)  is  the  2d  Laplacian 


A  (m) 


&  J_  1  d  .  3_  _  m2 

dr 2  r2  502  sin6  dQ  dd  sin2  9 


(2.13) 


As  stated  above  Vvv>  should  be  derived  from  the  static  potential  seen  by  the  scattered 
electron,  given  by  Eq.  2.5.  We  have  shown6  that  the  use  of  an  MCSCF  in  (2.14) 
leads  to  a  slightly  less  attractive  potential  than  one  obtained  from  an  SCF  approximation. 
That  result  is  in  qualitative  accord  with  that  found  by  Rumble  et  al.8  Very  recently, 
however,  Meyer  et  al.9  have  carefully  examined  the  IIS  resonance  in  the  static  exchange 
approximation  and  have  found  that  the  resonance,  using  an  MCSCF  ground  state,  when  the 
consistent  exchange  terms  are  included,  corresponds  to  a  slightly  more  attractive  effect 
than  the  corresponding  consistent  exchange  approximation  with  an  SCF  ground  state. 
Their  conclusion  (which  we  find  quite  surprising,  but  do  not  question)  is  that  the  use  of 
SCF  orbitals  for  exchange  in  Eqs.  2.10  and  2.11  is  sufficiently  inconsistent  with  the  use  of 
an  MCSCF  wave  function  in  calculating  the  direct  term,  that  it  gives  the  opposite  effect 
from  what  is  obtained  by  using  the  MCSCF  wavefunction.9 

Nevertheless,  with  regard  to  our  present  calculation,  we  believe  this  effect  is  included 
by  the  way  polarization  has  been  incorporated.  To  our  static  potential,  we  add  a  polar¬ 
ization  potential,  as  described  in  Ref.  6: 


V(r: ;  R)  Vstatic(r;  R )  +  Vpol(r: ;  R)  (2.14) 

where 

R)  =  (l  -  e-W'->!)  v£P(r-,R)  (2.15) 

Vpo/  ’  is  calculated  from  a  quasi-  ab  initio  polarized  orbital  derivation  of  the  polarization 
potential4,  but,  as  seen  in  (2.15),  the  latter  is  diminished  in  magnitude  by  a  tuning  factor, 
(1  —  exp[(r/r0)]2),  in  which  r0  was  adjusted  to  give  the  IIff  resonance  at  the  correct  (i.e. 
experimental)  energy. 

We  found  specifically  in  Ref.  6  that  a  value  of  r0  =  2.430  was  required  when  using 
V8tatic(MCSCF),  whereas  ro  =  2.934  was  required  for  an  SCF  wavefunction.  (There  is 
an  unfortunate  typographical  error  in  Ref.  6,  which  read  ro  =  2.394.)  What  this  implies, 
in  agreement  with  Meyer  et  al.9,  is  that  although  (symbolically) 

(Static  d“  Vtxch)\fcsCF  ^  (Static  "h  Vexch) SCF  ^  (2.16) 

because 


(VPol) 


MCSCF 


>  OM 


SCF  ’ 


(2.17) 
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it  is  perfectly  possible  that 


(y static  "l"  Vcxch  4"  Vpo/) MCSCF  ^  (^static  4"  Vtxch  4"  Vpol) SCF  ’  (2.18) 

[To  repeat,  the  above  equations  culminating  in  Eq.  2.18  are  intended  as  a  measure  of  “size” 
in  an  average  sense;  in  detail,  both  the  r  and  R  dependence  of  each  side  will  be  noticeably 
different  from  the  other.  We  expect  therefore  that  the  left  hand  side,  which  derives  from 
the  better  (MCSCF)  wave  function,  does  represent  a  considerable  improvement,  and  that 
is  what  is  used  here.] 

Thus,  from  the  pragmatic  point  of  view  of  this  calculation,  we  consider  the  criticism 
of  Ref.  9  to  have  been  overcome;  in  fact,  we  believe  that  the  results  we  shall  present  axe, 
in  the  region  of  the  center  of  the  resonance  (1.6  <  k2  <  3eV  and  v  <  2),  the  most  accurate 
that  have  thus  far  been  calculated. 

Moreover,  the  conclusion  of  Ref.  9  notwithstanding,  it  is  completely  possible  that 
if  one  had  included  the  effect  of  static  polarization,  exchange,  and  correlation  simultane¬ 
ously,  (part  of  such  effects  would  be  described  by  the  exchange-polarization  terms  in  a  full 
polarized  orbital  treatment10,  for  example),  and  if  one  had  then  isolated  the  piece  labelled 
Static  +  Vexch  separately,  then  it  might  have  led  to  the  result  of  Weatherford  et  al.6,  that 
the  MCSCF  result  would  be  less  attractive  than  the  corresponding  SCF  result  [i.e.  the 
reverse  of  (2.16)]. 


3.  RESULTS  AND  COMPARISONS 

To  the  one  resonant  partial  wave  (Jlg),  we  add  the  four  most  important  non-resonant 
partial  waves  (1Eg,1  II,,,1  Aj),  calculated  in  the  fixed-plus  adiabatic-nuclei  approxi¬ 

mations,  as  discussed  in  Ref.  6.  Let  us  first  show  and  discuss  the  totcil  cross  section  ax 
(the  sum  of  all  energetically  allowed  vibrational  channels,  summed  and  averaged  over  ro¬ 
tational  states  in  the  usual  way).  The  resonance  with  its  famous  substructure11,  is  usually 
compared,  as  it  is  here,  with  the  experimental  result  of  Kennerly.12  The  theoretical  curves 
shown  are  our  previous  10  state  result7  (dashed  curve),  our  present  15  state  result  (solid 
curve),  compared  to  experiment  and  the  Schwinger  multichannel  calculation  of  Huo  et 

al.13’14 

With  reference  to  our  calculations,  the  comparison  of  the  10  and  15  state  results  gives 
a  good  idea  of  the  convergence  of  the  close  coupling  expansions:  we  would  say  that  our 
calculations  are  well  converged  to  just  beyond  the  first  resonance  (k2  <  1.95eV),  reasonbly 
well  converged  to  just  beyond  the  second  (~  2.2eV),  and  approximately  converged  to 
«  2.5eV.  Of  particular  note,  therefore,  is  the  fact  that  ours  is  the  only  calculation  which 
describes  the  magnitude  of  the  first  peak  (ax  ~  27d2  at  k2  =  1.95eV)  and  gives  the  ratio  of 
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the  first  two  peaks  accurately.  The  calculation  of  Huo  et  al.14  is  almost  as  satisfactory,  while 
their  adjoining  paper13  gives  the  whole  sequence  of  vibrational  excitation  cross  sections  in 
remarkable  accord  with  experiment.15 

The  ability  of  such  theories  as  in  Refs.  13  and  14  to  achieve  such  an  elaborate 
overview  of  the  entire  resonance  substructure  goes  back  to  the  physical  ideas  underlying 
the  boomerang  model  (cf.  Ref.  17),  which  have  been  given  their  most  rigorous  justification 
in  the  R-matrix  theory  and  calculations  of  Schneider,  LeDourneuf,  and  VoKyLan.18  More 
discussion  of  these  theories  will  be  included  in  the  latter  parts  of  this  paper. 

In  this  paper,  we  shall  concentrate  on  angular  distributions.  In  Fig.  2,  we  show  the 
elastic  differential  cross  section  at  1.50  eV,  just  below  the  onset  of  the  II9  resonance.  Also 
shown  are  various  other  theoretical2,14’19  and  experimental20-22  results.  Both  similarities 
and  differences  are  evident.  The  results  continue  into  the  heart  of  the  resonance  region 
( k 2  «  2.1eV),  in  Figs.  3  and  4.  We  have  divided  those  results  into  two  parts:  in  Fig. 
3,  we  show  our  present  results  at  three  energies  surrounding  2.10  eV;  one  sees  how  sig¬ 
nificantly  the  shape  varies  over  0.1  eV.  This  is  particularly  relevant  to  the  fact  that  the 
latest  experiment22  only  claims  an  energy  uncertainty  of  just  that  amount.  Nevertheless, 
the  discrepancy  between  all  these  results  and  the  original  hybrid  theory  calculation2  is 
clear.  In  Fig.  4,  the  same  experimental  information  is  compared  to  the  present  calculation 
and  that  of  Huo  et  al.14  at  2.10  eV.  Here  the  similarity  of  the  calculated  results,  both  of 
which  provide  absolute  values,  is  the  most  striking  feature.  It  should  be  noted  that  an 
SCF  target  representation  was  used  in  Ref.  14.  The  agreement  between  theories  suggests 
that  the  recent  experiment22  is  dominated  by  the  particular  energy  in  the  composite  beam 
which  gives  the  dominant  cross  section  at  a  particular  angle. 

The  comparison  of  experiment  and  theoty  at  3.0  eV  is  shown  in  Fig.  5.  The  similarity 
of  the  recently  calculated  results  continues,  but  —  referring  back  to  Fig.  1  —  we  emphasize 
that  at  3  eV,  the  energy  is  definitely  pressing  the  outer  edge  of  reliability  of  the  present 
calculation.  The  other  interesting  feature  of  the  latest  experiment22,  noted  there  as  well,  is 
that  the  absolute  magnitude  in  the  directions  6  <  110°  favors  the  original  hybrid  result.2 

In  the  next  three  figures  (6,7,8),  we  show  differential  cross  sections  in  the  excitation 
transition  to  the  first  vibrational  state.  At  the  lowest  energy  (1.50  eV),  our  present  calcu¬ 
lation  (Fig.  6)  is  definitely  favored  by  the  experiment  of  Brennan  et  al.22,  whereas,  in  the 
vicinity  of  the  dominant  peak  (Fig.  7),  that  is  only  true  in  the  middle  of  the  angular  range. 
Note  that  a  different  set  of  experimental  results  from  an  Australian  group  has  also  been 
included.23  At  the  highest  energy,  3  eV,  the  graph  also  includes  results  of  an  R-matrix 
calculation24  and  yet  another  experiment.25  Here  the  similarity  of  the  present  results  with 
those  of  the  R-matrix  calculation  is  the  most  notable  feature,  particulary  at  middle  Jingles, 
where  they  agree  best  with  the  experimental  results  of  Ref.  23. 

Finally  we  show  in  Fig.  9,  the  angular  distribution  associated  with  the  excitation  of  the 
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second  vibrational  state  in  the  vicinity  of  the  dominant  peak.  Here  the  two  experimental 
results22,23  are  in  agreement  with  each  other  to  within  their  experimental  errors,  and  they 
agree  best  with  our  present  result  at  the  experimental  energy  2.10  eV,  rather  than  at  its 
fringes  (2.05,  2.15  eV). 


4.  DISCUSSION 


This  completes  the  presentation  of  the  results  of  this  calculation.  (More  results  are 
available  upon  request.)  Although  the  method  (the  hybrid  theory)  is  ab  initio  in  principle, 
it  contains  here  one  phenomenological  parameter,  the  polarization  cut-off  r0  [cf.  Eq.  2.15]. 
The  main  practical  purpose  for  doing  this  is  to  provide  the  best  cross  sections  for  (several 
space)  applications.  It  will  be  recalled  that  our  (hybrid)  theory  was  developed  in  the 
context  of  SAR  arcs.2  The  numerical  results  of  that  calculation  were  collected  as  a  NASA 
document.26  The  present  (more  accurate)  cross  sections  are  intended  for  the  understanding 
of  secondary  electron  flux  in  the  lower  F  region  of  the  ionosphere.  The  specific  question 
concerns  whether  or  not  there  is  a  dip  in  the  electron  distribution  function.27  Relevant 
ionospheric  calculations  are  now  ongoing  at  Goddard  Space  Flight  Center28,  using  our 
cross  sections  and  those  of  others,  principally  Huo  et  al.13’14 

From  a  more  fundamental  point  of  view,  the  present  results  are  intended  to  be  a  more 
definitive  comparison  with  experiment  in  the  resonance  region,  particularly  with  angular 
distributions.  As  we  have  seen  from  the  comparisons  in  the  previous  section,  there  is  still 
insufficient  agreement  among  experiments  themselves  to  provide  a  definitive  check  at  this 
time.  In  addition,  future  experiments  will  require  an  even  finer  energy  resolution,  at  a 
finer  energy  grid,  to  be  compelling  in  this  regard.  This  having  been  said,  it  is  important 
to  acknowledge  that  great  progress  in  experimental  angular  distributions,  culminating  in 
the  recent  work  of  the  Australian  group22,23,  has  already  been  made. 

With  regard  to  calculational  methodology,  specifically  hybrid  vs.  R-matrix  and/or 
other  L2  basis  set  theories,  in  addition  to  what  was  said  above  and  elsewhere29,  it  is  clear 
that  the  latter  are  capable  of  giving  the  greater  overall  accuracy  at  the  present  time,  as  is 
exemplified  by  the  results  of  Refs.  (13,14,18).  (In  further  detail,  for  example,  the  R-matrix 
theory  shows  that  the  II5  resonance  is  dominated  by  short  range  correlations  rather  than 
the  long  range  polarizability.)  However,  if  insufficient  correlation  is  included,  the  R-matrix 
method  and  L2  methods  in  general,  can  yield  noticeable  inaccuracy  (cf.  Ref.  30,  for 
example). 

Furthermore,  the  basic  tenet  of  the  hybrid  theory  is  common  to  both  approaches;2,29 
it  is  the  fact  that  if  the  interaction  time  of  the  resonance  is  comparable  to  vibrational  time 
scales  of  the  target  molecule,  but  short  compared  to  rotation  time  scales,  then  in  one  way 
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or  another,  a  dynamical  treatment  of  the  N  + 1  electron  system  is  required.  That  treatment 
can  be  either  a  dynamical  coupling  of  the  incident  particle  with  the  vibrational  motion  of 
the  target,  or  a  recalculation  of  the  Bom-Oppenheimer  problem  of  the  ( N  4-  l)-electron 
system,  followed  by  a  calculation  of  the  vibrational  spectrum  in  the  (necessarily  complex) 
potential  energy  well  of  the  compound  system. 

The  drawback  of  the  present  methodology  is  due  to  the  slow  convergence  of  the  vi¬ 
brational  close  coupling  expansion.  In  principle,  that  can  be  overcome  by  going  to  a  3d 
approach  in  which  the  intemuclear  separation  (R)  becomes  the  third  dynamical  variable. 
Such  a  theory  has  already  been  outlined,31  and  the  non-iterative  pde  method1 2 3 4 5 6  has  been 
generalized  to  three  and  higher  dimensional  equations  and  applied  to  a  solvable  model.32 
The  implementation  of  this  program  is  already  in  progress,  however,  its  completion  will 
not  be  easy,  in  particular,  the  detailed  substructure  that  such  calculations  reveal,  will 
depend  critically  on  the  range  and  mesh  size  with  which  the  R  variable  can  be  covered. 

Perhaps  the  most  exciting  potential  application  of  the  3d  pde  technique,  also  dis¬ 
cussed  in  Ref.  32,  is  the  fact  that  it  can  be  applied  to  scattering  (in  principle)  from 
arbitrary  polyatomics  in  the  fixed-nuclei  approximation.33  Augmented  by  the  adiabatic- 
nuclei  approximation,33  such  a  method  would  be  an  invaluable  tool  in  studying  scattering 
processes  in  galactic  environments,  such  as  the  Orion  nebula,  where  it  is  known  that  exotic 
molecules  can  form  (cf.  Ref.  34,  for  example)  which  may  not  be  amenable  to  laboratory 
experimentation,  so  that  theoretical  calculation  provides  the  only  reasonable  alternative. 
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Fig.  1.  Total  integrated  cross  section  in  square  Angstroms  (A2):  Experiment:  (*  *  *) 

Ref.  12.  Theory:  ( - )  present  [15  state  results];  ( - )  our  10  state  results  (Ref. 

7);  (-o-o-o-)  Ref.  14. 
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e‘+N2  :  v(0-->0) 

DIFFERENTIAL  C/S 
1.50  eV 


Fig.  2.  dcroo/df2  at  1.50  eV:  Theory:  ( - )  present;  (-o-o-o-)  Ref.  14;  ( - )  Ref. 

2;  ( - )  Ref.  19.  Experiment:  ($$$)  Ref.  22;  (+++)  Ref.  20;  (AAA)  Ref.  21. 
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e"+N2 :  v(0— >0) 

DIFFERENTIAL  C/S 
2.05,2.10,  2.15  eV 


Fig.  3.  d<r0o/dJ2.  Theory:  ( - )  present  at  2.05  eV;  ( - )  present  at  2.10  eV; 

( - — )  present  at  2.15  eV;  ( - )  Ref.  2.  Experiment:  ($$$)  Ref.  22;  (+++) 

Ref.  20. 
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e"+N2  :  v(0~>0) 
2.10  eV 


Fig.  4.  dcroo/dfi  at  2.10  eV.  Theory:  ( - )  present;  (-o-o-o-)  Ref.  14.  Experiment: 

($$$)  Ref.  22;  (+++)  Ref.  20. 
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e"+N2  :  v(0— >0) 
3.00  eV 


Fig.  5.  d<j00/dO  at  3.00  eV.  Theory:  ( - )  present;  ( - )  Ref.  2.  Experiment: 

($$$)  Ref.  22;  (+++)  Ref.  20. 
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CROSS  SECTION  (10**(-18)cm**2/sr) 


Fig.  6.  d<7oi/dft  at  1.51 
($$$)  Ref.  22; 


e"+N2  :  v(0— >1) 
DIFFERENTIAL  C/S 
2.05,  2.10,  2.15  eV 


Fig.  7.  dffoi/dH:  Theory:  ( - )  present  at  2.05  eV;  ( - )  present  at  2.10  eV; 

( - )  present  at  2.15  eV;  ( - )  Ref.  2.  Experiment:  ($<£$)  Ref.  22;  (ooo)  Ref. 

23. 
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e"+N2 :  v(0->l) 

DIFFERENTIAL  C/S 
3.00  eV 


ANGLE  (deg) 


Fig.  8.  d<7oi/dft  at  3.00  eV.  Theory:  ( - )  present;  ( - )  Ref.  2;  ( - )  Ref.  24. 

Experiment:  ($$4*)  Ref.  22;  (ooo)  Ref.  23;  (•••)  Ref.  25. 
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